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CHAPTER  1 
INTRODUCTION 


1.1  Overview 


PATCHEf.-III  is  a  computer  program  for  the  analysis  of  general 
three-dimensional  composite  structures,  laminated  composite  structures, 
and  axisynmetric  composites.  It  determines  the  linear  elastic  response 
of  a  general  heterogeneous  anisotropic  solid  under  thermal  loads, 
mechanical  loads,  and  imposed  displacements.  The  analysis  model  is 
based  on  a  64-node  isoparametric  solid  finite  element  with  variable 
material  properties.  Itiis  element  efficiently  models  the  strain 
discontinuities  at  heterogeneous  material  boundaries,  as  well  as  the 
continuous  strains  at  homogeneous  boundaries.  The  program  constructs 
models  for  geometry  and  pnysical  data  using  piecewise  parametric  cubic 
modeling  of  lines,  surfaces,  and  volumes.  Ibis  approach  eliminates  the 
grid  point  modeling  restriction  of  most  finite  element  programs  by 
allowing  the  synthesis  of  continuous  models  for  both  geometry  and 
physical  data.  It  is  necessary  to  input  coordinates  for  at  most  a  few 
grid  points;  coordinates  for  all  other  grid  points  are  computed 
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internally,  as  are  the  parametric  coefficients  in  the  mathematical 
models  created  for  lines,  surfaces,  and  volumes.  There  are  over  two 
dozen  bulk  data  directives  available  to  create  the  geometry  in  addition 
to  the  basic  grid  point  input.  These  bulk  data  directives  may  cross 
reference  each  other  and  may  be  input  in  any  order.  PATCHES-III  will 
process  them  in  the  proper  sequence  and  provide  diagnostics  if  there  are 
unresolvable  ambiguities  in  the  data.  A  similar  set  of  options  is 
available  to  model  temperature  data  or  any  physical  parameter  defined 
over  a  line,  surface,  or  volume. 

The  PATCHES-III  modeling  system  for  geometry  is  the  same  as  that 
used  in  the  PDA/PATRAN-G  interactive  graphics  system.  Models 
constructed  using  PATRAN  (Phase  1)  can  be  directly  input  to  PATCHES. 
There  are  minor  syntax  differences,  but  at  the  data  level  a  patch  is  a 
patch  is  a  patch.  A  translator  program  is  available  to  convert  the 
Phase-1  neutral  file  output  by  PATRA.N  into  LINE,  PATCH,  and  HPAT  data 
for  PATCHES.  The  Case  Control  and  Bulk  Data  syntax  of  PATCHES  are  the 
same  as  those  for  NASTRAN  to  facilitate  use  of  the  program  by  the  many 
people  experienced  with  data  input  to  NASTRAN.  There  are  major 
differences  in  the  two  finite  element  modeling  systems,  but  the  input 
syntax  is  very  similar. 

The  User's  Manual  is  the  source  book  or  dictionary  for  all  the  Case 
Control  and  Bulk  Data  directives  available  in  the  program.  A  limited 
amount  of  background  material  on  how  to  model  with  parametric  cubics  and 
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a  few  demonstr.it' ion  problems  ore  provided.  Those  interested  in  n  more 
comprehensive  description  of  parametric  cubic  modeling  will  find 
References  1,  2,  and  3  helpful,  and  the  PDA/PATRAN-G  User's  Guide  has 

many  interesting  examples.  The  documentation  for  specific  application 
areas  such  as  laminated  nozzle  components  is  available  as  separate 
User's  Guides.  The  demonstration  problems  in  this  manual  illustrate  the 
stress  precision  of  the  finite  element  library,  xfaich  is  unusually  high. 
A  graphic  example  of  this,  Figure  1-1,  shows  a  comparison  of  a 
31-element  solution  with  photoelastic  test  results.  The  maximum 
computed  stress  concentration  factor  was  equal  to  the  experimental 
result  to  three  significant  digits. 


1.2  Parametric  Cubic  Modeling 

Models  based  on  parametric  cubic  polynomials  are  developed  in 
detail  in  Reference  1.  This  chapter  introduces  the  user  to  the 
t  iminology  and  a  few  basic  concepts  from  that  development.  Finite 
dimension  lines,  surfaces,  and  volutes  are  continuous  geometry  modeling 
elements  added  to  the  familiar  grid  point.  The  mappings  Z(£), 
Z(c !»  ,  and  Z(£lf  £3)  that  define  these  continuous  elements  are 

cubic,  bicubic,  and  tricubic,  respectively,  based  on  the  four  Hermite 
polynomials,  FjU),  often  called  the  beam  functions.  They  map  the  unit 
interval,  unit  square,  and  unit  cube,  respectively,  into  any  desired 
shape  (see  Figure  1-2).  These  new  modeling  elements  are  mathematically 
related  to  each  other  in  a  very  special  way  that  allows  a  more 
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complicated  element  to  be  uniquely  constructed  from  the  union  of  simpler 
elements  and  in  a  variety  of  constructions.  A  hyperpatch,  for  example, 
can  be  defined  uniquely  by  64  points,  16  lines,  or  4  patches,  'ftiis 
property  makes  possible  the  model  construction  procedures  described  in 
Chapter  2  of  this  manual. 

The  points,  lines,  patches,  and  hyperpatches  of  the  geometry  model 
are  given  by  their  components  with  respect  to  a  global  Cartesian  frame 
whose  basis  vectors  are  denoted  by  ej,  e2,  and  e^  (see  Figures  1-3  and 
1-4).  The  tangent  vectors  Z,^  define  the  vectors  a.  =  Z and  they 
are  used  to  construct  local  curvilinear  (parametric)  frames.  Surface 
constraints  and  loads  are  usually  given  by  their  components  with  respect 
to  a  local  parametric  frame,  as  Chapter  2  describes  in  detail.  Tbe  8 
corner  points  of  a  hyperpatch  are  called  grid  points  in  PATCHES-Iir  and 
are  used  to  define  the  connectivity  of  a  finite  element.  The  64  points 
associated  with  a  uniform  parametric  mesh  (i.e.,  the  1/3  points)  are 
termed  mesh  points,  and,  obviously,  8  of  these  points  are  also  grid 
points.  Mesh  points  normally  are  not  used  to  construct  models,  but  they 
are  the  recovery  points  for  all  element  data  in  PATCHES-III.  The 
program  automatically  determines  mesh  point  connectivity  from  the  grid 
point  connectivity. 
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Cigure  1-3.  PATCHES-III  Tangent  Vectors  for  Lines  and  Surfaces 
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CHAPTER  2 


PATCHES-I II  CAPABILITIES 


2.1  Overview 


PATCHES-I I I  will  compute  for  a  general  solid:  3  components  of 
displacement,  6  components  of  strain,  and  6  components  of  stress  at  64 
points  in  each  finite  element  used  to  model  the  body.  Hie  basic 
difference  in  capability  with  respect  to  other  three-dimensional  finite 
element  programs  is  the  construction  of  parametric  cubic  models  for  both 
geometry  and  physical  data  by  the  program.  This  allows  the  finite 
element  model  to  be  created  with  very  little  input,  and  the  model  itself 
efficiently  predicts  the  response  of  heterogeneous  structures.  These 
characteristics  are  essential  to  the  solution  of  monolithic 
three-dimensional  problems  with  the  resources  normally  available  for  a 
structural  analysis.  This  section  describes  techniques  for  creating 
parametric  cubic  models  and  the  matrix  solution  method  available  in 
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In  mathematical  terms  a  parametric  cubic  model  is  one  in  which  the 
interpolation  functions  are  multivariate  and  piecewise  cubic  in  the 
parameters,  Those  interested  in  a  mathematical  development  are 
referred  to  Reference  1.  This  discussion  is  for  the  structural  analyst 
familiar  with  finite  element  modeling  and  interested  in  using 
PATCHES-III  without  becoming  expert  in  parametric  cubics.  The  basic 
concept  is  construction:  lines  from  points,  surfaces  from  lines,  and 
hypersurfaces  from  surfaces.  Parametric  models:  Z  =  (Z1(£)), 
Z  =  (Z1(CJ,C?)),  and  Z  =  (Z  1(  t  ,  e;^,  £, J )  are  used  because  they  are 
convenient  for  numerical  work  and  independent  of  the  shape  of  the  line, 
surface,  or  hypersurface.  A  geometric  line  requires  12  parameters  for 
definition,  a  geometric  surface  patch  48  parameters,  and  a  geometric 
hypersurface  patch  (hyperpatch)  192  parameters.  These  parameters  may  be 


uniquely  given  in  3  different  but  mathematically  equivalent  ways:  (1) 
Algebraic,  in  which  the  coefficients  of  or  or  ate  given; 


(2)  Geometric,  in  which  the  derivatives  Z,_  or  Z, .  or  Z,  t  r  are 

-  si  -  sitj  - 

given  at  the  corners;  and  (3)  Point,  in  which  the  value  of  Z  is 


specified  at  4,  16,  or  64  points  and  where  certain  of  these  points  must 


be  interior  points.  How  the  user  directs  PATCHES-III  to  compute  these 


parameters  will  be  described  briefly. 
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2.2  Geometry  Modeling 

Construction  of  the  geometry  model  is  accomplished  through  the 
translation,  rotation,  segmentation,  interpolation,  and  scaling  of  space 
figures.  These  figures  may  be  points,  lines,  surfaces,  or 
hypersurfaces,  any  of  which  also  may  be  the  result  of  a  previous 
operation.  This  latter  feature  allows  the  synthesis  of  complex  models 
from  simple  input.  To  avoid  order  dependent  input,  PATCHES-III  contains 
a  queuing  algorithm  that  determines  the  processing  order  required  to 
synthesize  the  geometry.  Diagnostics  are  provided  when  bulk  data 
directives  produce  conflicting  data,  such  as  different  coordinates  at  a 
grid  point,  or  when  bulk  data  directives  contain  unresolvable 
ambiguities,  such  as  a  reference  to  coordinates  never  computed. 

Consider  the  construction  of  a  hyperpatch  for  one  segment  of  a 
thick-walled  circular  cylinder.  Two  quite  different  but  virtually  equal 
constructions  illustrate  the  procedure  for  a  simple  shape.  In  the 
first,  grid  points  1  and  3  are  input  (Figure  2-1),  and  a  straight  line 
connecting  them  is  created  with  a  LINEPC  card.  This  line  is  rotated 
about  the  e^  axis  through  90  degrees  to  form  one  quadrant  of  a 
cylindrical  surface  using  a  PATCHR  directive.  In  this  process,  grid 
points  7  and  5  are  automatically  created.  The  surface  1-3-7-5  is 
expanded  a  unit  amount  in  the  direction  of  its  normal  to  create  a 
hyperpatch  using  the  HPN  directive.  The  grid  points  2,  4,  6,  and  8  are 
automatically  created  in  this  process  along  with  all  the  hyperpatch 
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parameters.  The  new  grid  point  identification  numbers  are  determined  by 
the  element  connectivity  specified  with  a  CPDE3  directive.  Hie 
construction  of  this  192-parameter  hyperpatch  required  five  (5) 
directives  of  very  simple  format.  Now  consider  the  same  construction 
problem,  but  this  time  input  grid  points  1,  2,  3,  and  4  (Figure  2-2).  A 
quadrilateral  surface  is  created  with  a  PATCHQ  directive,  and  this 
surface  is  rotated  about  the  e^  axis  through  90  degrees  to  form  a 
hyperpatch  using  the  HPR  directive.  The  construction  this  time  required 
six  (6)  bulk  data  directives.  Although  the  two  figures  constructed  are 
nearly  equal  in  volume  and  shape,  the  two  hyperpatches  are  quite 
different  parameter i2at ions  of  the  same  figure.  In  the  case  of  method 
one,  the  surface  Zl (r  ,0)  is  1 -3-7-5,  but  for  method  two  the  surface 
Z  (f.  ,f,  ,0)  is  1-5-6-2.  The  order  of  the  parameters  in  a  hyperpatch 
is  the  parameterization  or  "sort*  of  the  hyperpatch.  Fortunately,  a 
hyperpatch  is  invariate  under  reparameterization  in  PATTHES-III  because 
multivariate  interpolation  functions,  F.(f  )F  (c,  )F  (r  ) ,  are  used.  The 

3  1  J  ?  k  3 

program  allows  each  finite  element  to  be  in  any  right-handed 
parameterization;  however,  all  data  input  or  data  created  for  an 
element  must  use  the  same  parameterization.  This  subject  will  be 
discussed  further  in  the  data  modeling  section. 
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The  two  constructions  just  presented  for  the  sample  figure  are  but 
two  of  many  that  could  have  been  made.  Experience  with  parametric  cubic 
modeling  indicates  the  average  user  quickly  develops  his  own  techniques 
for  creating  geometry.  As  a  rule  of  thumb,  simple  figures  (cubes,  etc.) 
should  be  used  for  finite  elements  vrfienever  possible.  Elements 
completely  interior  to  a  body,  for  example,  can  be  parallelepipeds 
without  changing  tne  surface  geometry.  This  reduces  the  cost  of 
integrating  the  stiffness  matrix  for  the  associated  finite  elements, 
sometimes  by  a  factor  of  five. 
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BULK  DATA  DIRECTIVES 


GRID,  1, ,  1.0 
GRID,  3,,  1.0,  1.0 
LINEPC,  1,  1,3 

PATCHR.20,  1,,,,  0.0,  90.0,3 
HPN,  1,  20,  1.0 


CPDE3,  1,  1,  3,7,5 . 2,  4,8,6 


STEP 

BULK  DATA  OPTION 

INPUT 

1 

GRID 
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PATCHR 
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HPN 

1 
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Figure  2-1.  Hyperpatch  Construction  -  Method  One  Example 
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Before  leaving  geometric  modeling,  let  us  examine  the  components  of 
a  space  curve  represented  in  parametric  form  to  help  describe  the  basic 
concept  behind  PATCHES  geometric  models.  It  may  be  helpful  at  first  to 
think  of  the  parametric  parameter  E  as  representing  time  or  arc  length 
or  any  parameter  that  increases  monotonically  as  we  progress  from  one 
point  on  the  space  curve  Z(t)  to  another,  say,  from  time  t  =  0  to  t  =  T. 
Each  coordinate  function,  Z1 (t) ,  is  a  well-behaved  function  on  the 
interval  0  <  t  <  T,  as  illustrated  by  the  helix,  Z(t)  =  (R  cos  t,  R  sin 
t,  t) ,  Figure  2-3.  Using  a  monotonic  parameter  like  t  avoids  all  the 
tangent  vector  singularity  problems  (removable)  associated  with  line 
equations  in  traditional  analytic  geometry.  The  tangent  vector  is 
simply  Z,.  =  (-R  sin  t,  R  cos  t,  1),  and  all  the  geometric  properties  of 
the  curve,  such  as  curvature  and  twist,  can  be  computed  easily. 

The  last  step  is  to  represent  each  coordinate  function,  Z 1  (t) ,  with 
a  cubic.  First  we  change  to  a  normalized  parameter,  0  <  £  <  1,  where 
t  =  T£.  Then  PATCHES  computes  PC  coefficients  such  that  each  cubic 
interpolates  each  coordinate  function.  The  algebraic  coefficients  for 
line  coordinate  functions,  S^,  in  the  case  of  a  helix  are 

R  cos  EX  -  S^3  +  SgC2  +  +  sj  (2.1.1) 

R  sin  EX  2  S^C3  +  S2S2  +  S3 E  +  S2  (2.1.2) 

ft  =  0-e3  +  0'52  +  s\e  +  0  (2.1.3) 
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Figure  2-3.  Space  Curve  Defined  in  Parametric  Form 
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hyperpatch  associated  with  the  referenced  finite  element.  This  is  the 
connection  between  data  and  geometry  that  makes  the  parameterization 
(ordering)  of  the  hyperpatch  important.  As  a  user  convenience,  there  is 
an  option  that  will  reorder  any  data  hyperpatch,  but  note  that  since  the 
data-geometry  relation  is  implicit,  any  ordering  of  the  input 
coefficients  is  correct.  Only  the  user  can  determine  if  the  input  data 
distribution  is  oriented  as  he/she  intended.  A  DRY  run  is  recommended 
if  there  is  any  doubt  about  input  data  distributions.  This  allows  the 
user  to  inspect  the  data  models  in  point  format.  Another  user 
convenience  is  the  data  patch  equivalence  option.  This  will  allow  the 
same  data  patch  to  be  used  on  many  different  surfaces,  and  it  is 
particularly  useful  for  constant  pressure  loads.  Note  that  since  the 
data  are  modeled  in  parametric  space,  a  constant  pressure  data  patch  is 
the  same  for  every  surface.  It  merely  gives  the  magnitude  of  the 
pressure;  the  direction  of  the  surface  normal  is  determined  from  the 
hyperpatch  for  the  finite  element  loaded  by  the  pressure.  There  are 
options,  such  as  the  DPATCH  directive,  provided  especially  for  direct 
input  data  modeling.  This  particular  option  could  be  used  for  a  problem 
in  which  the  temperature  is  a  quadratic  function  of  radius.  In  such 
cases,  it  is  more  convenient  to  model  the  data  directly,  rather  than  to 
synthesize  the  model  from  grid  point  values.  In  other  cases,  it  may  not 
even  be  possible  to  synthesize  the  model  from  grid  point  values.  As  a 
convenience  in  these  situations,  options  are  provided  that  are  divorced 
from  any  grid  point  data  sets. 
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PERCENT  ERROR  IN  THE  RflOIUS 
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NORMALIZED  PERCENT  ERROR  FOR  90  DEG  PC  ARC 
AS  A  FUNCTION  OF  ANGLE 

ERROR  =  (COMPUTED  RADIUS  -  REAL  ARC  RADIUS)»100-/RERL  RADIUS 


Figure  2-5.  LARCPC  Radial  Error  Distribution 
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In  addition  to  distributed  data,  the  user  can  input  temporal  data 
using  data  table  functions.  The  functions  are  again  piecewise 
parametric  cubics  in  one  variable  at  present.  Temperature-dependent 
material  properties  are  input  using  data  table  functions,  and  in  advance 
development  versions  of  the  program  properties  can  be  functions  of  other 
state  variables.  Data  table  functions  are  nondimensional  in  the  sense 
that  the  data  assume  physical  significance  only  when  they  are  referenced 
and  used  by  some  other  directive. 


2.4  Constraint  Modeling 

Most  imposed  displacement  boundary  conditions  are  specified  over  a 
surface,  not  simply  at  a  grid  point.  As  a  user  convenience,  both  zero 
and  nonzero  constraint  options  are  available  that  constrain  the  entire 
face  of  an  element  with  a  single  bulk  data  directive.  Either  the 

reference  frame,  e.,  or  the  local  surface  frame,  a  ,  a  ,  n,  (Figure  2-6) 

„  1  -  1  -  2  - 

rr  y  be  used.  In  the  case  of  zero  constraints,  only  the  corner  grid 
point  numbers,  frame  identification  number,  and  displacement  component 
numbers  are  input.  In  the  case  of  nonzero  constraints,  up  to  3  data 
patches  per  surface  may  be  specified  that  define  the  magnitude  of 
constrained  displacement  components  over  the  surface.  In  both  cases, 
PATCHES-I I I  computes  any  required  frame  transformation  for  all  16  points 
on  the  surface. 
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At  the  intersection  of  constrained  surfaces,  a  number  of  abstruse 
conditions  can  occur.  The  constrained  displacement  components  may  be  in 
different  frames  at  the  same  point;  the  constraints  may  be  redundant, 
and  when  the  user  makes  an  error  they  may  be  inconsistent.  The  program 
allows  multiple  frames  at  a  point  and  redundant  constraints,  and  it  will 
provide  diagnostics  when  the  user  specifies  inconsistent  constraints.  A 
general  vector  approach  is  used  that  synthesizes  a  local  frame  at  every 
constrained  point  from  the  linearly  independent  constraint  vectors.  The 
details  are  described  in  Reference  4.  It  is  possible  using  a  debug 
option  (PARAM  directive)  to  output  the  transformation  matrix  relating 
vector  components  in  the  e^  frame  to  components  in  the  local  constraint 
frame. 
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2.5  Constraint  Finite  Elements 


The  finite  elements  used  to  model  laminate  force-deformation 
behavior  are  a  family  of  linear  constraint  options  developed  for  Version 
9.0  and  higher  of  the  PATCHES-III  program.  The  need  for  low-cost 
modeling  in  regions  of  uniaxial  or  biaxial  strain  was  noted  in  earlier 
versions,  and  the  new  elements  in  Table  2-1  provide  this  capability. 
They  are  based  on  the  linear  constraints  defined  by 


P  .  =  T,  P 

1  1U  <1 


i  =  1,  2,  3, 
a  =  1,  4 


where  the  coefficients  Tia  are  simply 


T.  = 

ia 


1 


0 


2/3  ,  1/3 
1/3  ,  2/3 
0  ,  1 


(2.5-]) 


(2.5-2) 


The  same  coefficients  apply  to  all  three  parametric  coordinates  and,  in 
general , 


P  =  T  T  T  P 

ijk  in  ,16  k>  aBy 


(2.5-3) 


If  constraints  are  introduced  in  only  two  coordinates, 

P  =  T  T  5  P 

ijk  ia  JB  k2  a62 

and  for  only  one  constraint 

P  =  T  <S  6  P 

ijk  ia  J  2  km  u2.ni 


(2.5-4) 


(2.5-5) 
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Table  2-1 

PATCHES-III  CONSTRAINT  FINITE  ELEMENT 


Displacements* 

Nodes 

Geometry 

Properties 

LLL 

8 

CCC 

CCC 

LLC 

16 

CCC 

CCC 

LCC 

32 

CCC 

CCC 

*  Any  combination  of  L  and  C  is  available.  L  =  linear,  C  =  cubic. 
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Two  key  issues  affecting  the  development  of  the  new  family  of 
elements  are  how  to  efficiently  generate  their  stiffness  matrices  and 
how  to  connect  them  to  each  other.  After  some  early  confusion,  it  was 
determined  that  all  linear  constraints  can  be  applied  before  integration 
with  the  same  result  as  when  they  are  applied  after  integration.  This 
greatly  reduces  the  cost  of  generating  their  stiffness  matrices.  It  is 
simply  to  demonstrate  this  equivalence  in  one  dimension  vfriere, 
obviously, 

KaB  *  /Ti„  fVfJ  FJ<«  TJ8dt 
-Ti„  FJ<«  «TJ« 

-Ti«Ki3TJB  l2-5-6' 

But  in  higher  dimensions  interpolatory  quadrature  is  used  in 

PATCHES-III,  and  this  equivalence  is  tedious  to  prove. 

The  second  issue  was  resolved  by  automating  the  generation  of 
interface  constraints  between  elements  ot  different  dimension.  This 
allows  the  user  of  PATCHES-III  to  specify  linear  constraints  on  any 
element  or  group  of  elements  by  simply  placing  a  mnemonic  of  the  type 
listed  in  Table  2-1  on  the  connectivity  card  for  that  element.  The 
program  first  generates  all  explicit  mesh  point  constraints  and  then  on 
a  second  pass  generates  all  interface  or  implicit  constraints.  This 
requires  extensive  checking  for  conflicts,  and,  in  order  to  reduce  their 
incidence,  all  three  displacement  components  are  constrained  alike.  At 
every  constrained  mesh  point,  one  of  the  Equations  (2.5-3)  -  (2.5-4)  is 
automatically  generated  and  applied  to  all  affected  matrices. 
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A  family  of  finite  elements  based  on  axisynmetric  constraints  is 
also  available.  These  elements  require  hyperpatches  of  the  form 

zV,,  f.  =  r(f, r  ?)  sin  ^ 

7.  '(Lj,  f,  ,„Cj)  =  r  (tn,  f,  2)  cos  c3 
Z3(s1,  5  ?,c3)  =  Z(51,  e  2) 

which  can  be  generated  using  the  HPR  directive.  Note  the  convention 
here  that  associates  the  hoop  coordinate  8  with  the  third  parametric 
coordinate  Ky  A  generalized  axisynmetric  displacement  constraint  that 
allows  torsion  U1,^  =  0  results  in  bicubic  displacement  functions  for 
Up,  Up,  and  Uz.  The  axisynmetric  finite  elements  are  designated  CCX, 
etc.,  as  listed  in  Table  2-2,  and  they  reduce  the  dimension  of  an 
element  by  a  factor  of  4.  In  general,  the  number  of  varying  strain 
components  remains  6  because  of  the  torsional  response  mode.  This 
behavior  is  typical  of  axisynmetric  composites  of  involute  construction. 
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Table  2-2 

PATCHES-III  AXISYEMETRIC  FINITE  ELEMENTS 


Displacements* 

Nodes 

Geometry 

Properties 

LLX 

4 

CCC 

CCC 

LCX 

8 

CCC 

CCC 

CCX 

16 

CCC 

CCC 

*  Any  combination  of  L  and  C  in  the  first  two  positions  is  available 
with  X  in  the  third  position.  L  =  linear,  C  =  cubic,  X  =  axi symmetric. 
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2.6  Solution  Method 

PATCHES-III  utilizes  a  scaled  conjugate  gradient  solution  procedure 
which  has  been  shown  to  be  extremely  effective  for  three-dimensional 
problems.  Core  requirements  increase  very  slowly  with  increasing 
problem  size,  virtually  eliminating  the  "spill"  problem  associated  with 
direct  decomposition  on  many  computers.  This  is  primarily  the  result  of 
only  a  single  element  stiffness  matrix  being  in  core  at  any  one  time. 
Also,  since  the  global  stiffness  matrix  is  never  constructed,  bandwidth 
minimization  considerations  are  eliminated. 

The  solution  procedure  has  been  formulated  around  the  dot  product 
operation,  which  permits  the  program  to  take  full  advantage  of  recent 
developments  in  machine  architecture.  As  a  result,  the  cost  per 
iteration  is  extremely  low,  and  the  total  cost  to  convergence  is 
generally  substantially  lower  than  for  direct  decomposition.  For  most 
isotropic  problems,  engineering  accuracy,  if  not  convergence,  will  be 
obtained  in  approximately  N/4  iterations  where  N  is  the  number  of 
degrees  of  freedom,  called  NFSET.  In  orthotropic  problems,  the 
convergence  is  slower,  depending  on  the  E/G  ratio  of  the  material.  Most 
laminate  problems  converge  in  N  cycles,  but  laminates  with  rubber  plies 
mixed  with  high  modulus  plies  can  take  2N  cycles.  Models  requiring  more 
than  2N  cycles  are  very  ill  conditioned  and  often  indicate  a  modeling 
error  of  some  sort.  Figure  2-7  shows  slow  convergence  caused  by  two 
rubber  plies  in  an  aircraft  window  model  that  act  as  shear  strain 
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isolators.  There  is  a  checkpoint  feature  provided  to  allow  continued 
iterations  on  a  subsequent  restart  run  if  additional  iterations  are 
required.  Through  default  values  or  specific  overrides,  the  user  can 
control  the  iteration  processes  defining  convergence,  maximum  cycles, 
maximum  time,  cut-off  for  rigid  body  modes,  checkpoint,  restart,  and 
more.  Full  strain,  stress,  and  force  recovery  is  always  available. 
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CHAPTER  3 

A  GUIDE  TO  THE  USER  MANUAL 


3.1  Overview 

After  reading  the  introduction  and  Chapter  2,  the  user  should  have 
a  general  understanding  of  the  modeling  features  available  in 
PATCHES-III.  This  chapter  introduces  the  user  to  the  actual  operation 
of  the  program  and  describes  its  architecture.  Those  familiar  with 
NASTRAN  will  recognize  the  bulk  data  and  case  control  input  formats. 
The  syntax  is  virtually  the  same,  including  order  independent  bulk  data 
and  checkpoint-restart  processing  which  has  been  simplified  to  eliminate 
the  checkpoint  dictionary.  Chapter  4  provides  a  detailed  description  of 
the  bulk  data  directives  in  the  familiar  NASTRAN  format.  Chapter  5 
provides  the  same  information  for  the  case  control  data  directives,  and 
Chapter  6  describes  the  optional  executive  requests.  Chapter  7  shows 
the  input  and  output  for  several  sample  problems  which  illustrate  the 
basic  simplicity  of  modeling  with  parametric  cubics.  Chapter  8 
illustrates  job  control  options  (tape  request,  etc.),  and  Chapter  9 
describes  diagnostic  and  user  information  data  output  by  the  program. 
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3.2  Data  Preparation 

The  sequence  of  input  directives  necessary  for  execution  of 
PATCHES-III  is  similar  to  that  necessary  for  execution  of  NASTRAN.  The 
input  file  is  shown  schematically  in  Figure  3-1  for  a  single  run.  It  is 
possible  to  concatenate  files,  but  this  would  not  normally  be  done. 
Detailed  information  on  the  preparation  of  executive,  case  control,  and 
hulk  data  may  be  found  in  Chapters  4,  5,  and  6  of  this  manual. 
Ii  formation  on  computer  processing  of  the  deck  and  a  schematic  of  the 
PATCHES-III  load  map  are  presented  next. 
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Figure  3-1.  PATCHES-III  Input  File 
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3 . 3  Computer  Processing 

To  use  PATCHES -i I !  efficiently  it  is  necessary  to  understand  the 
basic  structure  of  the  program  and  the  relative  expense  of  individual 
modules.  This  allows  the  user  to  plan  the  dry  runs,  checkpoints,  and 
restarts  that  best  suit  his  application.  A  schematic  of  the  PATCHES-III 
load  map  is  shown  in  Figure  3-2  where  the  core  storage  each  link 
requires  is  indicated  by  the  length  of  the  line  for  that  link.  The 
program  is  loaded  on  CDC  computers  using  the  segmentation  loader.  The 
resident  link,  Link  0,  is  always  in  core.  It  contains  the  PATCHES-III 
executive  system  and  the  communication  data  blocks.  Each  of  the  other 
links  is  executed  as  needed  under  the  control  of  the  PATCHES-III 
executive.  A  brief  description  of  the  function  of  each  link  is  provided 
in  Table  3-1. 

At  the  end  of  execution,  a  user  information  table  is  printed  which 
shown  the  CP  time,  core  storage,  and  random  access  disk  requirements  of 
each  link  and  major  subregions  within  complex  links.  The  various 
PATCHEc-lII  links  do  not  require  a  fixed  amount  of  core  storage  but, 
rather,  an  "open  core"  concept  wherein  each  link  determines  and  assigns 
the  necessary  core  based  upon  the  problem  requirements.  This  allows 
PATCHES-III  to  process  very  large  problems  without  penalizing  those 
small  problems  that  can  be  solved  with  modest  storage  requirements 
(approximately  170000^).  To  check  all  input  data  and  the  geometry  model 
before  execution  on  large  problems,  use  the  Case  Control  option  DRY. 
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This  will  terminate  the  run  at  the  flag  labeled  DRY  in  Figure  3-2.  Less 
than  10  percent  of  the  cost  of  a  complete  solution  is  required  for  a  DRY 
run,  and  the  percentage  decreases  with  increasing  problem  size. 

There  are  two  volumes  of  data  that  may  be  saved  on  two  separate 
files  of  a  checkpoint  tape  for  later  use  on  a  restart  run.  The  first 
volume  contains  the  element  stiffness  matrices,  and  the  second  contains 
the  current  iterate  and  direction  vector  from  the  iterative  solution. 
One  or  both  volumes  may  be  created  on  a  checkpoint  run  using  the  case 
control  options  CHKPNT,  ELEMENT  or  CHKPNT,  CG  or  CHKPNT,  ELEMENT,  CG. 
On  a  subsequent  restart  run,  any  or  all  of  the  element  matrices  can  be 
modified  using  the  option  RESTART,  ELEMENT,  LIST.  All  restart  runs 
require  the  complete  case  control  and  bulk  data  decks  to  be  present, 
even  on  a  CHKPNT,  CG  run.  This  is  necessary  because  all  data  not  on  the 
checkpoint  tape  (geometry  data,  etc.)  must  be  regenerated.  A  restart 
run  can  create  an  additional  GG  checkpoint  file,  thereby  allowing  the 
user  to  monitor  intermediate  results  in  a  large  analysis. 
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Table  3-] 


PATCHES-I I I  LINK  DESCRIPTIONS 


Link 

Region  Name 

Function 

0 

MAIN 

Executive  control  and  common  storage 

1 

BEGIN 

Initialize  PATCHES-III  system 

2 

INPCN 

Input  control  -  construct  geometry  and  data 
models 

3 

GET 

Initialize  integration  tables 

4 

MATCN 

Material  property  card 

13 

LOADS 

Generate  element  load  vectors 

7 

BIGMSH 

Generate  element  mesh  point  connectivity 

8 

IMDISP 

Imposed  displacement  model 

9 

EM  PC 

Element  mesh  point  constraints 

5 

*EKIJ 

Generate  element  stiffness  matrices 

10 

TTEKTD 

Transform  element  matrices  to  analysis 
coordinates 

13 

EXSCE1 

Single  point  constraint  eliminator 

14 

MCE1 

Multipoint  constraint  eliminator 

16 

*SSG 

Static  solution  generator 

17 

SDR 

Stress  recovery 

11 

SUBCOM 

Subcase  combinations 

Links  that  use  majority  of  the  CP  time 
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BULK  DATA  OPTIONS 


4.1  Overview 

The  primary  source  of  input  to  PATCHES-III  is  a  bulk  data  file 
whose  syntax  is  the  same  as  that  used  by  NASTRAN,  Reference  2.  This 
data  defines  the  geometry,  physical  properties,  boundary  conditions,  and 
loading  conditions  for  the  finite  element  model  of  the  structure. 
Certain  execution  parameters,  such  as  grid  point  tolerances,  can  also  be 
input  with  bulk  data  directives,  all  of  which  may  be  submitted  in  any 
order.  PATCHES-III  preprocesses  the  bulk  data  file  to  determine  the 
order  in  which  directives  must  be  processed  to  account  for  data 
hierarchies.  The  input  data  may  be  in  one  of  two  formats:  free-form  or 
fixed-form.  Fixed-form  is  identical  in  format  to  standard  NASTRAN  bulk 
data  cards.  Ihe  preferred  format  Cor  PATCHES  Control  Case  and  Bulk  Data 
Directives  is  free-form  because  of  the  string  notation  syntax  which  is 
convenient  for  input  via  a  CRT  terminal. 
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4.1.1  Free-Form  Input 

1.  Definition:  A  directive  is  assumed  to  be  of  free  form  if  a 
comma  exists  anywhere  on  the  first  physical  card  (line)  of  that 
directive.  Continuation  cards  (lines)  are  assumed  to  be  of  the 
same  format.  Colimn  positions  and  all  blanks  are  ignored  by 
the  free-form  processor. 

2.  Syntax:  The  data  input  to  PATCHES  is  broken  into  individual 
fields.  Each  rield  represents  a  particular  type  of  data.  A 
field  may  consist  of  any  number  of  characters  and  digits  and  is 
terminated  by  a  comma.  As  an  example,  consider  the  following 
directive: 

HPR,  7,  10,,,,  0.0,  90.0,  3 

As  in  a  NASTRAN  bulk  data  card,  the  first  field  is  the 
name  of  the  directive,  in  this  case,  HPR.  The  second  field 
identifies  by  number  the  item  to  be  created,  in  this  case, 
HYPERPATCH  number  7.  The  third  field,  10,  identifies  the  patch 
to  be  used  in  the  construction  option.  The  next  three  fields 
are  nul  for  this  option,  indicated  by  blanks  between  the 
commas.  This  is  similar  to  DMAP's  syntax.  The  last  three 
fields  on  this  directive  define  rotation  angles  and  rotation 
axis.  The  primary  rule  to  note  is  that  a  comma  terminates  a 
field;  column  positions  are  of  no  importance. 
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After  the  Linal  field  in  the  above  example,  a  comma  is  rx>t 
required.  If  the  last  nonblank  character  on  a  card  (line)  is  a 
comma,  the  next  card  (line)  is  assumed  to  be  a  continuation 
card  (line) .  Data  continues  on  the  next  card  as  if  the 
pnysical  card  were  in  excess  of  80  columns.  A  PATCHES  data 
directive  can  have  up  to  eight  continuation  cards,  but  an 
individual  field  cannot  cross  card  boundaries. 

3.  In-Line  Lists:  There  are  circumstances  in  which  a  list  can 
occur  within  a  field.  Such  a  field  can  be  input  using  a  slash 
,/,  in  place  of  the  comma  to  separate  items.  An  example  would 
be 

PATCHO,  11,  5/3,  4 

As  before,  a  field  cannot  extend  over  card  boundaries. 


4.1.2  Fixed-Form  Input 

1.  Definition:  A  directive  is  assumed  to  be  of  fixed  form  if 
there  are  no  commas  on  the  first  physical  card  (line)  of  that 
directive.  Continuation  cards  are  assumed  to  be  of  the  same 
format. 

2.  Syntax:  The  format  of  the  directive  is  identical  to  that  of  a 
NASTRAN  bulk  data  card.  Fields  exist  totally  within  8-  or 
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16-column  blocks.  Input  is  similar  to  free-form  input  if  you 
imagine  a  cornua  between  blocks.  The  following  two  cards  result 
in  identical  definitions: 

LIST,  5,  8,  12,  16  THRU  20,  34,  4,  3,  9,  40 
LIST  f,  8  12  16THRU20  34  4  3  9  +CONT 

+C0NT  40 

where  eight  columns  must  occur  in  each  field  of  the 

fixed-format  cards  indicated  by  underlining. 

Restrictions  and  use  of  fixed  format  are  otherwise 
identical  to  those  of  free  format.  The  16-character  input 
field  format  is  specified  by  appending  an  asterisk,  *,  to  the 
mnemonic  in  field  1  or  by  placing  an  asterisk  in  column  1  of  a 
continuation  card. 
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4.2  Available  Options 

There  are  over  70  bulk  data  directives  available  to  model  th° 
structure,  its  properties,  and  its  environment..  These  are  grouped  into 
six  functional  categories  and  are  catalogued  by  mnemonic  in  Table  4-1 
for  easy  reference.  Few  of  these  directives  are  mandatory  on  any  given 
run.  PATCHES-III  analyzes  the  input  bulk  data  and  provides  diagnostics 
if  there  are  errors  or  omissions.  Every  attempt  is  made  to  analyse  the 
entire  bulk  data  set  independent  of  the  number  of  errors  that  may  be 
found.  'bn  most  runs  a  geometry  model  will  be  created  and  output,  even 
when  th ere  are  fatal  errors.  Cross-referencing  between  bulk  data 
directives  is  by  explicit  reference  to  an  identification  number.  The 
only  exception  is  the  CPDF3  directive,  which  implicitly  requires  the 
property  card  and  hyperpatch  card  to  have  the  same  identification  number 
as  the  finite  element.  The  number  of  elements  is  limited  to  512,  which 
cor r  -.-'ponds  to  a  range  of  stiffness  matrix  dimensions  from  5,000  to 
1 00,000,  depending  on  the  finite  element  type. 
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Table  4-1 


BULK  DATA  OPTIONS 


1 .  GEOMETRY* 

GRID 

LARCPC 

PATCH 

HPATCH  SCALP 

LINE 

PATCHGR 

HPHEX  SCALPH 

LINECS 

PATCHL 

HPL  WOVE 

LINEGR 

PATCHO 

HFN 

LINEPC 

PATCHQ 

HPR 

PATCHR 

HP2PAT 

PATCH4L 

HP4PAT 

HP6PAT 

2.  ELEMENT  AND 

PROPERTIES 

CPDE3 

MATAL 

MAT1 

MATTA 

PPDE3 

MATC 

MATTAT 

MATE 

MATTC 

MATOR 

MATTE 

MATTD 

MATT] 

3.  CONSTRAINTS 

MPE1 

SDC10 

SDC1 

SPC1 

MPE2 

SDC20 

SDC2 

SPC2 

4.  LOADS 

FORCE 

FORCEL3 

FORCET 

PLOAD3 

TEMP 

(Table  continued  on  following 

page.) 

*  Mnemonic  Suffixes: 

A  =  Algebraic 

,  B  =  Geometric,  CS  =  Cubic  spline,  HEX  =  Hexahedra, 

L  =  Line,  N  = 

Normal ,  P 

=  Point,  PC  =  Parametric  Cubic, 

Q  =  Quadrilateral,  R  =  Rotation. 
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Tablo  4-1  (Continued) 
BULK  DATA  OPTIONS 


5.  DATA  MODELING* 


DATAG 

DLINCS 

DPATA 

DHPAT 

DLINE 

DPATCH 

DHPHEX 

DLINP 

DP ATE Q 

DHPL 

DLINPC 

DPATL 

DHPSORT 

DPATQ 

DHP2P 

DPAT4L 

DHP4P 

6 .  MISCELLANEOUS 

PARAM  MTRX-ID 

$comment  MTRX-CID 


*  Mnemonic  Suffixes: 

A  =  Algebraic,  B  =  Geometric,  CS  =  Cubic  spline,  HEX  = 
L  =  Line,  N  =  Normal,  P  =  Point,  PC  =  Parametric  Cubic 
Q  =  Quadrilateral ,  R  =  Rotation. 


ETCS 
Err  PC 


Hexahedra, 
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4.3  Bulk  Data  Directives 

This  section  details  in  alphabetical  order  for  each  of  the  bulk 
data  directives  their  input,  format,  restrictions  and,  where  necessary, 
additional  information  concerning  the  use  of  the  particular  directive. 
The  descriptions  are  in  the  nature  of  dictionary  information  with 
illustrated  examples  for  key  directives.  Table  4-2  shows  the  format 
used  to  document  each  directive. 
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Table  4-2 

DOCUMENTATION  FORMAT  FOR  BULK  DATA 


CATEGORY 


CONTENT 


Input  Directive: 

Descr ipt ion : 

Format: 

Example  Syntax: 

Field: 

Remarks*: 

Commentary*: 

Example  Application*: 


Mnemonic  for  directive. 

Dictionary  description  of  directive. 

Mnemonic  for  each  input  field. 

Typical  input  data. 

Dictionary  description  of  each  input  field. 
Restrictions  and  assumptions. 

Background  information  for  certain  directives. 
Typical  application  is  illustrated. 


*  This  item  of  documentation  is  not  available  for  all  directives. 
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BULK  DATA  INPUT 


Input  Directive: 
Description: 


Format : 

Example  Syntax : 
Remarks: 


$  Comment 

A  comment  vhich  the  user  may  input  to  annotate 
the  bulk  data  file  in  its  unsorted  form.  TTiis 
card  is  ignored  by  the  bulk  data  processor. 

$  Any  legitimate  characters  in  colnnns  2-80. 

$****CRAfflITE  PHENOLIC  PROPERTIES. 

1.  Comment  cards  cannot  be  used  with  the  current 
Version  9.5  for  the  VAX-11,/780. 
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rnput  Directive: 
Description: 

Format : 

Example  Syntax: 

Field 

EID 

Gl,  G2,...r  G8 

EM  PC 

Remarks : 


BULK  DATA  INPUT 


GPDE3 

Connectivity  card  for  a  three-dimensional 
parametric  discrete  (finite)  element.. 

CPDE3,  EID,  Gl,  G2,  G3,  G4,,,,  EMPC,  G5,  G6, 
G7,  G8 

CPDE3,  7,  1,  3,  13,  12,,,,  OCL,  2,  4,  14,  5 


Contents 


Element  identi l ication  number. 

Corner  grid  point  identification  numbers  in  the 
sequence  given  by  Figure  4-1. 

Mesh  point  constraint  mnemonic  for  elements 
reduced  to  less  than  cubic  in  one  or  more  of  the 
parametric  directions.  The  above  example,  CCL, 
generates  an  element  with  cubic-cubic-linear 
displacement  functions.  The  default  is  CCC. 


1.  The  element  identification  number  must  be  the 
same  as  the  hyperpatch  identification  number. 

2.  See  Figures  4-1  and  2-6  for  the  parameteriza¬ 
tion  and  element  surface  orientation  genera 
ted  by  the  CPDE3  card. 

3.  Any  element  connected  to  a  reduced  element 
will  automatically  have  the  same  constraints 
(if  any)  on  their  common  suiface. 
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Input  Directive: 
Description: 

Format : 

E>  .mple  Syntax : 

Field 

EID 

Gl,  G2 , . . . ,  G8 

EM  PC 

Remarks: 


BULK  EATA  INPUT 


CPDE3  (Axisymnetric  model) 

Connectivity  card  for  <1  three-dimensional 
parametric  (finite)  discrete  dement. 

CPDE3,  EID,  Gl,  G2,  G3,  G4,,,,  EMPC,  G5,  G6, 
G7,  G8 

CPDE3,  7,  1,  3,  13,  12,,,,  CCX,  101,  103, 
113,  112 


Contents 

Element  identification  number. 

Corner  grid  point  identification  numbers  in  the 
sequence  given  by  Figure  4-1. 

Mesh  point  constraint  mnemonic  for  elements 
reduced  to  less  than  cubic  in  parametric  oi ruc¬ 
tions  one  or  two  and  axisymnetric  in  the  third 
parametric  direction.  The  above  example,  CCX, 
generates  an  element  with  cubic-cubic-constant 
displacement  functions.  The  default  is  CCC; 

i.e.,  not  axisymnetric. 


1.  Hie  axisymnetric  constraint  is  designed  by 
an  X  and  can  only  be  used  in  the  third 
parametric  coordinate  direction. 

2.  All  remarks  concerning  a  general  CPDE3 
element  apply. 

3.  If  the  PPDE3  card  for  an  axisymnetric  element 
references  a  nonisotropic  material,  then 
variable  Euler  angle  data  should  be  used  to 
generate  an  axisyrmietric  material  for  the 
element. 

4.  Surface  displacement  constraints  must  not  be 
input  for  mesh  points  not  on  face  five. 
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Commentary: 

The  connectivity  bulk  data  card,  nw.3,  causes  the  liyjjerpatch  for  ■> 
discrete  element  to  have  the  parameterization  shown  in  Figure  4-1, 
independent  of  how  it  was  constructed.  Any  reparameterization  that  may 
be  required  takes  place  automatically  and  prior  to  all  output.  H.e 
CPDE3  card,  for  example,  results  in  the  G1-G2-G3-G4  surface  beinq  face  S 
and  the  G5-G6-G7-G8  surface  beinq  face  A  when  the  hyperpatch  is  in 
geometric  format.  The  parameter i ration  is  defined  as  "3"  sort  in  that 
face  5  is  associated  with  f =  0  and  face  6  with  ,  -  1.  All  qeometrie 
hyper {xi t  dies  in  PATCHEG-III  are  in  "3"  sort.  At  the  present  time,  data 
hyperpatches  are  not  effected  by  the  CPDF3  card,  and  the  user  must 
ensure  that  their  pnrameterizations  are  consistent  with  the  geometric. 
Consider,  for  example,  a  temperature  hyperpatch  created  by  a  DHP2P  card 
using  temperature  patches  on  the  G1-G4-G8-G5  surface  and  G?.~G 3-G7-GH 
surface  for  face  5  and  face  6,  respectively.  This  data  would  be  in  "2" 
sort,  referring  to  Figure  4-1,  and  a  DHPSORT  card  would  be  required  to 
change  from  "2"  sort  to  "3"  sort. 

The  convention  fur  defining  an  element  surface  by  giving  the  corner 
grid  points  is  also  established  by  the  CPDE3  card.  This  convention, 
detailed  on  Figure  4-1,  may  be  thought  of  as  the  "left-hand  rule,"  in 
which  the  sequence  always  proceeds  from  the  origin  clockwise  about  an 
axis  in  the  a.  direction  where  the  surface  is  associated 

i 

with  r  =  constant.  This  convention  must  be  used  when  defining 
constraint  surfaces  and  data  patches  over  an  element  surface. 

:>  I 


BULK  DATA  OPTIONS 


PAGE  4-14 


G5  G6 


Face 


5 

G1 - G2-G3-G4 

"  Z(Cr  K2,  0) 

6 

C5-G6-G7-G8 

n  Z(Cr  1) 

3 

G1-G4-G8-G5 

=  Z(Cr  0,  ?3) 

4 

G2-G3-G7-G6 

-  7(C1 ,  1 ,  C3) 

1 

G1 -GF-G6-G2 

-  7(0,  03) 

2 

G4-G3-G7-G3 

r  z(i,  ?2,  c3) 

Figure  4-1.  Grid  Point  Conventions  for  Element 
Connectivity  and  Element  Surface 
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1  npnt  Directive: 
Description : 
Format : 

Example  Syntax : 

F  i  el  d 
DSID 

GIDK 

DK 

Remarks: 


BULK  DATA  INPUT 


LYiTAG  General  spatial  data  input  :  .y  grid  point 

Scalar  input  data  at  a  set  of  grid  [joints. 

DATAG,  IF, ID,,  GID1,  Dl,  GID2,  D2,  GIDl,  D3, 
GID4,  D4,  GID5,  D5,...,  GIDN,  EN 

QATAG,  11,,  3,  -4.0,  4,  -3.6,  3,  -3.2,  6,  -3.0 


Contents 

Identification  number  for  th<»  data  .sit  defined  by 
this  card.  Pressure,  temperature,  and  force  com¬ 
ponents  are  typical  data  types  that  form  an 
individual  data  set.  1  DPID%2  12. 

Identification  number  of  the  jrid  [joint  that 
locates  the  input  data  in  space. 

Data  value  at  grid  point  GIDK. 


1.  The  data  set  identif ication  number  must  be  an 
integer  one  through  twelve. 

2.  Data  at  up  to  29  grid  points,  can  lx>  input 
with  a  single  directive.  Additional  data 
in  the  set  can  be  input  with  additional 

d i recti ves . 

3.  Data  at  grid  points  in  a  data  set,  DBID,  can 
also  be  created  with  data  line  directives 
without  the  use  of  DATAG  input. 
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Input  Pi  recti ve: 
Description: 

Format : 

Example  Syntax 

Field 

ID 

FORMAT 

MTRX-ID 

Remarks: 


BULK  DATA  INPUT 


DHPAT  Data  Hyperpatch. 

Direct  input  of  a  one-component  data  hyperpatch 
in  any  PC  format. 

DHPAT,  ID,  FORMAT,  MTRX-ID 

DHPAT,  3,  P,  117 


Contents 

Data  hyperpatch  identification  number. 

A  or  S  Algebraic  coefficients. 

B  Geometric  coefficients. 

P  Point  coeff icients. 

G  Gaussian  coefficients. 

Matrix  identification  number  containing  the 
coefficients. 


1.  The  default  FORMAT  is  point  format  where 
the  mathematical  definitions  are  provided  in 
Chapter  5. 

2.  The  64  hyperpatch  coefficients  in  any  format 
are  input  in  the  same  sequence  that  a  triply 
subscripted  array  is  stored  in  FORTRAN; 
namely,  Pill,  P211,  P311,  P411,  P121,  P221, 
...  P344 ,  P444 . 
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Input  Directive: 
Description : 

Format : 

Example  Syntax: 

Field 

DUPID 

ns  id 

G1-G8 

Remarks : 
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BUI.K  DATA  INPUT 


DHPHEX  Data  hyperpatch  generator 

One-component  data  hyperpatch  created  from  the  8 
corner  values. 

DHPHEX,  DHPID,  ESID,  GI ,  G2,  G3,  G4,,,,,  GS,  Gfi, 
G7 ,  G8 

DHPHEX,  5,  2,  1,  2,  4,  6,,,,,  II,  3,  7,  8 


Contents 

Data  hyperpatch  identification  number. 

Identification  number  of  the  data  set  that  defines 
data  values,  Dl ,  at  the  grid  (joints. 

Corner  grid  points  that  locate  the  data  hyperpatch 
in  space  in  the  sequence  given  by  Figure  4-1. 


i.  Tiie  data,  Dl ,  are  interpolated  trilinearly. 

.1.  Tiie  data  in  a  data  set,  CHID,  are  normally 
input  with  a  DATAG  directive. 
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Description : 

Format : 

Example  Syntax: 

Field 

DHPL 

DL1-DL12 

CROSSF 

Remarks: 
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BULK  DATA  INPUT 


DHPL  Data  hyperpatch  generator 

One-component  data  hyperpatch  created  from  data 
lines. 

DHPL,  DHPID,  DL1 ,  DL2,  DL3,  1X4,  DL5,  DLfi, , 
CROSSF,  DL7,  DL8,  DL9,  DL10,  DL11,  DL12 

DHPL,  10,  3,  6,  7,  2,  1,  5,,  0,  8,  11,  12,  13, 
21,  18 


Contents 


Data  hyperpatch  identification  number. 

Data  line  identification  numbers  in  the  sequence 
specified  for  the  HPL  directive. 

0  =  Set  cross  derivatives  to  zero. 

1  =  Interpolate  cross  derivatives  linearly. 


1.  TT>e  data  lines  are  checked  for  conrnon  values 
at  the  corner  nodes. 

2.  The  data  lines  must  all  reference  the  same 
data  set. 

3.  The  CROSSF  parameter  controls  the  local 
surface  warping.  A  zero  indicates  a  data 
surface  with  zero  twist  at  the  four  corners. 
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Input,  Directive: 
Description: 

Format : 

Example  Syntax : 

Field 
DHPI D 
INSORT 

OUTSORT 

NCOMP 

Roma  i>;:: : 


BULK  DATA  INPUT 

DUPSORT  Resort  a  data  hyperpatch 

Transform  a  data  hyperpatch  from  the  inpjt  or 
construction  parameterization  to  the  output 
parameterization  requested. 

DHPSORT ,  DHPID,  INSORT,  OUTSORT,  NCOMP 

DHPSORT,  10,  -2 


Contents 

Data  hyperpatch  ID. 

Input  sort  format. 

TN50RT  =  ■ i  where  i  =  1,  7,  or  1  relative 

to  the  geometric  hyperpatch  implicitly 
associated  with  the  data  hyperpatch. 

Output  sort  format.  Default  -  3. 

Number  of  components  in  the  data  hyperpatch. 
Default  -  1. 


i.  Resetting  a  data  hyperpatch  is  a  user 

convenience  feature.  The  data  for  an  element 
may  be  easier  to  create  in  an  order  different 
from  that  used  for  the  geometry. 

A  minus  sign  reverses  the  two  faces 
associated  with  a  particular  sort. 

INSORT  =  -3,  for  example,  reverses  faces  5 
and  6. 
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Input  Directive: 
Description: 

Format  : 

Example  Syntax: 

Field 
Df.'PID 
DPI,  DP2 

Remarks: 


BULK  DATA  INPUT 


DHP2P  Data  hyperpatch  generator 

One-component  data  hyperpatch  created  from  two 
data  patches. 

CHP2P,  DHPID,  CPI,  DP2 

DHP2P,  3,  11,  12 


Contents 


Data  hyperpatch  identif ication  number. 

Data  patch  identification  numbers  for  F,  =  0 

and  F  ^  -  1 . 


1.  The  data  patches  must  both  reference 
(directly  or  indirectly)  the  same  data  set. 

2.  The  data  are  linearly  interpolated  between  the 
two  referenced  data  patct.es. 


3.  The  data  hyperpatch  is  in  "3"  sort  by 

construction.  Use  DHPSORT  if  another  sort 
is  desired. 
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Input  Directive: 
Description: 

Format : 

Example  Syntax: 

Field 

ID 

DPI,  2,  3,  4 

Remarks: 


BULK  DATA  INPLTT 


DHP4P  Data  hyperpatch  generator 

Data  hyperpatch  generated  from  four  patches  in 
one  parametric  direction. 

DHP4P,  ID,  DPI,  DP2,  DP3,  DP4 

DHP4P,  1,  1,  2,  8,  7 


Contents 


Data  hyperpatch  identification  number. 
Data  patch  identification  numbers. 


1.  This  construction  assumes  the  input  order 
of  the  four  data  patches  is  associated  with 
53  =  0,  1/3,  2/3,  1  surfaces  of  the  geometry 
hyperpatch. 
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Input  Directive: 
Description : 

Format : 

Example  Syntax : 

Fi_el(i 

DLID 

DSID 

Cl,  G2, . . . ,  O'! 

Remar  ks : 


BULK  DATA  IN  PITT 


DLINCS  Data  line  generator 

Scalar  data  interpolated  over  a  set  of  grid  points 
using  cubic  splines. 

DLINCS,  DLID,  DSID,  Gl,  G2,...,  GN 

DLINCS,  10,  3,  4,  2,  1,  5,  16 


Contents 

Identification  number  for  the  one-component  data 
line  created  by  this  card. 

Identification  number  of  the  data  set  that  defines 
data  values,  DI,  at  grid  points  in  the  data  line 
DLID. 

The  list  of  grid  points  that  locate  the  data 
interpolation  points  in  space.  The  maximum 
number  is  25. 


1.  The  data  DI  are  interpolated  using  a 
piecewise  cubic  spline  with  the  lineal  arc- 
length  between  grid  points  serving  as 
abscissa.  Data  must  exist  at  Gl,  GN,  and  at 
least  one  other  grid  point  and  may  exist  at 
all  grid  points. 

2.  The  data  in  data  set,  DGID,  are  usually 
input  with  a  DATAG  directive. 
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Input  Directive: 
Description: 
Format : 

Example  Syntax: 

Field 

DLID 

DSID 

FORMAT 

PI,  P2,  P3,  P4 
Gl,  G2,...,GN 

Remarks: 


BULK  DATA  INPUT 


DLINE  Data  line  generator 

Direct  input  of  a  data  line  in  any  format. 

DLINE,  DLID,  DSID,  FORMAT,  PI,  P2,  P3,  P4,, 

Gl,  G2, . . . ,  GN 

DLINE,  5,  1,  B,  1.,  12.,  0.,  -.5,  103  thru  120 


Contents 


Data  line  identification  number. 

Data  set  identification  number. 

A  or  S  Algebraic  coefficients. 

B  Geometric  coefficients. 

P  Point  coefficients. 

G  Gaussian  coefficients. 

Coefficients  of  the  line  in  the  specified  format. 
Grid  points  at  which  data  will  be  generated. 


1.  The  default  FORMAT  is  point  format,  P(0), 
P(l/3),  P(2/3),  P(l) .  The  sequence  is  the 
same  for  all  formats  as  described  in  Chapter 

5. 

2.  The  data  line  extends  from  Gl  to  GN  with 
intermediate  values  determined  by  the  distance 
between  points. 
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Input  Directive: 
Description : 

Format : 

Example  Syntax : 

Field 
DLID 
DSID 
Gl,  G2 

D ( 1/3 ) ,  D(2/3) 


BULK  DATA  INPUT 


DLINP  Data  line  generator 

Scalar  data  for  a  data  line  in  point  format 
between  two  grid  points. 

DLINP,  DLID,  DGID,  Gl,  G2,  D(l/3),  D(2/3) 

DLINP,  10,  5,  3,  6,  .4,  -.2-1 


Contents 
Data  line  ID. 

Data  set  ID  that  defines  the  data  at  Gl,  G2. 

Grid  points  that  locate  the  ends  of  the  data  line 
in  space. 

Data  values  at  the  one-third  points  in  parametric 
space.  The  spatial  location  of  these  points  is 
determined  by  the  discrete  element  whose  edge  is 
associated  with  this  data  line. 
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BULK  DATA  INPUT 


Input  Directive: 

DLINPC  Data  line  generator 

Description: 

Scalar  data  interpolated  over  a  set  of  grid  points 
as  piecewise  linear. 

Format : 

DLINPC,  DLID,  GRID,  Gl,  G2,...,  GN 

Example : 

DLINPC,  7,  4,  14,  3,  1,  5,  6 

Field 

Contents 

DUD 

Identification  number  for  the  one-component  data 
line  created  by  this  card. 

DU  ID 

Identification  number  of  the  data  set  that  defines 
data  values,  DI,  at  grid  points  in  the  data  line 
DLID. 

Gl,  G2, . . . ,  GN 

The  list  of  grid  points  that  locate  the  data 
interpolation  points  in  space.  The  maximum 
number  is  25. 

Remarks : 

1.  The  data,  DI,  are  interpolated  as  piecewise 
linear  with  the  lineal  arclength  between  grid- 
points  serving  as  abscissa.  Data  must  exist 
at  Gl  and  Gm'  and  may  exist  at  all  grid  points. 
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Input  Directive 
Description: 

Format: 

Example  Syntax : 

Field 

DPATID 

CJ1 ,  2,  3,  4 

Sll ,  S21  ,S3l , . . . , 


BULK  DATA  INPUT 


DPATA  Data  patch  generator 

Algebraic  format  data  patch  on  one  surface  of 
an  element. 

DPATA,  DPATID,  Gl,  G2,  G3,  G4,,,,  Sll,  S21,  S31, 
S41,  S12,  S22,  S32,  S42,  S13,  S23,  S33,  S43,  S14, 
S24 ,  S34 ,  S44 

DPATA,  10,  1,  3,  16,  4,,,,,,  -1.3,  2.4,,,  1.0,  1.0 


Contents 


Data  patch  ID. 

Corner  grid  points  that  locate  the  data  patch  in 
space. 

S44  Data  patch  in  algebraic  format  where  a  blank  is 

equivalent  to  a  zero  SIJ.  The  spatial  location  of 
these  data  points  is  determined  by  the  discrete 
element  whose  surface  uses  the  data  patch. 


F(^,c2)  =  e,,i) 

"sn 

S12 

S13 

S14~ 

( 

\  2 

S21 

S22 

S23 

S24 

V 

l  7 

S31 

S32 

S33 

S34 

)  6 

1  7 

S41 

S42 

S43 

S44_ 

\  1 

04 
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Input  Directive: 
Description: 

Format : 

Example  Syntax: 

Field 

ID 

FORMAT 

MTRX-ID 

Gl,  G2,  G3,  G4 

Remarks : 


BULK  DATA  INPUT 


DPATCH  Data  patch 

Direct  input  of  a  one-component  data  patch  in 
any  format. 

DPATCH,  ID,  FORMAT,  MTRX-ID,  Gl,  G2,  G3,  G4 
DPATCH,  4,  B,  22 


Contents 

Data  patch  identification  number. 

A  or  S  Algebraic  coefficients. 

B  Geometric  coefficients. 

P  Point  coeff icients. 

G  Gaussian  coefficients. 

Matrix  identification  number  containing  the 
coefficients. 

Data  patch  spatial  location  given  by  corner  grid 
points. 


1.  The  default  FORMAT  is  point  format. 

2.  The  16  patch  coefficients  are  input  in  the 
sequence  Pll,  P21,  P31,  P41,  PI 2,  P22,  ... 
for  all  formats  as  described  in  Chapter  5. 
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BULK  EATA  INPUT 


Input  Directive: 

DPATEQ  Data  patch  generator 

Description: 

Creates  a  data  patch  equal  to  a  reference  data 
patch  but  located  on  a  different  surface. 

Format : 

DPATEQ,  DPATID,  REFID,  Gl,  G2,  G3,  G4 

Example  Syntax: 

DPATEQ,  3,  6,  1,  2,  5,  11 

Field 

Contents 

DP AT ID 

ID  of  the  data  patch  to  be  created. 

REFID 

ID  of  the  reference  data  patch. 

Gl,  G2,  G3,  G4 

Corner  grid  points  of  the  surface  on  which  data 
patch  DPATID  will  be  located. 
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Jlngul  Directive: 
Description : 

Format : 

Example  Syntax : 

Field 

DPID 

Gl,  G2,  G3,  G4 
DL1 ,  DL2,  DL3,  DL4 

Hemar  k$: 


BULK  DATA  INPUT 


DPATL  Data  patch  generator 

One-component  data  patch  created  from  data  liner. 

DPATL,  DPID,  Gl,  G2,  Gl,  G4,,,,,  DI.l ,  DL2,  1X3, 
DL4 

DPATL,  14,  5,  6,  10,  12,,,,,  3,  4,  7,  2 


Contents 

Data  patch  identification  number. 

Corner  grid  point  identification  numbers. 
Data  line  identification  numbers. 


1.  The  data  lines  are  checked  for  common  values 
at  the  corner  nodes. 

2.  The  data  lines  must  all  reference  the  same 
data  set. 

1.  DATAG  input  is  not  mandatory  if  the  data  lines 
are  direct  input.  The  grid  point  values  will 
be  determined  from  the  data  lines. 

4.  The  grid  point  and  line  sequencing  is  as 
shown  for  the  PATCHL  directive. 
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BULK  DftTA  INPUT 


Input  Directive: 

DPATQ  Data  patch  generator 

Description: 

One-component  bilinear  data  patch  created  from 
corner  values. 

Format : 

DPATQ,  DPID,  DSID,  Gl,  G2,  G3,  G4 

Example  Syntax: 

DPATQ,  2,  6,  6,  11,  7,  3 

Field 

Contents 

DPID 

Identification  number  for  the  one-component  data 
patch  created  by  this  card. 

DSID 

Identification  number  for  the  data  set  that  has 
data  values,  DI,  at  the  corner  grid  points. 

Gl,  G2,  G3,  G4 

Corner  grid  points  that  locate  the  data  patch  in 
space. 

Remarks: 

1.  See  PATCHL  directive  for  grid  point  sequence 
illustration. 
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Input  Pi  recti ve : 
Description: 

Format: 

Example  Syntax: 

Field 

DPIP 

Gl,  2,  3,  4 
ULI /S EG 

Remarks : 


BULK  LATA  INPUT 


DPAT4L  Data  patch  generator 

One-component  data  patch  generated  from  single 
line  segments  of  four  data  lines  in  o<'e  parametric 
direction. 

DPAT4L,  DPID,  Gl ,  G2,  G3,  G4,,,,,  DL1/SEG, 

DL2/SEG,  DL3/SEG,  DL4/SEG 

DPATL4 ,  6,  5,  4,  7,  9,,,,,  5,  19,  21/3,  22/4 


Contents 


Data  patch  identification  number. 

Corner  grid  point  identif ication  numbers. 

Data  line  number  for  line  I,  segment  number  SBG. 
(Default  SEG  is  1.) 

1.  See  PATCHL  directive  for  grid  point  sequence 
illustration. 
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Input  Directive: 
Description: 

Format : 

Example  Syntax : 

Field 

ID 

Tref 

Ti 

f  (Ti) 

Remarks: 


BULK  DftTA  INPUT 

DTCS  Data  table  function 

Defines  a  scalar  function  from  tabular  data 
using  piecewise  cubic  spline  interpolation. 

DTSC,  ID,  Tref,  TI,  f(T 1),  T2,  f (T2) ,  T3,  f(T3) 

7*,,  100.0,  10.1E6,  200.0,  9.8E6,  250.0,  8.4E6 

Contents 

Table  identification  number.  If  followed  by  an 
asterisk  the  function  is  normalized  such  that 
f*  (Tref)  =  1.0. 

Reference  abscissa,  TI  ^  Tref.<  TO.  The  default 
value  is  Tref  =  TI. 

Abscissa  values  in  increasing  order,  Ti  >  Ti-1. 
Value  of  the  function  f(T)  at  T  =  Ti . 


1.  The  function  f(T)  is  interpolated  by  a 
piecewise  cubic  spline.  Off  table  values  are 
set  to 

f (T)  2  f (TI)  for  T  <  TI 
f  (T)  5  f  (TO)  for  T  >  TO 

2.  Attempts  to  normalize  a  function  with 
f(Tref)  =  0.0  are  fatal  errors. 
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Input  Directive: 
Description: 

Format : 

Example  Syntax: 

Field 

ID 

Tref 

Ti 

f  (Ti) 

Remarks: 


BULK  EATA  INPUT 


DTPC  Data  table  function 

Defines  a  scalar  function  from  tabular  data 
using  piecewise  linear  interpolation. 

UTPC,  ID,  Tref,  TI,  f (TI) ,  T2,  f(T2),  T3, 
f(T3),...,  TO,  f (TN) 

CfTPC,  6*,  160,  60.,  30.0E6,  300.,  26.0E6,  1000., 
10.0+6,...,  5000,  1.0+6 


Contents 

Table  identification  number.  If  followed  by  an 
asterisk  the  function  is  normalized  such  that 
f*  (Tref)  =  1.0. 

Reference  abscissa,  Tl  <  Tref  <  TO.  The  default 
value  is  Tref  =  Tl. 

Abscissa  values  in  increasing  order, 

Ti  >  Ti-1. 

Value  of  the  function  f(T)  at  T  =  Ti . 


1.  The  function  f (T)  is  interpolated  as  piece- 
wise  linear.  Off  table  values  are  set  to 


f  (T)  =  f  (Tl)  for  T<  Tl 

f  (T)  =  f  (TO)  for  T  >TO 

2.  Attempts  to  normalize  a  function  with 
f(Tref)  =  0.0  are  fatal  errors. 
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BULK  DATA  INPUT 


Input  Directive: 
Description: 

Format: 

Example  Syntax : 


FORCE  Static  load 

Defines  a  static  load  at  a  grid  point  by  Cartesian 
vector  components. 

FORCE,  ID,  GID,  FI,  F2,  F3 

FORCE,  5,  21,  500.,,  600. 


Field  Contents 

ID  Load  set  identification  number  (integer,  1  to 

999). 


GID 


Grid  point  identification  number. 


FI,  F2,  F3 


Components  of  the  load  vector  in  the  reference 
Cartesian  coordinate  directions  eir  e2,  and  e3. 
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Input  Directive: 
Description: 
Format : 

Example  Syntax: 

Field 

ID 

EID 

DL1,  2,  3 
SI,  2,  3 

Remarks: 


BULK  DATA  INPUT 

F0RCEL3  Line  load  data 

Defines  a  line  load  on  one  edge  of  a  3-D  element. 
FORCED 3,  ID,  EID,  DL1,  DL2,  DL3,  SI,  S2,  S3 
F0RCEL3,  6,  10,  2,,  4,,  2.0 


Contents 


Load  set  identif ication  number. 

Element  identification  number. 

Data  lines  for  components  FI,  F2,  and  F3  where 
F  =  Fiai  and  the  as  are  unit  tangent  vectors. 

Scale  factors  applied  to  the  data  lines  DLl,  2,  3. 
Default  values  are  1.0. 


1.  The  data  lines  must  reference  the  same 

grid  points,  and  these  are  used  to  locate  the 
edge  of  element  EID  on  which  the  line  load 
acts.  A  blank  field  for  a  data  line  ID 
indicates  that  component  of  the  load  is 
identically  zero. 
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Input  Directive: 
Descr  iption: 

Format: 

Example  Syntax: 

Field 

ID 

EID 

DPI,  2 

SI,  2 

Remai  ks : 


BULK  DATA  INPUT 


PORCET  Surface  force  data 

Defines  a  surface  traction  on  one  surface  of  a 
3-D  element. 

PORCET,  ID,  EID,  DPI,  DP2,  SI,  S2 
PORCET,  6,  10,  11,,  -2.3 


Contents 

Load  set  identification  number. 

Element  identification  number. 

Data  patches  for  components  T1  and  T2  of  the 
traction  T  =  Tlt^  +  T2tz  where  tj,  t,  are  unit 
tangent  vectors T 

Scale  factors  applied  to  the  data  patches  DPI,  2. 
Default  values  are  1.0. 


1.  The  data  patches  must  reference  the  same 
grid  points,  and  these  are  used  to  locate  the 
face  of  element  EID  on  which  the  traction 
acts.  A  blank  field  for  DPI  (DP2)  indicates 
T1  (T2)  is  identically  zero. 

2.  The  sense  of  the  unit  tangent  vectors  is  given 
in  Figure  2-6. 
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BULK  DATA  INPUT 
GRID  Grid  point 

Defines  the  coordinates  of  a  grid  point. 

GRID,  GID,  CID,  Zl,  Z2,  Z3 
GRID,  5,  1,  0.0,  2.0,  4.05 

Contents 

Grid  point  identification  number. 

Coordinate  type.  1-Cartesian,  2-Polar 
(default  =  1,  integer) . 

Location  of  the  grid  point  in  coordinate  system 
CID. 

Remarks:  1.  The  input  sequence  for  polar  coordinates 

is  r  ,  ft  ,  7. . 


input  D i recti ve : 

Description: 

Format: 

Example  Syntax : 

Field 

GID 

CID 

7.1,  7.2,  73 
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Input  Directive: 
Description: 

For  it.  at : 

Example  Syntax: 

Field 

ID 

FORMAT 

Ml,  P,  3 
TID 

Remarks: 


BULK  DATA  INPUT 


HPATCU  Hyperpatch  generator 

Direct  input  of  a  hyperpatch  in  any  format. 

HPATCH ,  ID,  FORMAT,  Ml,  M2,  M3,  TID. 

HPATCH ,  1,  P,  101,  102,  103. 


Contents 

Hyperpatch  identification  number. 

A  or  S  Algebraic  coeff icien*'s. 

B  Geometric  coefficients. 

P  Point  coefficients. 

G  Gaussian  coefficients. 

MTRX  identi  f ication  numbers  for  7.1,  7.P,  and  73. 

Transformation  identification  number  to  be 
applied  to  the  hyperpat.ch,  if  any. 


1.  The  default  FORMAT  is  point  format. 

2.  The  64  coefficients  in  any  format  for  each 
hyperpatch  coordinate  function  are  input  in 
the  sequence  Pill,  P211,  P311,  P411, 

P121,  .... 

3.  The  mathematical  definition  of  the  different 
format  coefficients  is  defined  in  Chapter  a. 
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BULK  DATA  INPUT 


Input  Directive: 

HPHEX  Hyperpatch  for  a  linear  hexahedron 

Description : 

Generates  a  hyperpatch  from  the  8  corner 
grid  points. 

Format: 

HPHEX,  HPAT1D,  Gl,  G2,  G3,  G4,,,,,  G5,  G6,  G7,  G8 

Example  Syntax: 

HPHEX,  1,  1,  2,  3,  4,,,,,  5,  6,  7,  8 

Field 

Contents 

HP AT ID 

Hyperpatch  identification  number. 

G1 

The  grid  point  identification  number  for  the 
(0,  0,  0)  node. 

G2 

The  grid  point  identification  number  for  the 
(0,  1,  0)  node. 

G  i 

The  grid  point  identification  number  for  the 
(1,  1,  0)  node. 

G4 

The  grid  point  identification  number  for  the 
(1,  0,  0)  node. 

G5 

Tiie  grid  point  identification  number  for  the 
(0,  0,  1)  node. 

G6 

The  grid  point  identification  number  for  the 
(0,  1,  1)  node. 

G7 

The  grid  point  identification  number  for  the 
(1,  1,  1)  node. 

G8 

The  grid  point  identification  number  for  the 
(1,  0,  1)  node. 

Remarks: 

1.  The  trilinear  hexahedron  is  not  required 

to  have  flat  surfaces. 
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Commentary: 

Generation  of  a  hyperpatch  from  the  eight  corner  points. 


6 


Method : 

First  generate  the  lower  and  upper  surface  patches  B1  and  B2  using 
PATCH 1 ,  and  then  use  the  technique  outlined  in  the  description  of  HP2PAT 
to  compute  B3  and  B4  as 

[B3]  =  [B4]  =  [B2]  -  [Bl] 
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Lnput  Directive; 
Description: 
Format : 

Example  Syntax: 

Field 
HPATID 
LI  -  L12 
CROSS F 


Remarks : 


BULK  DATA  INPUT 

HPL  Gridline  definition  of  a  hyperpatch 

Generates  a  hyperpatch  from  the  12  edge  lines. 

HPL,  HPATID,  LI,  L2,  L3,  L4,  L5,  L6,,  CROSSF,  L7, 
L8,  L9,  L10,  Lll,  LI 2 

HPL,  1,  4,  12,  8,  9,  2,  11,,,  6,  10,  1,  3,  7,  5 


Contents 


Hyperpatch  identification  number. 

Line  identification  numbers,  as  shown  below. 

Cross  derivatives  flag.  0  for  0.  cross 
derivatives,  1  for  linearly  interpolated  cross 
derivatives. 


LI  =  Z(0,f,2,0) 
L2  =  Z(?,,1,0) 
L3  =  Z(1,5?,0) 
L4  =  Z(f>1,0,0) 
L5  =  Z(0,C2,1) 
L6  =  Z  ( F,  j ,  1 , 1 ) 


L7  =  Z(l,52,l) 
L8  =  Z  ,0,1) 
L9  =  Z(0,0,^,) 
L10  =  Z  (0,1,  £3) 
Lll  =  Z(l,l,f.j) 
L12  =  Z(1,0,F.3) 


1.  The  CROSSF  parameter  controls  the  local 
surface  warping.  A  zero  indicates  a  data 
surface  with  zero  twist  at  the  four  corners. 
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Method : 

Generate  surface  patch  PI  from  the  lines  LI,  L2,  L3,  L4. 

Generate  surface  patch  P2  from  the  lines  L5,  L6,  L7,  L8. 

Generate  surface  patch  P3  from  the  lines  L4,  L12,  L8,  L9. 

Generate  surface  patch  P4  from  the  lines  L2,  Lll,  L6,  L10. 

Generate  surface  patch  P5  from  the  lines  L9,  L5,  L10,  LI. 

Generate  surface  patch  P6  from  the  lines  L12,  L7,  Lll,  L3. 

Use  the  procedure  described  for  HP6PAT  to  compute  Bl,  B2,  B3,  B4 
from  six  patches  after  checking  all  lines  for  common  end  points  at  the 
corners. 


80 


BULK  DATA  OPTIONS 


Input  Directive: 
Description: 

Format: 

Example  Syntax: 

Field 

HPATID 

PATID 

THK 

Remarks: 
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BULK  DATA  INPUT 


HPN  Hyperpatch  for  thick  shells 

Generates  a  hyperpatch  from  a  patch  and  a  given 
thickness  in  the  normal  direction. 

HPN,  HPATID,  PATID,  TOK 

HPN,  1,  7,  1.675 

Contents 

Hyperpatch  identification  number. 

Base  patch  identification  number. 

Thickness  of  the  hyperpatch. 


1.  When  used  with  degenerate  patches,  the  normals 
near  the  degenerate  node  may  have  some  dis¬ 
tortion. 
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Commentary: 

Generation  of  a  hyperpatch  from  a  patch  and  a  given  thickness  in 
the  normal  direction. 


Thickness 


Method : 

Transform  B  to  point  format,  [Z]. 

Compute  the  unit  normals  at  16  points  on  the  patch,  [N] . 

If  t  >  0  define 

B1  =  [Z],  B2  :  [Z]  +  t  [N] 

If  t  <  0  define 

B1  MZ]  +  t  [N],  B2  5  [Z] 

Use  the  procedure  described  for  HP2PAT  to  compute  B3  and  B4. 


32 


BULK  DATA  OPTIONS 


PAGE  4-45 


Input  Directive: 
Description: 

Format : 

Example  Syntax : 

Field 

HPATID 

PATID 

Zl,  Z2,  Z3 

THETAB,  TOETAE 


BULK  DATA  INPUT 

HPR  Hyperpatch  for  a  body  of  revolution  segment 

Generates  a  hyperpatch  by  rotating  a  planar  patch 
about  a  coordinate  axis. 

HPR,  HPATID,  PATID,  Zl,  Z2,  Z3,  THETAB,  THETAE, 
IDRA 

HPR,  7,  3,  0.0,  1.0,  -2.5,  15.0,  40.0,  2 
Contents 

Hyperpatch  identification  nunber. 

Identification  number  of  the  patch  to  be  rotated. 

Coordinates  for  shifting  the  origin. 

Beginning  and  ending  angles  through  vshich  the 
patch  will  be  rotated  (real,  degrees). 


IDRA 


Axis  of  rotation  (integer,  -3  to  3). 

+1  *  +X,  -1  -  -X,  +2  *  +Y,  +3  -  +Z,  etc. 
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Commentary: 

Generation  of  a  hyperpatch  by  rotating  a  planar  patch  about  a 
coordinate  axis, 


Method : 

Given  the  planar  patch  Z(£  ,?2),  rotate  the  lines  Z(0,£2), 
Z(1/3,C2),  Z(2/3,S2),  and  Z(1,S2)  to  form  tha  patches  61,  62,  63,  and  64 
using  PATCHR. 

Then  compute  the  patch  coefficients. 

[Bl]  =  [81]  ,  [B2]  =  [84] 

[B3]  =  -5.5  [81]  +  9 [82]  -4.5{63]  +  [84] 

[B4]  =  -[81]  +4.5[B2]  -9[83]  +5.5[84] 


34 


BULK  DATA  OPTIONS  PAGE  4-47 

The  parameterization  generated  is  as  follows:  If  a*  and  a*  are  the 
tangent  vectors  of  the  original  patch  and 

(1)  If  ei  x  Z (1/2, 1/2)  •  N (1/2, 1/2)  >  0  then  a*  =  a3,  a*  =  a; 

(2)  If  e,.  x  Z (1/2, 1/2)  •  N (1/2, 1/2)  <  0  then  a*  =  a  ,  a*  =  a 

~ 1  ~  -  -  1  .  ^  .  9  . 9 
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Input  Directive: 
Description: 

Format: 

Example  Syntax : 

Field 

HPATID 

PAT1 

PAT  2 

Remarks: 
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BULK  DATA  INPUT 

HP2PAT  Hyperpatch  for  interpolating  two  surfaces 

Generates  a  hyperpatch  from  two  surface  patches 
with  linear  interpolation  of  the  hyperpatch 
parameters  between  faces  five  and  six. 

HP2PAT,  HPATID,  PAT1,  PAT 2 

HP2PAT,  1,  11,  12 


Contents 


Hyperpatch  identification  number. 

Surface  patch  identification  number  of  the 
C3  =  0  patch. 

Surface  patch  identification  number  of  the 
5 3  =  1  patch. 


1.  Ttie  hyperpatch  is  in  "3"  sort  by  construction. 
It  cannot  be  resorted  automatically. 
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Commentary: 


Cicnot.it  ion  of  i  hyjtorpafeh  I  rom  two  surf  ace  patches  wit  h  linen; 
interpolation  of  tiio  hyperpatch  parameters  between  the  faces. 


Method: 


71 

72 

(Z,.  )1 
'  1 

(Z,  )  2 

Patch  1,  [Bl 1  = 

Z4 

73 

(7.f  r  )  4 

( Z ,  f  )  3 

(Z'f 

)1 

1 

(Z,,  ) 2 

M 

(Z,  . 

i ,  ■ 

)4 

i 

(Z,r  ) 3 
f  1 

'  Zh 

Zf> 

(Z,,  )h 

(7,  )h 

7,  . 

Patch  2,  [B2)  = 

78 

7.1 

(Z,  .  ) 8 

(7  /  )  7 

'■? 

<S 

)c> 

(Z»  )6 

■  1 

(Z.f  .  )‘> 

'  !'■;> 

(Z,..  ,  )h 

<S, 

)8 

(Z.  )  7 

f  1 

N 

.fi 

**> 

N3 

00 

(Z,  )  7 

To  make  the  hyperpatch  linear  between 
'/  ( i  [ ,  1 )  compute 


faces  7  ( t  , 


0)  and 


[B3]  =  [B4]  =  [B2]  -  [Bl] 
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BULK  DATA  INPUT 


Input  Pi  recti ve: 
Description: 

Format : 

Example  Syntax: 


HP4PAT  Hyperpatch  generator 

Hyperpatch  generated  from  four  director  patches 
in  parametric  direction  three. 

HP4PAT,  ID,  PI,  P2,  P3,  P4,  TID 

HP4PAT,  1,  3,  4,  8,  7 


Field 

ID 

PI 

P2 

P3 

P4 

TID 


Contents 

Hyparpetch  identification  number. 

Patch  identification  number  for  f,3=  0. 

Patch  identification  number  for  £3  =  1/3. 

Patch  identification  number  for  £•,=  2/3. 

Patch  identification  number  for  S3=  1. 

Transformation  identification  number,  if  any, 
to  be  applied  to  the  hyperpatch. 
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BULK  DATA  INPUT 


Inpi  it  1)1  rectlve: 
Description: 


Format : 


Example  Syntax: 


HPRPAT  Hyperpatch  for  a  curvilinear  hexahedron 

Generates  a  hyperpatch  from  6  patches  which 
enclose  a  volume. 

HP6PAT,  HPATID,  Pi,  P2,  P3,  P4,  P5,  P6,  CROSSF 
HP6PAT,  1,  1,  2,  3,  5,  6,  4,  0 


Field 


Contents 


HPATID 


CROSSF 


Hyperpatch  identification  number. 

Surface  patch  identification  number  of  the 
=  0  patch. 

Surface  patch  identification  number  of  the 
f  =  1  patch. 

Surface  patch  identification  nunber  of  the 
=  0  patch. 

Surface  patch  ident if ication  number  of  the 
f,  =  1  patch. 

Surface  patch  identification  number  of  the 
Cj  =  0  patch. 

Surface  patch  identification  number  of  the 
Lj  =  1  patch. 

Cross  derivatives  flag.  0  for  0.  cross  deriva¬ 
tives  and  1  for  linearly  interpolated  cross 
derivatives. 


Remarks; 


1.  The  CROSSF  parameter  controls  the  local 
surface  warping.  A  zero  indicates  a  geometric 
surface  with  zero  twist  at  the  four  corners. 
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Commentary: 

Generation  of  a  hyperpatch  from  6  patches  which  enclose  a  volume. 


Method : 

Given  the  patches  Pi,  P2,  P3,  P4,  P5,  and  P6,  form  the  hyperpatch 
coefficients  Bl,  B2,  B3,  B4  as  follows: 

[Bl]  =  [PI] 

[B2]  =  [P2] 

Coefficients  for  B3  and  B4  are  extracted  from  P3,  P4,  P5,  and  P6 

after  checking  all  patches  for  common  lines  on  intersecting  edges. 
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Input  Directive: 
Description: 

Format: 

Example  Syntax: 

Field 

LID 

R 

71,  7.2,  73 
THETAB,  THFTAE 

I  DR/'. 

G ;  ,  O  GN 

Remarks : 


BULK  DATA  INPUT 


LARCPC  A  circular  arc  with  N-l  line  segments 

Generates  a  piecewise  parametric  cubic  line 
from  a  circular  arc  with  N-l  line  segments. 

LARCPC,  LID,  R,  Zl,  Z2,  Z3,  TOETAB,  THETAE,  IDRA, 
Gl,  G2,  G3, . . . ,  GN 

LARCPC,  21,  .5,  0.0,  3.2,  .006,  30.0,  60.0,  3,  2, 
5  THRU  10,  12 


Contents 


Line  identi f ication  number  for  the  generated  line 
( integer ,  1  to  180) . 

Radius  of  the  arc. 

Coordinates  of  the  center  of  the  circle. 

Beginning  and  ending  angles  of  the  sweep  of  the 
arc  (real,  degrees). 

Axis  of  rotation  (integer,  -3  to  3). 

+1  =  +X,  -1  =  -X,  +2  =  +Y,  +3  =  +7,  etc. 

Gl  is  the  first  grid  point  on  the  arc,  GN  is  the 
last  grid  point  on  the  arc,  and  the  interior  items 
in  the  list  are  the  equally  spaced  interior  grid- 
points. 


1.  See  PATCHGR  commentary  for  an  illustration  of 
the  construction  method. 
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Input  Pi  recti ve: 
Description: 

Format: 

Example  Syntax : 

Field 

LID 

FORMAT 

Ml,  2,  3 
TID 

Gl,  2,...,  N 

Remarks: 


BULK  DATA  INPUT 


LINE  Direct  input  of  a  parametric  cubic  line. 

Direct  input  of  a  parametric  cubic  line  in  any 
format.  The  line  is  subdivided  into  N-l  line 
segments  with  N  grid  points. 

LINE,  LID,  FORMAT,  Ml,  M2,  M3,  TID,  Gl ,  G2,  G3, 
G4, . . . ,  GN 

LINE,  7,  B,  201,  202,  203,,  21,  24,  25,  30,,,  100 


Contents 


Line  identification  nunber. 

A  or  S  Algebraic  coefficients. 

B  Geometric  coefficients. 

P  Point  coefficients. 

G  Gaussian  coeff icients. 

Matrix  identification  nunbers  for  Zl,  Z2,  and  Z3. 

Transformation  identification  nunber,  if  any,  to 
apply  to  the  line. 

Grid  point  identification  nunbers  for  the  N  grid- 
points  on  the  line. 


TTie  default  FORMAT  is  point  format.  The 
mathematical  definition  of  the  input  coefficients 
for  each  format  is  in  Chapter  5. 
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Input  Directive: 
Description: 

Format : 

Example  Syntax : 

Field 

LID 

GI  ,  G  2 , .  •  •  ,  GN 

Remarks: 


BULK  DATA  INPUT 


LINECS  Parametric  spline  line 

Generates  a  piecewise:  parametric  cubic  line  from 
the  input  grid  points  and  spline  constraints. 

LINECS,  LID,  Gl,  G2,  G3,...,  GN 

LINECS,  5,  2,  4  THRU  10  EXCEPT  7 


Contents 


Line  identif ication  number  for  the  generated  line 
(integer,  1  to  180). 

List  of  grid  points  through  which  the  spline  is  to 
be  passed. 


There  must  be  at  least  three  grid  points,  N  3*  3, 
to  use  this  line  card.  There  can  be  a  maximnn 
of  25. 
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Input  Directive: 
Description: 

Format : 

Example  Syntax: 

Field 

LID 

GP 

ZA1,  2,  3 

ZB1,  2,  3 

TID 

GAMMA,  GAMMO 

Gl,  2,...,  N 

Remarks: 


BULK  DATA  INPUT 


LINEGR  A  line  from  a  general  rotation 

Generation  of  a  piecewise  parametric  cubic  line 
from  the  general  rotation  of  a  point  about  an 
arbitrary  axis. 

LINEGR ,  LID,  GP,  ZA1,  ZA2,  ZA3,  ZB1,  ZB2,  ZB3, 
TID,  GAMMA,  GAMMO,  Gl,  G2,...,  GN 

LINEGR,  6,  21,  0.,  0.,  0.,  1.,  1.,  -2.5,,  30., 
0.,  1  THRU  6,  50 


Contents 

Line  identification  number. 

Grid  point  to  be  rotated. 

Coordinates  of  the  base  of  the  rotation  vector 
(axis)  . 

Coordinates  of  the  head  of  the  rotation  vector 
(axis)  . 

Transformation  identification  number,  if  any,  to 
be  applied  to  the  line. 

Subtended  angle,  initial  offset  angle.  Same  as 
PATCHGR. 

The  grid  point  identification  numbers  for  the  N 
grid  points  on  the  line. 


1.  The  default  offset  angle  is  zero. 

2.  See  the  PATCHGR  directive  for  an  illustration 
of  the  construction. 
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Input  Directive: 
Description: 

Format : 

Example  Syntax : 

Field 

LID 

G1 ,  G2, . .  • ,  GN 

Remarks: 


BULK  DATA  INPITT 


LINEPC  Generate  a  straight  line  with  N  -  1  line 
segments 

Generates  a  piecewise  parametric  cubic  line  with 
N  -  1  line  segments. 

LINEPC,  LID,  Gl,  G2,  G3,...,  QJ 

LINEPC,  5,  1,  5,  8  THRU  16  EXCEPT  14 


Contents 


Line  identification  nunber  for  the  generated  line 
(integer,  1  to  180). 

Gl  is  the  starting  grid  point  for  the  line.  GN  is 
the  ending  grid  point  for  the  line. 


1.  The  coordinates  of  the  intermediate  grid 
points  are  automatically  computed  and  are 
uniformly  spaced. 

2.  The  end  points  Gl  and  must  be  input  or 
constructed  by  another  directive. 
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BULK  DATA  INPUT 


Input  Directive: 
Description: 


Format : 


MATAL  Thermal  expansion  coefficients 

Matrix  identification  numbers  for  the  three- 
dimensional  thermal  expansion  coefficients  at  N 
points  in  an  element. 

MATAL,  MID,  FRAME,  POINTS,  AL1,  AL2, ALN 
MATAL,  7,  1,  8,  101  THRU  108 


Field  Contents 

MID  Material  identification  number. 

FRAME  =  1,  Thermal  expansion  coefficients  £  are  in  an 

orthonormal  Cartesian  frame. 

=  2,  Thermal  expansion  coefficients  J  are  in  the 
normalized  parametric  frame. 

POINTS  =  1,  Constant  properties. 

=  8,  Trilinear  variation  in  properties. 

=  64,  Tricubic  variation  in  properties. 

AL1,  AL2, . . . ,  ALN  Matrix  identification  numbers  for  the  thermal 

expansion  coefficients  a  at  the  interpolation 
points.  If  N  =  64  a  single  entry  is  used  to 
identify  a  matrix  containing  the  64  identifica¬ 
tion  numbers. 


Remarks:  I.  The  alJ  are  input  in  the  sequence 

all,  a22,  a33,  a!2,  «13,  <»23. 


96 


BULK  DATA  OPTIONS 


PAGE  4-59 


Input  Directive: 
Description: 

Format: 

Example  Syntax: 

Field 

MID 

FRAME 


POINTS 

Cl,  C2, . . . ,  CN 


Remarks: 


BULK  DATA  INPUT 


MATC  Materials  stiffness  matrix 

Matrix  identification  numbers  for  the  three- 
dimensional  stress-strain  equations  at  N  points  in 
an  element. 

MATC,  MID,  FRAME,  POINTS,  Cl,  C2,...,  CN 
MATC,  7,  2,  8,  6  THRU  12,  16 


Contents 


Material  identification  number. 

=  1,  C  Matrix  is  in  an  orthonormal  Cartesian 
frame. 

=2,  C  Matrix  relates  contravariant  physical 

components  of  stress  to  covariant  physical 
components  of  strain  in  the  parametric 
frame.  (Useful  with  orthotropic  bodies  of 
revolution.) 

=  1,  Constant  material  properties  (N  =  1). 

=  8,  Trilinear  variation  of  material  properties 

(N  =  8). 

=  64,  Tricubic  variation  of  material  properties 
(N  =  64 ) . 

Matrix  identification  numbers  for  the  stress- 
strain  coefficient  matrices  at  the  interpolation 
points.  If  N  =  64  a  single  entry  is  used  to 
identify  a  matrix  containing  the  64  CID's.  A 
special  partitioning  format  is  used  to  input  fC] 
for  user  convenience,  as  shown  on  MTOX-CID 
directive. 


1.  The  WTRX-OID  directive  defines  the  input 
sequence. 
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Input  Pi rective: 
Description: 

Format: 

Example  Syntax : 

Field 

MID 

FRAME 


POINTS 

Cl,  C2, . . • ,  CN 


BULK  DATA  INPUT 


MATE  Materials  compliance  matrix 

Matrix  identification  numbers  for  the  three- 
dimensional  strain-stress  equations  at  N  points  in 
an  element. 

MATE,  MID,  ERAME,  POINTS,  Cl,  C2,...,  CN 
MATE,  7,  2,  8,  6  THRU  12,  16 


Contents 

Material  identification  number. 

=  1,  E  =  C  1  Matrix  is  in  an  orthonormal 
Cartesian  frame. 

=  2,  E  =  C-1  Matrix  relates  contravariant  physical 
components  of  strain  to  covariant 
physical  components  of  stress  in  the 
parametric  frame.  (Useful  with 
orthotropic  bodies  of  revolution.) 

=  l,  Constant  material  properties  (N  =  1). 

=  8,  Trilinear  variation  of  material  properties 

(N  =  8). 

=  64,  Tricubic  variation  of  material  properties 
(N  =64). 

Matrix  identification  numbers  for  the  strain- 
stress  coefficient  matrices  at  the  interpolation 
points.  If  N  =  64  a  single  entry  is  used  to 
identify  a  matrix  containing  the  64  CID's.  A 
special  partitioning  format  is  used  to  input  [E] 
for  user  convenience,  as  shown  on  MTOX-CID 
directive. 


98 


BULK  DATA  OPTIONS 


PAGE  4-61 


Input  Directive; 
Description: 

Format: 

Example  Syntax : 

Field 

MID 

FRAME 

POINTS 

01 ,  02, ,  ON 


Remarks: 


BULK  DATA  INPUT 


MATOR  Orthotropic  material  definition 

Defines  the  material  properties  for  a  linear, 
temperature  independent,  orthotropic  material  from 
engineering  constants. 

MATOR,  MID,  FRAME,  POINTS,  01,  02,...,  ON 
MATOR,  3,  1,  8,  2,  6  THRU  12 


Contents 


Material  identification  number. 

=  1,  Properties  are  in  an  orthonormal  Cartesian 
frame. 

=  2,  Properties  are  in  the  normalized  parametric 
frame  of  the  element.  Assumes  parametric 
frame  is  quasi-cylindrical  or  spherical. 

=  1,  Constant  material  properties  (N  =  1) . 

=  8,  Trilinear  variation  of  material  properties 

(N  =  8). 

=  64,  Tricubic  variation  of  material  properties 
{N  =  64) . 

Matrix  identification  number  for  the  orthotropic 
material  engineering  constants  at  the  interpola¬ 
tion  points.  If  POINTS  is  equal  to  64  a  single 
eptry  is  used  to  identify  a  matrix  containing  the 
64  OID's. 


1.  The  engineering  constants  are  entered  in 
sequence  Ell,  E22,  E33,  NU12,  NU13,  NU23,  G12, 
G13,  G23  on  the  MTOX  matrix  card(s). 

2.  The  Air  Force  Design  Guide  convention  for 
Poisson  ratios  is  used; 

l.e.,  E..  NU^  =  E^  NUir 
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Input  Directive: 
Description: 

Format : 

Example  Syntax : 

Field 

MID 

MPID 

DTI 

Remarks: 


BULK  DATA  INPUT 


MATTA  Temperature-dependent  properties 

Specifies  the  data  table  functions  that  define 
the  temperature  dependence  of  the  thermal 
expansion  moduli  relative  to  a  reference  test 
temperature. 

MATTA,  MID,  MPID,  DTI,  DT2,  DT3,  DT2,  OT5,  DT6 
MATTA,  3,,  1,  1,  1,  2,  2,  2 


Contents 


Material  identif ication  number.  The  same  MID  must 
appear  on  a  MATAL  card  if  normalized  data  table 
functions  are  used. 

Identification  number  of  the  matrix  that  defines 
the  meshpoints  for  which  these  data  apply.  The 
default  is  all  meshpoints  in  the  element.  TTiis 
may  be  delimited  by  other  MATTA  cards. 

Data  table  function,  OTCS  or  ETPC,  for  component  I 
of  the  a  vector.  The  default  for  a  blank  field 
is  the  or,  given  by  the  MATAL  card  if  present  and 
zero  if  none  is  defined. 


1.  When  the  reference  temperature  for  a  data 
table  function  is  not  equal  to  the  ambient 
temperature  for  the  case,  then  u(T)  is 
computed  as 

a(T)  =  (f(T)(T-Tref)-f(TA)(TA-Tref))  /  (T-T A) 

to  ensure  zero  thermal  strain  at  ambient 
temperature. 
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BULK  DATA  INPUT 


Input  Directive: 
Description: 


Format : 

Example  Syntax : 

Field 

MID 

MPID 


DTI 


Remarks : 


MATTAT  Temperature-dependent  properties 

Specifies  the  data  table  functions  that  define  the 
temperature  dependence  of  the  thermal  expansion 
moduli  relative  to  a  reference  test  temperature 
using  strain  data  directly. 

MATTAT,  MID,  MPID,  ET1,  DT2,  OT3,  OT4,  DT5,  ET6 
5,,  11,  14,  3 


Contents 


Material  identification  number.  The  same  MID  may 
appear  on  a  MATTA,  but  the  meshpoints  referenced 
by  MPID  must  be  different. 

Identification  number  of  the  matrix  that  defines 
the  meshpoints  for  which  these  data  apply.  The 
default  for  a  blank  field  is  all  meshpoints  in  the 
element.  This  may  be  delimited  by  other  MATTAT 
cards. 

Data  table  function,  DTCS  or  DTPC,  for  component  I 
of  the  strain  vector  t.  The  default  for  a  blank 
field  is  zero. 


1.  When  the  reference  temperature  for  a  data 
table  function  is  not  equal  to  the  ambient 
temperature  for  the  case,  then  a(T)  is  com¬ 
puted  as 

a(T)  =  (f  (T)  -  f  (Ta)  )  /  (T  -  Ta) 

to  ensure  zero  thermal  strain  at  ambient 
temperature. 

2.  Since  there  is  no  temperature- independent 
property  card  based  on  thermal  strain  data, 
the  data  table  functions  cannot  be  normalized. 
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Input  Directive: 
Description: 

Format : 

Example  Syntax: 

Field 

MID 

MPID 

DTIJ 

Remarks: 


BULK  OATA  INPUT 


MATTC  Temperature-dependent  stiffness  properties 

Specifies  the  data  table  functions  that  define 
the  temperature  dependence  of  each  CIJ  component. 

MATTC,  MID,  MPID,  OT 11,  DT12,  DT13,  VT?2,  OT23, 
DT33 

MATTC,  7,,  10,  11,  11,  10,  11,  10 


Contents 

Material  identification  number.  The  same  MID  must 
appear  on  a  MATC  card  vhen  using  normalized  data 
table  functions. 

Identification  number  of  the  matrix  that  defines 
the  meshpoints  for  which  these  data  apply.  The 
default  for  a  blank  field  is  all  meshpoints  in  the 
element.  This  may  be  delimited  by  other  MATTC 
cards. 

Identification  number  of  the  data  table  function 
for  component  CIJ  in  the  sequence  defined  on  the 
MTRX-C ID  card.  The  default  for  a  blank  field  is 
the  CIJ  given  by  the  MATC  card  if  present  and 
zero  if  none  is  defined. 


1.  Normalized  data  table  functions  result  in 
CIJ  (T)  -  CIJ  f*  (T) 

where  CIJ  is  defined  by  a  MATC  card.  This 
allows  one  data  tablp  function  to  model 
several  CIJ  (T)  when  components  have  similar 
temperature  dependence. 
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Input  Pi  recti ve: 
Description: 

Format : 

Example  Syntax: 

Field 

MID 

MPID 

rm  j 

Remarks : 


BULK  DATA  INPITT 


MATTE  Temperature-dependent  compliance  properti.-; 

Specifies  the  data  table  functions  that  define 
the  temperature  dependence  of  each  component, 

EIJ  =  CIJ  inverse. 

MATTE,  MID,  MPID,  DT11 ,  DIL  2,  I7T13,  DT22,  OT23, 

Dr  33 

MATTE,  5,  4,  11,  12,  12,  11,  12,  13 


Contents 

Material  identification  number.  The  same  MID  must 
appear  on  a  MATE  card  when  using  normalized  data 
table  functions. 

Identif ication  number  of  the  matrix  that  defines 
the  meshpoints  for  which  these  data  apply.  The 
default  for  a  blank  field  is  all  meshpoints  in  the 
element.  This  may  be  delimited  by  other  MATE 
cards. 

Identification  number  of  the  data  table  function 
for  component  DIJ  in  the  sequence  defined  on  the 
MTRX-CID  card.  The  default  for  a  blank  field  is 
the  EIJ  given  by  the  MATE  card  if  present  and  zero 
if  none  is  defined. 


1.  Normalized  data  table  functions  result  in 
EIJ  (T)  =  DIJ  f*  (T) 

where  EIJ  is  defined  by  a  MATE  card.  This 
allows  one  data  table  function  to  model 
several  EIJ  (T)  when  components  have  similar 
temperature  dependence. 


103 


BULK  DATA  OPTIONS 

Input  Directive: 
Description: 

Format : 

Example  Syntax: 

Field 

MID 

MPID 

DTI 

Remarks : 


PAGE 


BULK  DATA  INPUT 


MATTO  Temperature  dependent  orthotropic 
properties 

Specifies  the  data  table  functions  that  define  the 
temperature  dependence  of  the  nine  elastic 
constants  for  an  orthotropic  material. 

MATTO,  MID,  MPID,  DTI,  DTP,  I7T3,...,  I7T9 

MATTO,  1,,  1,  2,  3,  4,  S,  S,  6,  7,  3 


Contents 

Material  identification  number.  The  same  MID 
must  appear  on  a  MATOR  directive  when  using 
normalized  data  table  functions. 

Identification  number  of  the  matrix  that  defines 
the  mesh  points  for  which  these  data  apply. 

The  default  for  a  blank  field  is  all  mesh  points 
in  the  element (s)  referencing  this  MID.  The  list 
may  be  delimited  by  other  MATTO  directives. 

Identification  number  of  the  data  table  function 
for  Ell,  E22,  E33,  NU12,  NU13,  NU23,  G12,  G13, 
G23  in  that  order.  The  default  for  a  blank 
field  is  the  value  given  that  field  by  the  MATOR 
if  present  and  zero  if  none  is  defined. 


1.  Normalized  data  table  functions  result  in 
Eli (T)  =  Eli  •  f*  (T) 

where  Ell  is  defined  by  a  MATOR  directive.  This 
feature  allows  one  data  table  function  to  model 
several  EII(T)  when  several  modclii  have  the  same 
temperature  dependence. 
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Input  Directive: 
Description: 

Format : 

Example  Syntax: 

Field 

MID 

E 

G 

NU 

Remarks: 
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BULK  DATA  INPUT 


MAT1  Material  property  definition 

Defines  the  material  properties  for  linear, 
temperature- independent,  isotropic,  and  certain 
orthotropic  materials. 

MAT1,  MID,  E,  G,  NU 

MAT1,  17,  3. +7,  1.9+7 

Contents 

Material  identification  number  (integer,  1  to  50). 
Young's  Modulus  (real,  >  0  or  blank). 

Shear  Modulus  (real,  >  0  or  blank). 

Poisson's  Ratio  (real  or  blank). 


1.  If  all  three  elastic  constants  are  input 
the  material  is,  in  general,  orthotropic. 

2.  Any  two  of  the  elastic  constants  can  be  input 
to  define  an  isotropic  material.  TTie  blank 
field  value  will  be  automatically  computed. 
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Input  Directive: 
Description: 

Format : 

Example  Syntax : 

Field 

MID 

MPID 

DTI 


BULK  DATA  INPUT 


MATT1  Temperature  dependent  properties 

Specifies  the  data  table  functions  that  define 
the  temperature  dependence  for  an  isotropic 
material . 

MATT1 ,  MID,  MPID,  OT1,  OT2,  DT3 
MATT1,  3,,  1,,  1 


Contents 


Material  identification  number.  The  same  MID  must 
appear  on  a  WT1  card  when  using  normalized  data 
table  functions. 

Identification  number  of  the  matrix  that  defines 
the  meshpoints  for  which  these  data  apply.  The 
default  for  a  blank  field  is  all  meshpoints  in  the 
element.  This  may  be  delimited  by  other  MATT] 
cards. 

Identification  number  of  the  data  table  function 
for  E,  G,  and  NU,  respectively.  The  default  for  a 
blank  field  is  the  moduli  given  by  the  MAT1  card 
when  present.  If  no  MAT1  card  is  present,  at 
least  two  data  table  functions  must  be  present. 
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Input  Directive: 
Description: 

Format : 

Example  Syntax : 

Field 

SID 

COMPS 

EID-P 

IJK-P 

EID-S 

IJK-5 

PSI,  THETA,  FHI-P 

PSI,  THETA,  EHI-S 


BULK  DATA  INPUT 


MPE1  Mesh  point  equality  constraint 

Imposes  displacement  equality  constraints 
at  one  or  more  mesh  points  in  Cartesian 
coordinates. 

MPE1,  SID,  COMPS,  EID-P,  IJK-P,  EID-S,  IJK-5, 
PSI-P,  THETA-P,  FHI-P,  PSI-S,  THETA-S,  PHI-S 

MPE1,  10,  23,  1,  111,  2,  114,  -2.2,  90,  90, 
-32.2,  90,  90 


Contents 


Constraint  set  identification  number. 

Constrained  displacement  components  Ul,  U2,  U* 
identified  by  any  combination  of  the  digits  1, 
2,  and  3. 

Element  identification  number  associated  with 
the  primary  or  independent  mesh  points. 

Defines  the  primary  mesh  points  using  the  mesh 
point  convention  described  in  the  Commentary 
for  this  directive. 

Element  identification  number  associated  with 
the  secondary  or  dependent  mesh  points. 

Defines  the  secondary  mesh  points  using  the 
mesh  point  convention  described  in  the 
Commentary  for  this  directive. 

Euler  angles  in  the  3,  1,  3  sequence  that 
rotate  the  reference  Cartesian  coordinate 
frame  to  the  constraint  frame  at  a  primary 
mesh  point. 

Euler  angles  in  the  3,  1,  3  sequence  that 
rotate  the  reference  Cartesian  coordinate 
frame  to  the  constraint  frame  at  a  secondary 
mesh  point. 
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Input  Pi  recti ve: 


Remarks: 


MPE1  Mesh  point  equality  constraint 
(Continued) 


1.  The  primary  and  secondary  mesh  points  can 
come  from  the  same  element. 

2.  Conflicts  between  SDC  and  MPE  constraints 
are  not  allowed. 

3.  The  number  of  primary  and  secondary  mesh 
points  must  be  equal  unless  there  is  only 
one  primary. 

4.  Constraint  frame  synthesis  continues  just 
as  with  SPC  constraints.  The  same  mesh 
point  may  have  both  SPC  and  MPE  con¬ 
straints  if  they  are  linearly  independent. 

5.  Equality  constraints  can  be  used  for 
periodic  boundary  conditions  and  slidinq 
interfaces. 
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Commentary: 


The  constraint  of  entire  surfaces  or  strings  of  mesh  points  with  a 
single  MPE  directive  is  provided  as  a  user  convenience.  To  support  this 
option,  the  following  convention  for  IJK  definition  of  mesh  points  is 
adopted . 


IJK 

ELEMENT  SURFACE 

MESH 

POINTS 

1 

Z(0,  i,,  £,) 

111, 

121, 

131,... 

,  134, 

144 

2 

Z(l,  Co,  C3) 

411, 

421, 

431,... 

,  434, 

444 

3 

z(clf  0,  ;3) 

Z(fM,  1,  C3) 

Z(Ct,  C2,  0) 

Z(C3,  C2,  1) 

111, 

211, 

311,... 

,  314, 

414 

4 

141, 

241, 

341,... 

,  344, 

444 

5 

111, 

211, 

311,... 

,  341, 

441 

6 

114, 

214, 

314,... 

,  344, 

444 

11 

Z (1/3,  f2,  C3) 

211, 

221, 

231,... 

,  234, 

244 

12 

z (2/3,  e„  f3) 
Z(Cj,  1/3,  C3) 
ZtCj,  2/3,  f3) 
ztclf  c2,  i/3) 
Z(f.w  C?,  2/3) 

311, 

321, 

331,... 

,  334, 

344 

13 

121, 

221, 

321,... 

,  324, 

424 

14 

131, 

231, 

331,... 

,  334, 

434 

15 

112, 

212, 

312,... 

,  342, 

442 

16 

113, 

213, 

313,... 

,  343, 

443 

111 

N 

O 

O 

o 

111 

211 

z (1/3,  0,  0) 

211 

444 

Z(l,  1,  1) 

444 

1001 

List  dependent 

MTRX-1 

Input  mesh 

point  list 

1002 

List  dependent 

MTRX-2 

Input  mesh 

point  list 

10ID 


List  dependent 


f*rTRX-ID  Input  mesh  point  list 
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The  IJK  convention  for  individual  mesh  points  is  identical  to  the 
point  format  input  convention.  The  convention  for  surfaces  is 
illustrated  in  Figure  4-2  for  the  six  external  surfaces  or  faces  of  an 
element  and  the  six  internal  parametric  surfaces.  The  use  of  an  MTOX 
directive  to  input  a  list  of  IJK  mesh  points  is  convenient  when  several 
elements  are  similarly  constrained. 
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BULK  DATA  INPUT 


Input  Directive: 
Description: 


MPE2  Mesh  point  equality  constraint 

Impsses  displacement  equality  constraints  at 
one  or  more  mesh  points  in  local  surface 
coordinates. 


Format : 


MPE2,  SID,  COMPS,  EID-P,  IJK-P,  EID-S,  IJK-5, 
PSI-P,  TOETA-P,  PHI-P,  PSI-S,  THETA -S,  PHI-S 


Example 


MPE2,  10,  123,  1,  5,  2,  6 


Field 


Contents 


COMPS 


EID-P 


IJK-P 


Constraint  set  identification  number. 

Constrained  displacement  components  Ul,  U2,  U3 
identified  by  any  combination  of  the  digits  1, 
2,  and  3  where  u  =  Ultj  +  U2t?  +  U3m. 

Element  identification  number  associated  with 
the  primary  or  independent  mesh  points. 

Defines  the  primary  mesh  points  using  the  mesh 
point  convention  described  in  the  Commentary 
for  the  MPE1  directive. 


f.id-s 


PSI,  THETA,  PHI-P 


PSI,  THETA,  FHI-S 


Element  identification  number  associated  with 
the  secondary  or  dependent  mesh  points. 

Defines  the  secondary  mesh  points  using  the 
mesh  point  convention  described  in  the 
Comnentary  for  the  MPE1  directive. 

Euler  angles  in  the  3,  1,  3  sequence  that 
rotate  the  reference  Cartesian  coordinate 
frame  to  the  constraint  frame  at  a  primary 
mesh  point. 

Euler  angles  in  the  3,  1,  3  sequence  that 
rotate  the  reference  Cartesian  coordinate 
frame  to  the  constraint  frame  at  a  secondary 
mesh  point. 
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Input  Directive: 


Remarks: 


MPE2  Mesh  point  equality  constraint 
(Continued) 


1.  Hie  primary  and  secondary  mesh  points  can 
come  from  the  same  element. 

2.  Conflicts  between  SDC  and  MPE  constraints 
are  not  allowed. 

3.  'ftie  number  of  primary  and  secondary  mesh 
points  must  be  equal  unless  there  is  only 
one  primary. 

4.  Constraint  frame  synthesis  continues  just 
as  with  SPC  constraints.  The  same  mesh 
point  may  have  both  SPC  and  MPE  con¬ 
straints  if  they  are  linearly  independent. 

5.  Equality  constraints  can  be  used  for 
periodic  boundary  conditions  and  sliding 
interfaces.  The  example  syntax  is  from 

a  run  using  periodic  boundary  conditions. 
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Input  Directive: 
Description: 

Format: 

Example  Syntax : 

Field 

MTRX-ID 

V. 

] 

Remarks: 


BULK  DATA  INPUT 

MTTRX-ID  Input  of  a  matrix  or  table 

Input  of  a  matrix  or  table  of  values  which  is 
referenced  by  1  or  more  bulk  data  cards. 

MTRX-ID,  VI,  V2,  V3,  V4 ,  V5,  V6,  V7,  V8 

MTRX-5,  1.0,  2.0,  3.0,  4.0,  1.0,  2.0,  3.0,  4.0 
MTRX-6,  2  (1T4) 


Contents 


The  NTIRX  identification  number.  The  ID  nunber  is 
a  free-form  integer  within  the  last  4  positions 
of  the  mnemonic  field.  The  dash  (or  any  other 
character)  is  optional. 

Input  values  in  floating  point  format.  Up  to  64 
values  may  be  defined. 


1.  Multiple  values  may  be  indicated  by 

parentheses;  e.g.,  4(20,  5,  8).  Parentheses 
may  not  be  nested;  i.e.,  4(2,  3(4),  3,  5). 
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Input  Directive: 
Description: 

Format : 

Example  Syntax: 

Field 

CIJ 

Remarks: 


BULK  DATA  INPUT 

NTTRX-C  ID 

Stress-strain  coefficient  matrix  associated  with 
CID. 

MTRX-CID,  Cll,  C12,  C13,  C22,  C23,  C33,  C44,  C45, 
C46,  C55,  C56,  C66 

MTRX-6,  1.0E7,  0.3E7,  0.3E7,  2.0E7,  0.2E7,  1.0E7, 
0.0,  0.0,  0.5E7,  0.0,  0.5E7 

Contents 


Cll, 

C12, 

C13, 

C22, 

C23, 

C33 

(Data 

items 

1-6) 

C44, 

C45, 

C46, 

C55, 

C56, 

C66 

(Data 

items 

7-12) 

C14 , 
C34, 

C15, 

C35, 

C16, 

C36 

C24, 

C25, 

C26, 

(Data 

items 

13-21) 

1.  If  the  material  has  cubic  symmetry,  locations 
13-21  will  be  zero  and  need  not  be  entered. 


"oil" 

"Cll 

C12 

C13 

C14 

C15 

C16 

£11 

o22 

C21 

C22 

C23 

C24 

C25 

C26 

£22 

o33 

C31 

C32 

C33 

C34 

C35 

C36 

£33 

ol2 

C41 

C42 

C43 

C44 

C45 

C46 

£12 

al3 

C51 

C52 

C53 

C54 

C55 

C56 

£13 

o23 

C61 

C62 

C63 

C64 

C65 

C66 

£23 

—  — 

— 

_  — 

BULK  DATA  OPTIONS 


PAGE  4-78 


BULK  DATA  INPUT 


Input  Directive:  PARAM  Special  parameters 

Description:  Special  program  parameter  definition  card  for 

flags  and  other  constants. 

Format:  PARAM,  VARl,  VALUE1,  VAR2,  VALUE2, . . . ,  VARN, 

VALUEN 

Example  Syntax :  PARAM,  ACCG,  6,  DEBUG,  1,  DISK,  13 


Field  Contents 

VAR  i  Name  of  the  variable  to  be  defined. 

VALUE  i  Value  assigned  to  VAR  i. 


Currently  the  following  parameters  are  recognized: 

ACCG  Number  of  significant  digits  of  accuracy  to  be 

used  in  the  test  for  consistency  between  generate 
grid  point  coordinates.  Default  value  =  5. 

ORDERG/ORDERD  Changes  the  geometric  (G)  or  data  (D)  line  inter¬ 

polation  order  to  cubic  first;  i.e.,  if  the  same 
grid  point  appears  on  two  line  cards  without  input 
values,  the  (first)  cubic  line  will  be  used  to 
interpolate  the  point  rather  than  the  (first) 
linear  line  for  VALUE  =  1.  Default  value  =  0. 

DEBUG  Debugging  flag  to  dump  certain  data  during 

execution.  Default  value  =  0.  See  Chapter  9. 

AXY  The  subtended  angle  in  a  OCX  axisymmetric  model. 

(Continued  on  following  page.) 
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Input  Directive: 

Field 

ITER 

NSTP 

ERRORS 

FLLINK 

CONVRG 

EKREL 


PARAM  Special  parameters  (continued) 


Contents 


Maximum  number  of  iterations  to  be  executed  in  the 
conjugate  gradient  solution.  Default  =  110%  of 
NFSET. 

Maximum  number  of  steepest  descent  moves  before  a 
CG  solution  is  terminated.  Default  =  5. 

Maximixn  number  of  errors  that  the  system  will 
permit  beyond  the  DRY  point.  Default  =  0. 

Field  length  override.  Value  is  n,  m  vhere  n  is 
the  link  number  and  m  is  the  octal  number  of 
thousand  words  of  core;  e.g.,  12,  16  requests  link 
12  to  begin  with  16K  octal  words. 

CG  solution  convergence  value.  Convergence  is 
said  to  have  occurred  vrfaen  the  difference  term  is 
unchanged  to  CONVRG  significant  places. 

Default  =  8. 

Relative  factor  in  element  stiffness  generation 
which  defines  how  small  a  number  in  algebraic 
format  must  be  relative  to  the  largest,  value  in 
the  element  before  it  is  set  to  0.0. 

Default  =  l.E  -  11  on  UNIVAC,  1.3  -  9  on  CDC. 
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} 

>  BULK  DATA  OPTIONS 


Input  Pi rective: 

Description: 

Format: 

Example  Syntax : 

Field 

PID 

FORMAT 

MTRX-1,  2,  3 
TID 

Remarks: 
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BULK  DATA  INPUT 


PATCH  Patch  generator 

Direct  input  of  a  patch  in  any  format. 

PATCH,  PID,  FORMAT,  MTOX-1 ,  MTRX-2,  MTRX-3,  TID 
PATCH,  4,  S,  101,  102,  103 


Contents 


Patch  identification  nunber. 

A  or  S  Algebraic  coefficients. 

B  Geometric  coefficients. 

P  Point  coefficients. 

G  Gaussian  coefficients. 

Matrix  identification  number  for  Zl,  Z2,  and  Z3. 

Transformation  identification  number,  if  any,  to 
be  applied  to  the  patch. 


1.  The  default  FORMAT  is  point  format. 

2.  The  16  coefficients  input  in  any  format 

for  each  patch  coordinate  function  are  input 
in  the  sequence  Pll,  P21,  P31,  P41,  P21, 

P22, . . . ,  P44. 
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BULK  DATA  IN  PITT 


Input  Directive: 
Description: 


Format : 


Example 


PATCHGR  Patch  generated  by  general  line  rotation 

Generates  a  bicubic  patch  for  the  surface  created 
by  rotating  a  PC  line  about  a  general  axis  of 
rotation  through  ganma  degrees. 

PATCHGR,  ID,  LID/SEG,  ZA1,  ZA2,  ZA3,  ZB1,  ZB2, 
ZB3,  TID,  GAMMA,  GAMMO 

PATCHGR,  5,  3,  1.5,  0.,  -3.0,,  2.5,  0.3,,  25 


Field 

ID 

LID,  SEG 


ZAI ,  ZB I 

TID 


GAMMA,  GAMMO 


Contents 


The  identification  number  to  be  given  the  patch 
generated  from  line  LID,  segment  nunber  SEG. 

The  line  number,  LID,  and  segment  number,  SEG, 
that  identify  the  PC  line  to  be  rotated.  A  blank 
SEG  defaults  to  one. 

Coordinates  of  two  points  that  define  the  rotation 
axis  directed  from  ZA  to  ZB. 

Transformation  ID,  if  any,  that  defines  a  geo¬ 
metric  transformation  to  be  applied  to  the  PC 
line  before  rotation.  The  line,  LID,  does  not 
change . 

The  PC  line  is  rotated  through  GAMMA  degrees 
starting  GAMMO  degrees  from  the  initial  position 
of  the  line.  The  sense  of  rotation  is  determined 
by  the  right-hand  rule  and  the  directed  line 
(vector)  from  ZA  to  ZB. 
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Commentary: 

Generation  of  a  surface  patch  by  rotating  a  PC  line  about  a  general 
axis  of  rotation.  Four-sided  and  degenerate  three-sided  patches  are 
constructed. 


Method : 

Given  the  PC  line  Z(\),  rotate  the  points  Z(0),  Z(l/3),  Z (2/3 )  r  and 
Z(l)  to  form  the  surfact  patch  using  LINEGR.  Hie  axis  of  rotation  is 
the  vector  from  ZA  to  Z^. 


120 


BULK  DATA  OPTIONS 

Input  Directive: 
Description: 

Format : 

Example  Syntax : 

Field 

ID 

Gl,  G2,  G3,  G4 
LI,  L2,  L3,  L4 

Remarks: 
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BULK  DATA  INPUT 


PATCHL  Define  a  patch  with  boundary  lines 

Generates  a  bicubic  patch  from  four  grid 
line  segments.  Hie  orientation  of  the  segments 
is  checked,  as  is  the  presence  of  referenced 
grid  points. 

PATCHL,  ID,  Gl,  G2,  G3,  G4,,,,,  LI,  L2,  L3,  L4 
PATCHL,  2,  1,  2,  3,  4,,,,,  4,  5,  8,  12 


Contents 


The  identification  to  be  given  to  the  patch. 

Grid  point  identification  nunbers  of  the  corner 
points  of  the  patch. 

Line  identification  numbers  for  the  lines  that 
bound  the  patch. 


1.  Hie  grid  point  and  line  sequencing  is  shown 
below. 
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Input  Directive: 
Description: 

Format : 

Example  Syntax : 

Field 

ID1 

BLID,  SEG 

OLID,  SEG 

TID 

FRAME 

ID2,  3,...,  N 


BULK  DATA  INPUT 


PATCHO  Outline  patch (es) 

Generates  patch (es)  by  moving  an  outline  curve 
along  a  base  curve  with  a  fixed  orientation  of  the 
outline  curve  in  the  global  frame  or  in  the  local 
Frenet  frame  of  the  base  curve. 

PATCHO,  ID1,  BLID/SEG,  OLID/SEG,  TID,  FRAME, 

ID2, . . . ,  IDN 

PATCHO,  5,  6/1,  3,,  F,  11  THRU  15 


Contents 


Patch  identification  number  for  first  patch 
generated . 

Baseline  identification  number  and  segment  number, 
SEG.  If  the  SEG  is  not  specified,  the  entire  line 
BLID  is  used. 

Outline  curve  identification  number  and  segment 
number,  SEG.  If  the  SEG  is  not  specified,  the 
entire  line  OLID  is  used. 

Transformation  ID,  if  any,  that  defines  a  rotation 
matrix  to  reorient  the  outline  curve  relative  to 
the  base  curve.  The  outline  curve  is  always 
translated  to  the  first  grid  point  of  the  baseline 
independent  of  TID. 

E  for  fixed  outline  orientation  with  respect  to 
the  Cartesian  frame.  F  for  fixed  outline  orienta¬ 
tion  with  respect  to  the  local  Frenet  or  binormal 
coordinate  frame  of  the  base  curve. 

List  of  identification  numbers  to  be  given  the 
second  and  subsequent  patches  generated,  if  any. 
This  sequence  proceeds  from  the  second  line  seg¬ 
ment  of  OLID  to  the  last  and  then  repeats  from 
the  first  segment  for  the  next  segment  of  BLID. 
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Input  Directive: 
Remarks : 


PAGE 

PATCHO  Outline  patch (es)  (Continued) 

1.  The  Frenet  frame  is  defined  by 

T  =  Z,  ^  '  y=?'55'  B=TXN 


-85 
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PAGE  4-86 

BULK  DATA  INPUT 

Input  Directive: 

PATCHQ  Bilinear  patch  from  four  corner  points 

Descr ipt ion : 

Defines  a  parametric  bilinear  patch  from  four 
previously  defined  corner  grid  points. 

Format : 

PATCHQ,  ID,  Gl,  G2,  G3,  G4 

Example  Syntax: 

PATCHQ,  5,  1,  2,  7,  6 

„  leld 

Contents 

ID 

Patch  identification  number  (integer,  1  to  50). 

G 1 1  G2,  G3,  G4 

Grid  point  identification  numbers  of  the  four 
corner  points  in  the  order  (0,  0),  (0,  1),  (1,  1), 
(1,  0)  (integer,  1  to  9999). 
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Input  Directive: 
Description: 

Format: 

Example  Syntax : 

Field 

ID1 

LID 

SEG 

7.1,  7.2,  Z3 
1 1  iETAB ,  THETAE 

IDRA 

ID2,  ID3, . . . ,  IDN 

Remarks: 


BULK  DATA  INPUT 


PATCHR  Patch  generation  by  1 ine  rotation 

Generates  a  bicubic  surface  of  revolution  patch 
or  set  of  patches  by  the  rotation  of  a  line  seg¬ 
ment  or  set  of  line  segments  through  a  specified 
arc  about  a  coordinate  axis. 

PATCHR,  ID1,  LID/SEG,  Zl,  Z2,  Z3,  THETAB,  TOETAE, 
IDRA,  ID2,  ID3,  ID4,...,  IDN 

PATCHR,  3,  5/2,  0.0,  5.0,  -12.2,  30.,  60.,  -3 


Contents 


The  identification  nimber  to  be  given  to  the  first 
patch  to  be  generated.  This  is  the  patch  genera¬ 
ted  from  line  LID,  segment  SEG  if  SEG  is  defined 
and  segment  1  if  SEG  is  not  defined  (integer,  1  to 
50). 

Tfae  identification  nunber  of  the  line  to  be  rota¬ 
ted. 

If  specified,  the  segment  within  line  LID  to  be 
rotated.  If  not  specified,  all  segments  of  line 
LID  will  be  rotated. 

Coordinates  for  shifting  the  origin. 

The  beginning  and  ending  angles  through  which  the 
segment  or  segments  will  be  rotated  (degrees) . 

Axis  of  rotation  (integer,  -3  to  3). 

+1  =  +X,  -1  =  X,  +2  =  +Y,  +3  =  +Z,  etc. 

List  of  identification  nunbers  to  be  given  to  the 
second,  third,...,  nth  patch  to  be  generated. 

The  list  may  include  any  of  the  list  operators 
such  as  THRU,  EXCEPT,  etc. 


1.  See  the  PATCHGR  commentary  for  a  description 
of  the  method. 
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Input  Directive: 
Description : 

Format : 

Example  Syntax : 

Field 

PID 

LI,  2,  3,  4 

SEG 

TID 


BULK  DATA  INPUT 
PATCH4L  Patch  generator 

Generates  a  patch  from  four  lines  in  parametric 
direction  tv*>. 

PATCH4L,  PID,  Ll/SEG ,  L2/SEG,  L3/SEG,  L4/SEG,  TID 
PATCH4L,  5,  1/3,  17/6,  9,  12 

Contents 

Patch  identification  number. 

Line  identification  number  for  lines  Z(£  ,  0), 
Z(5j,  1/3),  Z(Clf  2/3),  and  Z(Cj,  1).' 

The  segment  number.  (Default  SEG  is  1.) 

Transformation  identification  number,  if  any,  to 
be  applied  to  the  patch. 
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Input  Directive: 
Description: 

Format : 

Example  Syntax : 

Field 

ID 

EID 

DPID 

Scalar 

Remarks : 


BULK  DATA  INPUT 


PLOAD3  Pressure  load  data 

Defines  a  normal  pressure  load  on  one  surface 
of  a  3-D  element. 

PL0AD3,  ID,  EID,  DPID,  Scalar 

PLQAD3,  6,  2,  10 


Contents 


Load  set  identification  number. 

Element  identification  number. 

Data  patch  identification  number. 

The  magnitude  of  the  data  in  DPID  will  be  scaled 
by  this  number.  Hie  default  value  is  1.0. 


1.  The  grid  points  referenced  by  the  data  patch 
are  used  tc  locate  the  face  of  element  EID  on 
which  the  pressure  acts.  A  positive  pressure 
is  directed  in  the  positive  direction  of  the 
normal  to  the  surface. 

2.  Refer  to  Figure  2-6  for  definition  of  normal 
vectors  on  the  hyperpatch  faces.  WARNING:  A 
positive  pressure  is  not  in  on  all  faces! 
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Input  Directive: 
Description: 

Format : 

Example  Syntax : 

Field 

EID 

MID 

DCID 

PSI,  THETA,  FHI 


Remarks : 


BULK  DATA  INPUT 


PPDE3  Property  data  for  a  solid  element 

Property  card  for  a  three-dimensional  parametric 
discrete  element. 

PPDE3,  EID,  MID,  DCID,  PSI,  THETA,  PHI 
PPDE3,  10,  5,,  30.0,  0.0,  1302 


Contents 


Element  identification  number. 

Material  identification  number. 

Direction  cosine  matrix  identification  minber. 

If  blank  or  zero,  the  Euler  angle  option  will  be 
used  to  relate  the  material  axes  to  the  reference 

axes  e  . 

'  i 

Euler  angles  in  the  3,  1,  3  sequence  that  rotate 
the  reference  axes  to  the  material  axes  relative 
to  the  21,  22,  Z3  system.  Same  sequence  of  Euler 
angles  as  used  by  H.  Goldstein,  Classical 
Mechanics,  p.  107,  1959  edition.  If  an  Euler 
angle  varies  over  the  element,  its  data  hyperpatch 
ID  +  1000  is  entered  in  the  field  normally  used 
for  that  angle. 


1.  If  material  MID  is  given  in  the  parametric 
frame  (TYPE  =  2) ,  any  direction  cosine  or 
Euler  angle  data  will  be  ignored. 

In  the  above  example,  PHI  is  defined  by  data 
hyper patch  302. 
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Input  Directive: 
Description: 

Format : 

Example  Syntax : 

Field 

SI,  S2,  S3 

Z1(0),  Z2 (0) ,  Z3 (0) 
PATCH1,...,  PATCHN 


BULK  DATA  INPUT 


SCALP  Scaling  of  a  bicubic  patch 

Scales  a  bicubic  patch  or  set  of  patches  relative 
to  a  specified  scaling  origin. 

SCALP,  SI,  S2,  S3,  Z1(0),  Z2(0),  Z3(0),,,  PATCH1, 
PATCH2 , . . . ,  PATCHN 

SCALP,  1,  2.0,  1.5,,  5.,,,,  4,  15  THRU  40 
EXCEPT  22 


Contents 


Scaling  factors  in  the  Zl,  Z2,  and  Z3  directions 
(default  =  1.0) . 

Origin  of  scaling  (default  =  0). 

The  list  of  patch  identification  numbers  to  be 
scaled. 


Remarks: 


1.  Hie  origin  strongly  influences  the  scaling 
operation.  See  SCALPH  commentary. 


BULK  DATA  OPTIONS 

1 

1 

BULK  DATA  INPITT 

Input  Directive: 

SCALPH  Scaling  of  all  hyperpatches 

Description: 

Scales  all  hyperpatches  relative  to  a 
specified  scaling  origin. 

Format : 

SCALPH,  SI,  S2,  S3,  Zl(0),  Z2(0),  Z3(0) 

Example  Syntax: 

SCALPH,  1.,  2.0,  1.5,,  5. 

Field 

Contents 

SI,  S2,  S3 

Scaling  factors  in  the  Zl,  Z2,  and  Z3  directions 
(default  =  1.0) . 

7.1(0),  22  (0)  ,  Z3  (0) 

Origin  of  scaling  (default  =  0). 

Remarks: 

1.  The  origin  strongly  influences  the  scaling 
operation.  See  cornnentary. 

1 30 


BULK  DATA  OPTIONS 


PAGE  4-93 


Commentary: 

The  SCALP  and  SCALPH  cards  both  are  based  on  the  transformation 
ZS  =  [S]  (Z  -  ZO)  +  zo 

where  [S]  is  a  diagonal  matrix  of  coordinate  scale  factors,  ZO  is  the 
scaling  origin,  and  ZS  is  the  vector  of  scaled  coordinates.  Note  that 
ZS,r^  =  [S]  Z ,  ^ .  The  effect  of  ZO  on  scaling  is  illustrated  below  for  a 
parallelopiped  where  [S]  is  the  same  for  both  figures. 


8100609 
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Input  Directive: 
Description: 

Format : 

Example  Syntax : 

Field 

ID 

EID 

DPI,  2,  3 

SI,  2,  3 

Remarks: 


BULK  DATA  INPUT 

SDC1  Surface  displacement  constraint 

Imposes  a  surface  displacement  constraint  on 
components  in  the  reference  Cartesian  frame. 

SDC1,  ID,  EID,  DPI,  DP2,  DP3,  SI,  S2,  S3 

SDC1,  3,  16,  2,,  1,  3.4 

Contents 

Constraint  set  identification  number. 

Element  identification  number. 

Data  patches  for  the  components  Ui  of  the 
displacement  u  =  Ule  +  U2e  +  U3e  .  A  blank 

■v  ~  -L  -3 

field  indicates  the  component  is  unconstrained. 

Scale  factors  applied  to  the  data  patches. 

Default  values  are  1.0. 

1.  The  data  patches  must  reference  the 

same  grid  points,  and  these  are  used  to  locate 
the  constrained  surface  on  elem«.-.t  EID. 
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Input  Directive: 
Description: 

Format : 

Example  Syntax: 

Field 

ID 

EID 

COMPS 

Gl,  2,  1,  4 


BULK  DATA  INPUT 


SDC10  Surface  displacement  constraint  —  zero 
components 

Imposes  a  surface  displacement  constraint  of  zero 
on  components  in  the  reference  Cartesian  frame. 

SDC10,  ID,  EID,  COMPS,  Gl,  G2,  G3,  G4 

SDC10,  10,  15,  23,  11,  12,  16,  24 


Contents 


Constraint  set  identification  nimber. 

Element  identification  nimber. 

Components  UI  of  the  displacement  vector 
u  =  Ulej,  U2e?,  U3e3  that  are  to  be  constrained 
to  zero. 

Corner  grid  points  of  element  EID  that  define  the 
constrained  surface  on  that  element. 
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BULK  DATA  INPUT 


Input  Directive: 


Description: 


Format : 


Example  S' 


SDC20  Surface  displacement  constraint  —  zero 
components 

Imposes  a  surface  displacement  constraint  of 

zero  on  components  in  the  local  surface  coordinate 

directions. 

SDC20,  ID,  EID,  COMPS,  Gl,  G2,  G3,  G4 
SDC20,  5,  6,  31,  8,  9,  1,  3 


Fi el d  Contents 

ID  Constraint  set  identification  nunber. 


EID 


Element  identification  nunber. 


COMPS  Components  UI  of  the  displacement  vector 

u  =  Ult  +  U2t  +  U3n  in  the  local  surface  frame. 

~1  -  2 

The  constrained  components  are  identified  by  any 

combination  of  the  digits  1,  2,  and  3.  These  are 

physical  components  in  that  tj,  t?,  n  are 

normalized  to  unity  for  constraint  purposes. 

Gl,  2,  3,  4  Corner  grid  points  of  element  EID  that  define  the 

constrained  surface  on  that  element.  The  patch 
for  that  surface  is  extracted  from  the  hyperpatch 
EID  and  used  to  compute  tj,  t7,  and  r  on  that 
surface. 
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Input  Directive: 
Description : 

Format: 

Example  Syntax : 

Field 

ID 

GID 

UI 

DCFL 


DCID 

PSI,  THETA,  PHI 


BULK  DATA  INPUT 


SPC1  Single  point  constraint 

Imposes  a  displacement  constraint  at  a  grid  point 
for  components  in  an  orthonormal  frame. 

SPC1,  ID,  GID,  Ul,  U2,  U3,  DCFL,,,  DCID,  PSI, 
THETA,  PHI 

SPC1,  3,  5,,  0.0,,  1,,,,  30.0 

Contents 


Constraint  set  identification  number. 

Grid  point  identification  number. 

Value  of  constrained  displacement  component  I. 

A  blank  field  implies  that  component  is  uncon¬ 
strained. 

Direction  cosine  flag.  If  blank  or  any  value 
other  than  1,  the  components  UI  are  in  the 
reference  ej  frame,  and  the  continuation  card  can 

be  omitted.  A  value  of  1  indicates  the  components 
are  in  the  frame  defined  by  the  continuation  card. 

Direction  cosine  matrix  identif ication  number.  If 
blank  or  zero  the  Euler  angles  PSI,  THETA,  and  mi 
used  to  define  the  constraint  frame. 

Euler  angles  in  the  3,  1,  3  sequence  that  rotate 
the  reference  frame  ej  to  the  constraint  frame. 
Same  sequence  of  Euler  angles  as  used  by  H. 
Goldstein,  Classical  Mechanics,  p.  107,  1959 
edition. 
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Input  Directive: 
Description: 

Format : 

Example  Syntax : 

Field 

ID 

EID 

PIJK 

UI 

UCFL 


DC  ID 

Remarks: 


BULK  DATA  INPUT 


SPC2  Single  point  constraint 

Imposes  a  displacement  constraint  at  a  meshpoint 
of  an  element  for  components  in  an  orthonormal 
frame. 

SPC2,  ID,  EID,  PIJK,  Ul,  U2,  U3,  DCFL, ,  DCID, 

PS  I ,  THETA,  PHI 

SPC2,  10,  6,  223,,,  0.0 


Contents 


Constraint  set  identification  nunber. 

Element  identification  number. 

Meshpoint  identification  number  defined  as  the 
subscripts  in  point  format  Pjjk  for  one  of  the  UI 
meshpoints. 

Value  of  constrained  displacement  component  I. 

A  blank  field  implies  that  component  is  uncon¬ 
strained. 

Direction  cosine  flag.  If  blank  or  any  value 
other  than  1,  the  components  UI  are  in  the  e; 

frame,  and  the  continuation  card  can  be  omitted. 

A  value  of  1  indicates  the  components  are  in  the 
frame  defined  by  the  continuation  card. 

Direction  cosine  matrix  identification  number. 

If  blank  or  zero  the  Euler  angles  PSI,  THETA,  and 
PHI  in  the  3,  1,  3  sequence  (see  SPC1  card)  are 
used  to  define  the  constraint  frame. 


1.  If  the  meshpoint  is  on  (a)  surface(s) 

the  point  also  will  be  constrained  on  all 
connected  elements.  The  location  in  space 
of  the  constrained  point  is  available  from 
the  hyperpatch  for  element  EID. 
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PAGE  4-100 

rnput  Directive: 

TEMP  Temperature  data 

Description: 

Defines  the  temperature  over  the  volume  of  an 
element. 

Format : 

TEMP,  ID,  EID,  EHP 

Example  Syntax: 

TEMP,  10,  3,  6 

Field 

Contents 

ID 

Temperature  set  identification 

number. 

EID 

Element  identification  number. 

DHP 

Data  hyperpatch  identification 

number. 

Remarks: 

1.  The  ambient  temperature  is 
control  card  for  each  case. 

defined  by  a  case 
The  temperature 

difference,  6  =  T  -  Ta,  is  used  for  computing 

thermal  loads  and  strains  vrfiere  the  tempera¬ 
ture  T  is  defined  by  a  TEMP  card.  The  default 
value  for  ambient  temperature  is  zero. 
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Input  Directive: 
Description: 

Format: 

Example  Syntax: 

Field 

ID 

ZOI 

DCID 

PSI,  THETA,  PHI 

TI 

Remarks : 


BULK  DATA  INPUT 


TWOVE  Rigid  body  transformation 

Defines  a  transformation  that  moves  objects 
(lines,  patches,  hyperpatches)  as  rigid  bodies. 

TWOVE,  ID,  ZOI,  Z02,  Z03,  DCID,  PSI,  THETA,  PHI, 
Tl,  T2,  T3 

TWOVE,  1,  3.0,  0.,  3.,,  30.,  20.,  -10.,  3.0, 

0.,  4.0 


Contents 


Transformation  identification  number  (1  to  100). 

Defines  an  origin  for  rotation  of  the  object. 

Direction  cosine  matrix  identification  number. 

If  blank  or  zero  the  Euler  angles  PSI,  THETA,  and 
PHI  define  the  rotation  matrix. 

Euler  angles  in  the  3,  1,  3  rotation  sequence. 
Same  sequence  of  Euler  angles  as  used  by 
Goldstein,  Classical  Mechanics,  p.  107,  1959 
edition. 

Defines  a  translation  to  be  applied  after  the 
rotation. 


1.  The  complete  transformation  can  be 
defined  as 

Z*  =  R(Z  -  ZO)  +  T 
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CHAPTER  5 

CASE  CONTROL  OPTIONS 


5.1  Overview 


Case  control  in  PATCHES-III  consists  of  executive  data  that  directs 
execution,  case  definition  data  that  selects  loads  and/or  constraints, 
and  output  requests  that  select  the  data  for  output.  As  with  bulk  data, 
the  format  and  function  of  the  case  control  input  is  by  design  similar 
to  NASTRAN.  A  checkpoint-restart  editor  has  been  provided  for  more 
efficient  multiple  run  applications  and  to  help  protect  the  user  in  the 
event  of  abnormal  termination.  Unlike  NASTRAN,  no  checkpoint  dictionary 
is  required  for  restart.  The  case  control  data  is  checked  for  input 
errors  and  diagnostics  provided.  These  checks  i.  lude  cross-referencing 
the  bulk  data  to  ensure  all  referenced  load  and/or  constraint  sets  are 
present.  The  syntax  of  all  case  control  cards  is  free  form.  There  are 
no  fixed  fields;  the  data  may  be  placed  anywhere  on  the  card.  A  corona 
acts  as  the  delimiter  between  items  in  a  list. 


140 


CASE  CONTROL  OPTIONS  PAGE  5-2 

5.2  Available  Options 

There  are  over  forty  (40)  case  control  options  available.  These 
options  are  summarized  in  Table  5-1  where  they  are  grouped  by  function 
into  three  basic  categories.  Case  control  directives  may  be  input  in 
any  order  within  a  category,  but  the  categories  must  be  input  in  the 
order  shown.  All  output,  if  there  is  to  be  any,  must  be  selected  by  one 
or  more  of  the  case  control  options.  The  fundamental  model  for  the  most 
data  output  by  PATCHES-m  is  the  parametric  cubic.  This  representation 
is  available  for  lines,  surfaces,  and  volunes.  A  parametric  line 
representing  either  geometric  data  or  physical  data  can  be  output  in 
geometric  or  point  format,  as  shown  in  Table  5-2  for  a  one-component 
line.  Surface  patch  formats  are  also  shown  in  Table  5-2,  including  the 
algebraic  format.  Patch  coefficients  in  the  algebraic  format  do  not 
have  a  simple  expression  in  terms  of  coordinate  (data)  functions  or 
their  derivatives.  They  are  simply  printed  as  a  matrix  of  Sij 
coefficients.  Hyperpatches  can  also  be  output  in  geometric  format, 
Table  5-3,  or  point  format.  The  latter  output  is  the  same  as  that  for  a 
patch  except  that  there  are  four  surfaces:  =0,  L3  =  1/3, 

r,  =  2/3,  and  f.  3  -  1.  The  location  in  space  of  a  data  model  can  only 

be  determined  by  reference  to  the  geometry  of  the  associated  finite 
element.  The  standard  output  format  for  hyperpatches  is  point  format. 
Table  5-4,  and  this  is  the  only  format  used  for  finite  element  solution 
data. 
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Table  5-1 

CASE  CONTROL  DIRECTIVES 


execittive 

DATA 

TITLE 

CHKPNT 

AMBIENT 

TIME 

DRY 

RESTART 

CASE  DEFINITION 

LOAD 

SUBCASE 

SUBCOM 

SDC 

SUBTITLE 

AXY 

OUTPUT  REQUESTS 

OUTPUT 

SET 

EVERYTHING 

ALL 

GRID 

LINEB 

PATCHA 

HPB 

VOLUME 

EMESH 

LINEP 

PATCHB 

HPP 

DETJ 

LORDER 

PATCH? 

DATAG 

DLINEB 

DPATB 

DHPB 

I7TLINEB 

ODISP 

OLQAD 

DLINEP 

DPATP 

DHPP 

CTLINEP 

ELEMENT 

EDISP 

ES TRAIN 

ESTRESS 

EI.OAD 

EFORCE 

MSTRAIN 

MSTRESS 

PSTRAIN 

PSTRESS 
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Table  5-2 

LINE  AND  PATCH  OUTPUT 


Line 


Geometric  - 
Point 

Gaussian  - 
Algobraic*- 


Z(0),  Z  (1 ) ,  Z,r  (0),  Z,r(l) 

Z (0) ,  Z (1/3) ,  Z (2/3) ,  Z(l) 

Z(. 069432),  Z(. 330010),  Z(. 669991),  Z(. 930568) 
SI,  S2,  S3,  S4 


Patch 


Geometric  - 

^  (0,0), 

Z(0,1), 

2  r  £2  ^  9 

0), 

z,c2(o,i) 

i  Z  (1,0)  , 

z (1,1)  , 

z,52d, 

0), 

z,  £,2  (1,1) 

o 

o 

Z,£j(0,l), 

Z, C1E2 (®'®^  ' 

Z,£i£2(°,1) 

Z,fj(l,0), 

z,5l (1,1), 

Z^j^d'O), 

Z,£i£2  (1,D 

Point 

z (0,0)  , 

z (0,1/3), 

Z (0,2/3) , 

2(0,1) 

Z (1/3,0)  , 

Z (1/3, 1/3) 

,  Z  (1/3 

,2/3), 

Z (1/3,1) 

Z (2/3,0), 

Z (2/3, 1/3) 

,  Z (2/3 

,2/3), 

Z (2/3,1) 

Z(1,0), 

Z (1,1/3) , 

Z (1,2/3) 

2(1,1) 

Algebraic*- 

Sll,  S12 , 

S13,  S14 

t 

-9/2 

27/2  -27/2 

9/2 

S21 ,  S22, 

S31,  S32, 

S23,  S24 

S33,  S34 

NU  = 

9 

-11/2 

-45/2  18 

9  -9/2 

-9/2 

1 

S41 ,  S42, 

S43,  S44 

1 

• 

0  0 

0 

*  The  algebraic  coefficients  are  related  to  the  point  format 
coefficients  by  N;,  where  S [  =  Nj^Pj  and  Sjj  =  Nj^P^^Nj]. 
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Table  5-3 

HYFERPATCH  OUTPUT  -  GEOMETRIC  FORMAT 


Z(0,0,0) 

Z(0,1,0) 

z,^  (0,0,0) 

z,c  (0,1,0) 

Z(1,0,0) 

Z(1,1,0) 

z,?2 (1,0,0) 

Z,e? (1,1,0) 

Z,tl  (0,0,0) 

Z,5l (0,1,0) 

Z, ^  ^  (0,0,0) 

Z»(T  jfr2  (0,1,0) 

z,Ci (1,0,0) 

z,Cl (1,1,0) 

Z,^^^  (1,0,0) 

z,^3^2 (1,1,0) 

Z  (0,0,1) 

Z(0,1,1) 

Z,c  (0,0,1) 

z,?2  (0,1,1) 

z  (1,0,1) 

z  (1,1,1) 

z,C2 (1,0,1) 

Z,t? (1,1,1) 

z,tl  (0,0,1) 

Z,5l (0,1,1) 

Z'ClC2  (0'0,1) 

z,f>1  (1,0,1) 

Z, Cl  (1,1,1) 

Z'L^2 (1,0,1) 

z'j;1e;2  (1'1'1) 

Z,^  (0,0,0) 

Z,?3(0,1,0) 

^ f  ^2 £  3  ^ ^ 

(0,1,0) 

Z,  ^(1,0,0) 

Z,t  (1,1,0) 

Z'£2£3 

z»^2^3 (1,1,0) 

Z,5  5  (0,0,0) 

Z^lC?C3(0,0,0) 

Z'C1e2C3(0,1'0) 

Z,  ^(1,0,0) 

Z,£j£3 (1,1,0) 

Z'{lW3,W'0) 

Z*tlt2C3(3,1'0) 

z,f  3  (0,0,1) 

Z,5  (0,1,1) 

Z,, ^3(0,0, D 

t’Ms10'1'1* 

Z,  (1,0,1) 

Z,^  (1,1,1) 

Z,  (1,0,1) 

Z'^3(1'1'1) 

z,Ci£3  (0,0,1) 

Z,eiC3(0,1,1> 

Z,^1^2^3 (0,0,1) 

Z,Ci5253(0,1,1) 

Z'r,1f;3  ( 1 , 0 ' 1 ) 

z'e;1C3(1'1'1) 

Z,ci€ze3(1,0,1) 

Z'5i5253(1'1'1) 
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Table  5-4 

HYPERPATCH  OUTPUT  -  POINT  FORMAT 


[PI] 


[P2] 


[P3] 


[P4] 


2(0, 0,0) , 

Z(0, 1/3,0), 

Z (1/3, 0,0) , 

Z(l/3, 1/3,0), 

Z (2/3, 0,0)  , 

Z (2/3, 1/3,0), 

Z(1,0,0), 

Z (1, 1/3,0), 

z (0,0, 1/3) , 

Z (0,1/3, 1/3) , 

2(1/3, 0,1/3), 

Z (1/3, 1/3, 1/3), 

Z(2/3, 0,1/3), 

Z (2/3, 1/3, 1/3), 

Z(l, 0,1/3), 

Z (1 , 1/3, 1/3 ) , 

7.(0, 0,2/3), 

Z (0, 1/3, 2/3) , 

Z (1/3, 0,2/3) , 

Z (1/3, 1/3, 2/3) , 

Z (2/3, 0,2/3), 

Z (2/3, 1/3, 2/3) , 

Z (1, 0,2/3), 

Z (1,1/3, 2/3) , 

z (0,0,1), 

Z(0, 1/3,1), 

Z (1/3, 0,1) 

z (1/3, 1/3,1), 

Z (2/3, 0,1), 

Z (2/3, 1/3,1), 

z (1 , 0,1), 

Z(l, 1/3,1), 

Z (0, 2/3,0), 

Z(0,1,0) 

Z (1/3, 2/3,0), 

Z (1/3, 1,0) 

Z (2/3, 2/3,0), 

Z (2/3,1 ,0) 

Z (1, 2/3,0) , 

7.(1, 1,0) 

Z (0,2/3, 1/3), 

Z (0,1, 1/3) 

Z (1/3, 2/3, 1/3), 

Z(l/3, 1,1/3) 

Z (2/3, 2/3, 1/3)  , 

Z(2/3, 1,1/3) 

Z (1,2/3, 1/3) , 

Z (1, 1,1/3) 

Z (0,2/3, 2/3), 

Z (0,3 ,2/3) 

Z (1/3, 2/3, 2/3) , 

Z (1/3, 1,2/3) 

Z (2/3, 2/3, 2/3), 

Z (2/3, 1 ,2/3) 

Z (1,2/3, 2/3), 

Z (1,1, 2/3) 

Z (0, 2/3,1), 

Z  (0,1,1) 

Z (1/3, 2/3,1), 

Z (1/3, 1,1) 

Z (2/3, 2/3,1), 

Z (2/3, 1,1) 

Z(l, 2/3,1), 

Z(l,l,l) 
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5.3  Case  Control  Cards 

There  are  three  basic  categories  of  case  control  cards:  executive 
data,  case  definition,  and  output  requests.  The  following  sections 
detail  each  of  the  options  within  these  categories.  In  the  description 
of  the  format  for  individual  directives,  items  enclosed  in  parentheses 
are  optional  inputs.  The  item  "list"  represents  a  free-form  list  of 
integers,  separated  by  commas,  and  can  include  the  operators  ALL,  NONE, 
THRU,  and  EXCEPT. 

5.3.1  Execution  Control 

Tiie  executive  data  portion  of  the  case  control  deck  is  used  to 
control  the  executive  functions  CHKPNT  and  RESTART,  the  ay.isymmetric 
AXY,  as  well  as  the  global  parameters  TITLE,  DRY,  and  TIME.  It  is 
possible  to  use  RESTART  and  CHKPNT  options  on  the  same  run,  and  in  some 
instances  a  CHKPNT  file  is  automatically  generated.  The  management  of 
these  files  is,  of  course,  machine  dependent,  and  the  job  control 
language  (JCL)  directives  are  described  in  a  later  chapter. 
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5.3.2  Case  Definition 

The  case  definition  section  of  the  case  control  deck  is  used  to 

i 

!  activate  certain  sets  of  data  from  the  bulk  data  on  each  subcase  and  to 

!  define  the  subcase  title.  The  SUBCASE  card  defines  the  beginning  of  the 

next  subcase.  Within  any  subcase  at  least  one  of  the  set  identifiers 
LOAD  or  SDC  must  be  defined.  If  any  errors  are  detected  in  the 
definition  or  execution  of  a  subcase,  that  subcase  will  be  skipped.  Up 
to  15  subcases  may  be  defined,  including  the  zero  or  implied  subcase. 


5.3.3  Output  Requests 

The  output  requests  section  of  the  case  control  input  selects  the 
particular  data  to  be  printed.  Entry  into  this  region  is  initiated 
through  the  OUTPUT  card,  which  is  mandatory  if  any  output  requests  are 
to  be  made.  Up  to  25  sets  can  be  defined  using  SET  cards.  The  items  in 
these  sets  are  not  limited  to  any  particular  type  of  data  so  that  they 
can  serve  a  number  of  purposes.  The  individual  output  requests  are  made 
through  the  use  of  the  mnemonics  listed  in  Table  5-5.  The  majority  of 
these  cards  are  of  the  form  GRID,  set  where  set  can  be  a  list  of  SET 
identifiers. 


147 


CASE  CONTROL  OPTIONS 


PAGE  5-9 


Input  Directive: 
Description: 

Format: 

Example  Syntax: 
Remarks : 


CASE  CONTROL  INPUT 


AMBIENT  Temperature 

Defines  a  temperature  at  which  there  are  no 

thermal  strains. 

AMBIENT,  T 

AMBIENT,  70.0 

1.  If  not  specified,  the  ambient  temperature  is 
assumed  to  be  zero. 

2.  The  ambient  temperature  is  used  to  compute  the 
temperature  difference,  AT  =  T  -  TA,  for 
thermal  loads,  stresses,  and  strains. 
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Input  Pi  recti ve : 
Description: 

Format : 

Example  Syntax: 

Remarks : 


CASE  CONTROL  IN  PITT 


AXY  Axisymmetric  model 

Defines  a  model  as  axisymmetric  and  sets  the 
size  of  CCX  elements. 

.AXY,  A8 

AXY,  -15.0 


1.  The  subtended  angle  in  a  CCX  axisymmetric 
model  may  also  be  input  with  a  PARAM  directive 
in  the  bulk  data. 

2.  The  magnitude  and  sign  of  AO  must  agree  with 
the  hyperpatch  construction  for  the  OCX 
element  shape. 
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Inpot  Directive: 
Description: 

Format : 

Example  Syntax : 
Remarks: 


BULK  DATA  INPUT 


CHKPNT  Checkpoint  request 

Defines  the  restart  situations  for  which  a 

checkpoint  volume  will  be  written. 

CHKPNT,  checkpoint  1,  checkpoint  2 

CHKPNT,  ELEMENT,  03 

1.  Checkpoint  and  restart  may  not  be  requested 
on  the  same  run.  The  checkpoint  request  is 
ignored  in  this  situation. 

2.  A  file  or  tape  labeled  INPT  (CDC)  or  25. 
(UNIVAC)  must  be  requested  in  the  job  control 
deck. 

3.  CHKPNT,  ELEMENT  causes  all  element  stiffness 
matrices  to  be  saved  on  tape. 

4.  CHKFNT,  OG  causes  the  last  iterate  and 
direction  vector  to  be  saved  for  continued 
iteration  on  a  subsequent  restart  run. 

5.  Whenever  a  CG  solution  is  terminated  prior  to 
convergence  (due  to  a  TIME  insufficiency, 
maximum  iterations,  or  maximum  steepest 
descent  moves)  a  CG  checkpoint  is  appended  to 
the  restart  file.  Hius  the  file  or  tape  must 
have  write  permission. 
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BULK  DATA  INPUT 

Input  Directive: 

DRY  Dry  run  request 

Description: 

Automatically  halts  execution  after 
processing  of  all  modeling  data. 

Format : 

CRY 

Example  Syntax: 

DRY 

Remarks: 

1.  This  feature  allows  all  geometry, 

material,  load  and  boundary  condition  data 
to  be  checked  before  execution. 
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Input  Pi rective: 
Description: 

Format : 

Example  Syntax : 
Remarks: 


CASE  CONTROL  INPUT 


LOAD  Defines  active  load  set 

Defines  for  each  subcase  the  load  set  to  be 
used  for  analysis. 

LOAD,  ID 

LOAD,  10 

1.  "frie  LOAD  directive  selects  the  active 

load  set  from  all  those  present  in  the  bulk 
data. 

2.  Thermal,  as  well  as  mechanical,  load  sets 
are  controlled  with  this  directive. 
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Input  Directive: 
Description: 

Format : 

Remarks: 


BULK  DATA  INPUT 


OUTPUT  Output  selection 

Defines  the  beginning  of  the  output  requests 

reg  ion . 

OUTPUT 

1.  This  card  is  required  if  there  are  any  output 
requests . 

2.  All  output  flags  in  Table  5-5  can  be  selec¬ 
tively  turned  off.  The  modifier  OFF  is  used, 
for  example,  OLQAD  =  OFF,  to  prevent  output 
from  a  flag  activated  by  an  EVERYTHING 
request. 
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Input  Directive: 
Description: 

Format : 

Example  Syntax: 

Remarks: 


; 


CASE  CONTROL  INPUT 


RESTART  Restart  from  a  checkpoint  tape 

Specifies  the  checkpoint  volute  on  the  restart 
tape  to  be  used  on  the  restart  run.  Also  identi¬ 
fies  elements  to  be  modified. 

RESTART,  type,  (element  list) 

RESTART,  CG,  n  file,  STP 
RESTART,  ELEMENT 
RESTART,  ELEMENT,  6  thru  14,  18 
RESTART,  CG,  2,  S 

1.  The  element  list  identifies  those  elements 
that  are  to  be  modified  or  deleted. 

2.  All  case  control  and  bulk  data  cards  are  sub¬ 
mitted  on  a  restart. 

3.  In  a  CG  restart  request,  the  parameters  n  file 
and  STP  are  optional;  n  file  is  the  OG  restart 
file  nunber.  This  can  be  greater  than  1  for 
multiple  restarts  and  checkpoints. 

Default  =  1.  STP  is  request  for  the  restart 
to  begin  with  a  steepest  descent  move.  Any 
character  in  this  field  will  activate  the 
request.  The  default  is  no  steepest  descent 
move  in  vdiich  case  iteration  continues  from 
the  requested  restart  file  direction. 


4.  A  CG  restart  request  automatically  defines  an 
element  restart. 
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Input  Directive: 
Description: 

Format: 

Example  Syntax: 

Field 

n 

Remarks: 


CASE  CONTROL  INPUT 


SDC  Surface  displacement  constraint  selection 

Selects  the  surface  displacement  constraint  set 
to  be  applied  to  the  structural  model. 

SDC,  n 

SDC,  12 


Contents 


Set  identification  number  found  on  at  least  one 
surface  or  grid  point  displacement  card  within  the 
bulk  data  deck  (integer  >  0). 


1.  The  SDC  card  is  supplied  at  the  subcase  level. 

2.  The  total  load  applied  will  be  the  sum  of 
external  LOAD  and  SDC  loads. 

3.  At  least  one  LOAD  or  SDC  card  must  be  input 
for  each  subcase. 

4.  Single  point  constraint  sets  are  also  activa¬ 
ted  by  the  SDC  card. 

5.  Constraints  adequate  to  prevent  rigid  body 
motion  must  be  present. 
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Input  Directive: 
Description: 

Format: 

Example  Syntax: 

Field 

ID 

List 


CASE  CONTROL  INPUT 


SET  Set  definition 

Specifies  the  integer  identification  nimbers  of 
all  elements  in  the  set  being  defined. 

SET  ID,  List 

SET  7,  1,  2,  7,  10  THRU  25  EXCEPT  15  THRU  20 


Contents 


Identification  nisnber  of  the  set  being  defined. 

A  string  of  integers  and/or  the  list  operators 
ALL,  EXCEPT,  NONE,  THRU  that  define  the  elements 
in  the  set.  The  example  list  defines  SET  7  as 
containing  the  elements  1,  2,  7,  10,  11,  12,  13, 
14,  21,  22,  23,  24,  25.  All  specifications 
following  the  operator  EXCEPT  turn  off  items  that 
may  have  been  activated  earlier  on  this  card. 
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Input  Directive: 
Description: 
Format : 

Example  Syntax : 

Field 

n 


CASE  CONTROL  INPUT 

SUBCASE  Subcase  delimiter 
Delimits  and  identifies  a  subcase. 

SUBCASE,  n 
SUBCASE,  501 

Contents 

Subcase  identification  number  (integer  >  0). 
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CASE  CONTROL  INPUT 


SUBCQM  Subcase  combination 

Defines  a  linear  combination  of  subcases  for 

output  data  recovery. 

SUBCOM,  R^ ,  Rp,..,  Rn,/,  S] , . . ,  SM 

SUBCOM,  1.1,  -3.2,/,  1.5,,,  -6.7 

SUBCOM,  -1.047,  0.163 

1.  R^  is  the  load  factor  for  subcase  number 
i  generated  by  current  execution. 

2.  Sj  is  the  load  factor  for  subcase  number  i  on 
the  INDATA  file,  if  any. 

3.  There  are  no  default  values. 

4.  A  subcase  may  be  skipped  using  a  double  comma. 

5.  In  general,  to  recover  data  in  all  octants  of 
a  synmetry-syirmetry-symmetry  model  will 
require  selectively  changing  the  sign  of 
individual  displacement  components  before 
combination. 
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CASE  CONTROL  INPUT 


Input  Directive:  SUBTITLE  Output  subtitle 

Description:  Defines  a  subtitle  which  will  appear  on  the 

second  heading  line  of  each  page  of  PATCHES-III 
printer  output. 

Format :  SUBTITLE,  any  data  string 

Example  Syntax:  SUBTITLE,  PATCHES-III  SUBTITLE  FOR  SUBCASE  501. 

Remarks:  1.  SUBTITLE  will  title  output  for  the  subcase  in 

which  it  is  defined. 

2.  If  no  SUBTITLE  card  is  supplied,  the  SUBTITLE 
field  will  be  blank. 
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Input  Directive: 
Description: 
Format : 

Example  Syntax : 

Field 

Minutes 

Remarks : 


CASE  CONTROL  INPUT 

TIME  CP  time  estimate  in  minutes 

Defines  the  maximum  allowable  CP  time  in  minutes. 

TIME,  minutes 

TIME,  12 

Contents 

Integer  maximum  CP  time  in  minutes. 

1.  The  default  value  is  1  minute. 

2.  This  time  estimate  is  used  by  several  PATCHES 
utilities  to  check  the  time-to-go  before 
attempting  to  execute  a  module. 
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CASE  CONTROL  INPUT 


TITLE  Output  title 

Defines  a  title  which  will  appear  on  the  first 
heading  line  of  each  page  of  PATCHES-III  printer 
output . 

TITLE,  any  bed  data 

TITLE,  SAMPLE  TITLE  FOR  PATCHES-III 

Remarks :  1.  If  no  TITLE  card  is  supplied,  blanks  are 

assumed . 

2.  TITLE  information  is  also  written  onto  the 
restart  volume. 


CASE  CONTROL 


REQUEST 

ALL 

EVERYTH TNG 

GRID,  sot. 
DAT  AG ,  set 
LINEB,  set 
LINER,  set 
DLINEB,  set 
DLINEP,  set 
PATCHB,  set 
PATCHP ,  set 
PATCHA,  set 
DPATB ,  Set 
DP ATP,  set 
HPB,  set 
HPP,  set 
DHPB,  set 
DHPP,  set 
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Table  5-5 

CASE  CONTROL  OUTPUT  REQUESTS  AND  FLAGS 


EFFECT 


All  output  requests  are  turned  on  except  those 
subsequently  redefined  and  except  the  algebraic 
and  geometric  format  output.  Principal  frame  output 
and  several  diagnostic  data  are  also  turned  off. 

All  output  requests  are  turned  on  except  those 
which  are  redefined  subsequently. 

Grid  points. 

Data  grid  points. 

Geometric  lines  in  geometric  format. 

Geometric  lines  in  point  format. 

Data  lines  in  geometric  format. 

Data  lines  in  point  format. 

Geometric  patches  in  geometric  format. 

Geometric  patches  in  point  format. 

Geometric  patches  in  algebraic  format. 

Data  patches  in  geometric  format. 

Data  patches  in  point  format. 

Geometric  hyperpatches  in  geometric  format. 

Geometric  hyperpatches  in  point  format. 

Data  hyperpatches  in  geometric  format. 

Data  hyperpatches  in  point  format. 

(Table  continued  on  following  page.) 
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REQUEST 
DETJ,  set 

EDISP,  set 
ELEMENT,  set 

ESTRAIN,  set 
MSTRAIN,  set 
PSTRAIN,  set 
ESTRESS,  set 
MSTRESS,  set 
PSTRESS ,  set 
MATC,  set 
MATA,  set 
VOLUME,  set 
ELOAD,  set 

DTLINEB,  set 
DTLINEP,  set 
EFORCE ,  set 


Table  5-5  (Continued) 

CASE  CONTROL  OUTPUT  REQUESTS  AND  FLAGS 

EFFECT 

Determinate  of  geometric  hyperpatch's  Jacobians  at 
Gaussian  points. 

Element  Cartesian  displacema.ts  in  point  format. 

Element  results.  Simultaneously  activates  EDISP, 
ESTRAIN,  ESTRESS,  MSTRAIN,  MSTRESS,  PSTRAIN,  PSTRESS. 
ELEMENT,  ALL  is  assumed  unless  overridden. 

Element  Cartesian  strains  in  point  format. 

Material  frame  strains  in  point  format. 

Principal  frame  strains  in  point  format. 

Element  Cartesian  stresses  in  point  format. 

Material  frame  stress  in  point  format. 

Principal  frame  stresses  in  point  format. 

Element  elastic  constants  in  point  format. 

Element  thermal  expansion  constants  in  point  format. 

Volume  of  hyperpatches. 

Element  load  vectors  in  point  format  prior  to  trans¬ 
formation  to  analysis  coordinates. 

Data  table  lines  in  geometric  format. 

Data  table  lines  in  point  format. 

Element  node  forces  in  point  format. 

(Table  continued  on  following  page.) 
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Table  5-5  (Continued) 

CASE  CONTROL  OUTFUT  REQUESTS  AND  FLAGS 


FLAG 

EFFECT 

OLQAD 

Output 

the  global  load  vector. 

ODISP 

Output 

the  global  displacement  vector. 

EMESH 

Output 

from  the  matrix  assembly  module. 

LORDER 

Output  the  order  in  which  geometric  lines  were 
processed . 
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EXECUTIVE  CONTROL  OPTIONS 


6.1  Overview 


The  executive  control  deck  is  completely  optional  and,  in  fact,  is 
rarely  used.  The  executive  control  deck  provides  specialized  system 
level  control  over  the  execution  of  PATCHES-III.  This  section  should  be 
skipped  by  most  users  and  used  by  experienced  users  in  special 
circumstances . 


C  2  Executive  Control  Options 

Executive  control  cards  are  the  first  input  cards  in  an  execution 
and  always  start  with  PATCHES  as  the  first  nonblank  field  followed  by 
the  option  selections.  A  sample  executive  control  card  is  shown  below. 

PATCHES,  NODROP,  MAXFL  =  90112,  BREAKPT 
A  card  which  does  not  begin  with  the  letters  PATCHES  is  assumed  to  be 
the  first  card  of  the  case  control  deck.  The  current  executive  control 
options  are  described  below. 
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BRKAK1T 

.specifier,  that  this  run  is  to  create  random  access  breakpoint  files  at 
the  completion  of  each  execution  link  for  a  possible  restart.  7he  files 
are  stored  on  the  RASTUS  random  access  file  RNDM16.  This  option  should 
not  !x?  jon f used  with  the  case  control  breakpoint/restart  option.  An 
executive  control  restart  is  a  continuation  of  execution  from  the 
beginning  of  a  Specified  link  using  the  random  access  file  as  the  source 
ol  all  input  and  r  silts.  The  BREAKPT  request  adds  an  overhead  cost  of 
approximately  1  . 4  ;  -.coiicls  per  link  ili  to  60  seconds  for  a  typical 
execution' . 

MAXFL  =•  run 

specifies  a  maximum  field  length  above  which  the  program  will  not 
attempt  to  expand.  The  default  MAXFL  is  131000, or  377670R.  An 
attempt  by  t.jo  [.ron  s  to  exceed  this  limit  to  satisfy  an  open  core 
requirement  will  resu’t  in  a  fatal  diagnostic  and  termination  of 
execution.  The  form  is  MAXFL  -  nnn  where  rein  is  a  positive,  base  10 
integer . 

NURFt. 

does  not  permit  the  fiela  length  to  vary  during  execution.  Hie  field 
length  upon  entry  to  the  program  will  remain  in  effect  throughout  the 
run. 

RASTUS IN 

uses  the  existing  random  access  RASTUS  file  located  on  file  RNDM16  as 
the  basis  for  random  access  operations.  Otherwise  a  new  file  is 
generated. 
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EXECUTIVE  CONTROL  OPTIC' IS 


RESTART  =  n 

requests  that  execution  be  restarted  from  the  random  access  RASTUS  file 
RNDM16  from  the  beginning  of  link  n.  The  input  case  control  and  bulk 
data  decks  must  not  be  input,  and  only  the  EM)  DATA  card  should  be 
included.  All  options  and  data  will  be  restored  to  their  conditions  at 
the  time  of  the  BREAKPT  generation.  (With  one  exception:  If  a  case 
control  OG  restart  was  in  effect  and  completed  the  link  16  BREAKPT  then 
a  subsequent  RESTART  =  16  will  proceed  from  the  CG  restart  file 
generated  on  the  most  recent  execution.  In  other  words,  multiple  runs 
using  the  RESTART  =  16  and  BREAKPT  executive  control  cards  following  an 
original  RESTART  =  CG  case  control  request  will  result  in  continuous 
iteration  without  going  through  the  geometry  links,  etc.).  A  less 
complicated  example  would  be  if  a  run  with  BREAKPT  started  execution  of 
link  17  (displacement,  strain,  stress,  force  recovery)  and  terminated 
for  MAX  pages,  then  a  RESTART  =  17  would  restart  execution  from  the 
beginning  of  link  17. 

SCAFFOLD 

(Plot  system  or  generalized  postprocessor  only.)  Does  not  attempt  to 
read  the  PPDATA  file. 
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CHAPTER  7 

EXAMPLE  INPUT/OUTPUT 


7.1  Disk  Thermal  Stresses 


A  solid  disk  with  free  surfaces  is  heated  to  a  temperature 
distribution,  T  =  100  -  lOOOOr2  where  r  is  the  radius  in  inches.  This 
is  an  axisymmetric  problem  whose  solution  is  known  in  closed  form.  The 
PATCHED- II I  model  of  the  disk.  Figure  7-1,  illustrates  a  number  of 
important  modeling  features.  First,  element  number  one  has  a 
degenerated  surface  at  the  center  of  the  disk.  The  bulk  data  for  the 
model  demonstrate  that  no  special  input  are  required  to  either  create 
this  hyperpatch  or  identify  the  element  as  degenerate. 

Note  that  the  parameterization  of  both  hyperpatches  is  determined 
from  the  connectivity  data  (CPDE3  cards)  per  Figure  4-1  and  not  from  the 
HPR  card  parameterization.  Also,  since  the  through  the  thickness 
behavior  of  this  problem  is  well-behaved,  half  of  the  degrees  of  freedom 
of  the  model  can  be  eliminated  at  little  cost  in  accuracy  by 
constraining  the  element  to  be  linear  in  the  ^  direction  (the  OCL 
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parameter  on  the  CPDE3  card).  Second,  these  data  illustrate  the  input 
of  symmetry  boundary  conditions.  Note  that  the  centerline  is  fully 
restrained  in  the  (1,2)  plane  as  a  result  of  two  surface  normal 
constraints  (SDC20  cards)  on  surfaces  that  intersect  along  the 
centerline.  To  constrain  rigid  body  motion  in  the  direction,  grid 
point  one  has  component  number  three  constrained  via  an  SPC1  card.  The 
synthesis  of  constraints  at  this  grid  point  then  involves  both  surface 
and  grid  point  data  cards.  Third,  several  data  modeling  features  are 
illustrated  by  the  temperature  input.  The  algebraic  format  data  patch 
is  particularly  simple  in  this  case;  only  two  of  the  sixteen  parameters 
are  nonzero  for  data  patch  10  and  three  for  data  patch  20.  (See  lines  3 
to  10  of  the  bulk  data.)  Note  that  r  =  0.05^  for  element  one  and 
r  =  0.05 (^  +  1)  for  element  two.  Since  the  temperature  field  is 
axisyimetric,  the  data  patch  equivalence  option  can  be  used  to  create 
data  patches  30  and  40.  The  data  hyperpatches  are  created  using  DHP2P 
cards  where  the  parameterization  is  consistent  with  the  geometric 
hyperpatches.  Were  they  not  consistent,  the  DHPSORT  card  would  be  used 
to  reparameterize  the  data  hyperpatches.  The  stress  results  from 
PATCHES-III  are  compared  to  the  closed-form  solution  (Reference  3)  in 
Table  7-1.  The  maximum  difference  is  les=-  tan  one  percent  (1%)  with 
respect  to  the  peak  stress.  Note  that  '  .  r  the  outer  radius  is  a 
natural  boundary  condition  which  is  very  closely  approximated  with  only 
a  two-element  model.  Finally,  Table  7-1  shows  that  subtended  angle  has 
little  effect  on  modeling  error  in  this  problem.  The  90-degree 
subtended  angle  model  has  virtually  the  same  stress  accuracy  as  the 
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30-degree  model .  This  is  somevdiat  surprising  because  a  circular  arc 
cannot  be  modeled  exactly  with  a  parametric  cubic,  and  the  geometry 
modeling  error  for  the  90-degree  case  is  much  higher  than  the  30-degree 
case. 
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Table  7-1 

DISK  THERMAL  STRESS  COMPARISONS 


PATCHES-III  (PSI)  EXACT  SOLUTION  (PSI) 


= 

30° 

$  = 

90° 

lOr 

0 

a 

r 

-24.97 

°<f> 

-24.97 

a 

r 

-25.14 

°4> 

-25.13 

a 

r 

-25.00 

°<J> 

-25.00 

1/6 

-23.93 

-22.63 

-24.26 

-22.97 

-24.31 

-22.92 

1/3 

-21.81 

-16.29 

-22.08 

-16.55 

-22.22 

-16.67 

1/2 

-18.59 

-  5.99 

-19.17 

-  6.40 

-18.75 

-  6.25 

2/3 

-13.57 

8.64 

-13.72 

8.59 

-13.89 

8.33 

5/6 

-  7.32 

27.39 

-  7.91 

27.02 

-  7.64 

27.08 

1 

.24 

50.27 

.10 

50.58 

0.00 

50.00 
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f 

f 

. 


E  =  107  PSI 
a=  10  7/F° 


0100608 


Figure  7-1.  Disk  Thermal  Stress  Model 
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PATCHES  CONTROL  DIRECTIVES 

TITLE,  DISK  THERMAL  STRESS  PROBLEM 

TIME, 5 

RESTART,  CG 

SDC,  10 

LOAD,  10 

OUTPUT 

EVERYTHING 

BEGIN  BULK 
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EXAMPLE  INPUT /OUT PUT 
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CPDE3, 

PATCHES  DATA  DIRECTIVES 

1,  1,  1,  7,  2,,,,  CCL,  4,  4,  9, 

5 

CPDE3, 

- ,  2,  7,  8,  3 ,,,,  CCL,  5,  9,  10, 

f> 

DPATA , 

10,  1,  1,  7,  2,,,,  13(0),  -25.0,  0, 

100.0 

DPATA , 

20,  2,  7,  8,  3,,,,  13(0),  -25.0,  -50. 

0,  75.0 

UPATf'U, 

30,  10,  4,  4,  9,  5 

DPATKC, 

40,  20,  5,  9,  10,  6 

DHP2P, 

1,  10,  30 

DHP2P , 

2,  20,  40 

TEMP, 

10,  1,  1 

TEMP, 

10,  2,  2 

GRID, 

1,,  0.0,  0.,  0. 

GRID, 

2,,  0.05,  0.,  0. 

GRID, 

3,,  0.10,  0.,  0. 

GRID, 

4,,  0.0,  0.,  0.01 

GRID, 

5,,  0.05,  0.,  C.01 

GRID, 

6,,  0.10,  0.,  0.01 

PATCHQ, 

11,  1,  2,  5,  4 

PATCHQ , 

12,  2,  3,  6,  5 

HPR, 

1,  11,,,,  0.,  90.,  3 

HPR, 

2,  ±2, r , ,  0.,  90.,  3 

MATAL, 

1,  1,  1,  2 

MAT!  , 

1,  1 .  +7 , ,  0.3 

MTRX-2, 

1.0-7,  1.0-7,  1.0-7 

PPDE3, 

1,  1 

PPDE3, 

2,  1 

SDC20, 

10,  2,  3,  7,  8,  10,  9 

SDC20 , 

10,  2,  3,  2,  3,  6,  5 

SDC20 , 

10,  1,  3,  1,  7,  9,  4 

SDC20, 

10,  1,  3,  1,  2,  5,  4 

SPC1, 

10,  1,,,  0.0 

END  DATA 
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7.2  Interlaminar  Normal  Stresses 

One  of  the  few  three-dimensional  composite  laminate  problems  for 
which  corroborative  solutions  exist  is  a  four-ply  graphite-epoxy  plate 
under  uniaxial  load.  A  finite  difference  solution  of  the  elasticity 
equations  (Reference  4),  a  stress-function  discrete  element  solution  of 
the  elasticity  equations  (Reference  5),  an  analytic  solution  of  certain 
higher  order  laminated  plate  theory  equations  (Reference  6) ,  and  the 
PATCHES -I I I  displacement-function  discrete  element  solution  of  the 
elasticity  equations  all  agree  well  for  the  interlaminar  normal  stress. 
■Stresses  are  computed  for  the  0u/90o/90n/0°  laminate  shown  in  Figure  7-2 
under  a  uniform  imposed  displacement  in  the  axial  direction.  Taking 
full  advantage  of  symmetry,  a  simple  four-element  PATCHES-III  model  was 
created  from  46  bulk  data  directives.  There  are  three  planes  of 
symmetry  all  modeled  with  SDC10  cards.  The  imposed  axial  displacement 
which  constrains  only  the  ej  component  on  the  2(1,  £3)  surfaces  is 
modeled  with  SDC1  directives.  Note  that  blank  fields  for  the  e2,  e?, 
data  patches  on  the  SDC1  directives  allow  these  components  to  remain 
unconstrained.  The  material  property  modeling  required  only  one 
materials  matrix  since  all  lamina  are  the  same  except  for  orientation. 
The  Euler  angle  option  (PSI  =  90* ,  THETA  =  PHI  =  0°)  was  used  on  the 
PPDE3  property  directive  for  elements  1  and  2  and  to  obtain  the 
materials  matrix  for  a  90°  lamina.  Using  this  feature  it  would  be  a 
trivial  data  change  to  model  a  45°/-450/-45°/450  laminate.  The  nature 
of  the  solution  made  it  possible  to  use  two  tri-linear  elements  leading 
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into  the  linear-cubic— cubic  elements  at  the  free  edge. 

A  view  of  the  deformation  of  the  center  cross-section,  Figure  7-3, 
indicates  the  local  nature  of  the  distortion  caused  by  the  free  edge. 
The  interlaminar  normal  stress,  Figure  7-4,  peaks  at  the  free  edge  and 
damps  out  quickly  as  expected.  Comparison  with  the  reference  solutions 
in  this  figure  show  excellent  agreement.  The  printed  output  shows  that 
the  displacement  at  any  point  in  the  cross-section  is  independent  of  7.1 , 
as  assumed  in  Reference  4,  and  that  the  axial  strain,  e,j  ,  is  constant. 
Table  7-2  shows  the  interlaminar  normal  stress  for  three  models  that  use 
the  constant  strain  information  to  reduce  the  cost  of  the  analysis. 


Table  7-2 

INTERLAMINAR  NORMAL  STRESS  COMPARISONS* 


f/2h** 

CCC/CCC 

(428  D.O.F. ) 

LLL/LCC 
(92  D.O.F.) 

LLL/LCL 
(44  D.O.F. 

0 

+2.95 

+2.89 

+2.76 

1/3 

-  .26 

-  .27 

-  .31 

2/3 

-  .43 

-  .44 

-  .25 

1 

-  .16 

-  .05 

-  .03 

2 

-  .02 

-  .03 

-  .03 

3 

+  .01 

-  .01 

-  .00 

4 

-  .01 

+  .02 

+  .02 

♦Comparisons  at  midsurface  between  90-degree  plies. 
**Z7  =  b  -  r  ;  distance  from  free-edge. 
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Figure  7-2.  Model  of  a  Graphite  Epoxy  Laminate 


Figure  7-4. 


Interlaminar  Normal  Stress  Comparison 
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PATCHES  CONTROL  DIRECTIVES 

TITLE,  INTERLAMINAR  STRESS  PROBLEM 

TIME,  8 

SDC,  10 

OUTPUT 

EVKRYrHING 

BEGIN  BULK 
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PATCHES  DATA  DIRECTIVES 


CPDE3, 

1, 

14, 

8,  7, 

13,,,,  LLL 

CPDE3, 

2, 

S, 

2,  1, 

7, 

,,,  LCC 

CFDE3, 

3, 

16, 

10,  9, 

15,,,,  LLL 

CPDE3, 

4, 

10, 

4,  3, 

9,,,,  LCC 

DATAG, 

1» 

,  1, 

0.01,  3 

,  o. 

01,  5,  i 

0. 

01, 

15,  0.01 

,  17 

,  0.01 

DPATQ, 

10, 

If 

13,  15, 

9, 

7 

DPATQ, 

20, 

1, 

7,  9, 

3, 

1 

DPATQ, 

30, 

1, 

15,  17, 

11. 

9 

DPATQ, 

40, 

If 

9,  11, 

5, 

3 

GRID, 

15, 

f  8. 

,  0.,  0. 

GRID, 

9, 

,  8. 

,  6. ,  0. 

GRID, 

3, 

f  8. 

,  8.,  0. 

GRID, 

16, 

,  o. 

,  0.,  0. 

GRID, 

10, 

,  o. 

,  6. ,  0 

GRID, 

4, 

,  o. 

,  8.,  0 

HFN, 

1, 

If 

-.5 

HFN, 

2, 

2, 

-.5 

HPN, 

3, 

If 

.5 

HPN , 

4, 

2, 

.5 

PPDE3, 

1. 

1, 

,  90. 

PPDE3, 

2, 

If 

,  90. 

PPDE3, 

3, 

1 

PPDE3, 

4, 

1 

MATC, 

1, 

If 

1.  1 

MTRX-1, 

20. 

2+6, 

0.56+-6, 

0.56+6,  2.; 

0. 

0,  0 

.85+6,  0 

.0, 

0.85+6 

PATCHQ , 

If 

16, 

10,  9, 

15 

PATCflQ, 

2, 

10, 

4,  3, 

9 

GDC  10, 

10, 

1, 

1,  14, 

16, 

10,  8 

SDC10, 

10, 

1, 

2,  14, 

13, 

15,  16 

GDC10, 

10, 

If 

3,  14, 

8, 

7,  13 

SDC10, 

10, 

2, 

1,  8, 

10, 

4,  2 

SDC10, 

10, 

2, 

3,  8, 

2, 

1,  7 

SDC10 , 

10, 

3, 

1,  16, 

18, 

12,  10 

SDC10, 

10, 

3, 

2,  16, 

15, 

17,  18 

SDC10, 

10, 

4, 

1,  io, 

12, 

6,  4 

SDC1, 

10, 

1, 

10 

SDC1, 

10, 

2, 

20 

SDC1, 

10, 

3, 

30 

SDC1 , 

10, 

4, 

40 

END  DATA 
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7 . 3  Thick  Cylinder  Pressure  Stresses 

A  widely  used  three-dimensional  test  case  for  mechanical  loads  is  a 
thick-walled  cylinder  under  pressure.  An  exact  solution  for  this 
problem,  the  so-called  Lame  cylinder,  is  available  in  Reference  7.  A 
ono-element  PATCHES-III  model  of  the  cylinder,  Figure  7-5,  was  loaded  by 
an  internal  pressure  of  5  psi  and  an  external  pressure  of  10  psi.  The 
principal  new  modeling  feature  demonstrated  by  this  problem  is  the 
pressure  load  card,  PLGAD3.  Note  that  a  negative  scale  factor  was  used 
for  data  patch  30  because  the  positive  normal  on  this  surface  is 
directed  away  from  the  outer  surface  of  the  cylinder. 

The  stress  results  for  three  different  subtended  angles  are 
compared  in  Table  7-3  with  the  exact  solution.  Again,  as  with  the  disk 
problem,  there  is  little  loss  in  accuracy  even  in  the  90-degree  model. 
The  maximum  error  is  approximately  three  percent  (3%)  relative  to  the 
peak  stress  using  a  single  CCL  element  to  model  one  quadrant. 
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Table  7-3 

CYLINDER  PRESSURE  STRESS  COMPARISONS 


r 

°4> 

(psi) 

0 

r 

(psi) 

EXACT 

A<{>=10 

A<f>=45 

Aij>=90 

EXACT 

A<fr=10 

A4>=4  5 

A<(>=90 

i 

-18333. 

-18464. 

-18450. 

-18516. 

-5000. 

-5512. 

-5503. 

-5554, 

4/3 

-15417. 

-15340. 

-15319. 

-15276. 

-7917. 

-7757. 

-7741. 

-7763, 

5/3 

-14067. 

-14114. 

-14091. 

-14014. 

-9267. 

-9392. 

-9372. 

-9388 

2 

-13333. 

-13172. 

-13143. 

-12995. 

-10000. 

-9726. 

-9703. 

-9676. 

EXAMPLE  IN PUT /OUT PUT 
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PATCHES  CONTROL  DIRECTIVES 


TITLE,  LAME  CYLINDER 
TIME,  6 

CHKPNT,  ELEMENT,  CG 

SDC,  10 

LOAD,  10 

OOTPUT 

EVERYTHING 

BEGIN  BULK 
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PATCHES  DATA  DIRECTIVES 


CPDE3, 

1, 

1, 

5,  8,  4,,,, 

CCL, 

2, 

GRID, 

1, 

2, 

1.0,  0.0,  1.0 

GRID, 

2, 

2, 

1.0,  0.0,  2.0 

GRID, 

3, 

2, 

2.0,  0.0,  2.0 

GRID, 

4, 

2, 

2.0,  0.0,  1.0 

DATAG, 

1,, 

1, 

1000.,  2,  1000. 

,  6, 

1000 

7, 

1000 

.,  8,  1000. 

DPATQ, 

20, 

1. 

1,  2,  6,  5 

DPATQ, 

30, 

1, 

4,  3,  7,  8 

HPR, 

1, 

1». 

,,,  10.0,  3 

MATC, 

1, 

1, 

1,  1 

MTRX-1 , 

.40385+7,  .17308+7,  .17308+7,  . 

1.1538+6 

,  0.,  0.,  1.1538+6, 

0., 

PLQAD3, 

10, 

1, 

20,  5.0 

PLOAD3, 

10, 

1, 

30,  -10.0 

PATCHQ, 

1, 

1, 

4,  3,  2 

PPDE3, 

1. 

1 

SDC20, 

10, 

1, 

3,  5,  8,  7, 

6 

SDC20, 

10, 

1, 

3,  1,  4,  3, 

2 

SDC10, 

10, 

1, 

3,  2,  6,  7, 

3 

SDC10, 

10, 

1, 

3,  1,  5,  8, 

4 

END  DATA 

7,  3 
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7.4  .Symmetry-Asymmetry  Analysis  of  a  Filamentary  Composite 

The  efficient  modeling  of  structures  often  requires  using  any 
symmetries  or  asymmetries  present.  In  general,  this  requires  the 
superposition  of  several  cases  to  obtain  the  complete  solution.  The 
SUBCOM  case  control  card  automates  this  procedure,  and  its  use  is 
illustrated  using  the  fiber-reinforced  structure  shown  in  Figure  7-6. 
The  reinforcing  fibers  are  oriented  at  -45  degrees  in  the  one-two  plane. 
The  Z(f  ,  0,  £  )  boundary  is  restrained  in  the  e  direction,  and  a 
uniform  tension  is  applied  to  the  opposite  face  Z (£  ,  1,  f  ) .  This  load 
can  be  factored  into  the  sum  of  a  symmetr ic-symmetric  load  and  an 
asymmetric-asymmetric  load  in  the  material  axes,  e'j .  A  structural 
symmetry  model.  Figure  7-6,  was  first  analyzed  with  symmetry-symmetry 
boundary  conditions  and  the  results  stored  on  a  postprocessing  data 
file.  Then  on  a  second  run,  the  asymmetry-asymmetry  solution  was 
obtained  and,  using  the  SUBCOM  feature,  added  to  the  previous  solution 
to  obtain  the  complete  solution.  The  control  cards  for  this  procedure, 
shown  in  Section  8.2  of  this  manual,  are  very  simple.  Hie  results  from 
the  combined  solution  agree  to  six  places  with  a  single-element  model  of 
the  complete  cubical  structure.  The  shear  stresses,  which  should  be 

_p. 

identically  zero,  are  on  the  order  of  1.0  X  10  psi  in  the  full  cube 
model  and  1 .0  x  10  psi  in  the  symmetry  model.  Note  that  although  the 
n  displacement  component  is  zero  for  the  asymmetry  case,  this  will  not 
txi  true  generally  and  is  not  part  of  the  asymmetry  constraints. 
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|  The  strain  results  for  the  three  cases  are  summarized  in  Table  7-4. 

i 

t 


Table  7-4 

SYMMETRY-ASYMMETRY  MCDEL 


Case 

Symmetry 

Asymmetry 

Combined 


Ell 

3.8024-7 
-7.1429-7 
-  .3340-6 


E22 

3.8024-7 

7.1429-7 

1.0945-6 


E33 

-2.3708-7 

0.0 

-2.3708-7 
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8100606 


Figure  7-6.  Symmetry  Modeling  of  a  Filamentary  Composite 


EXAMPLE  INPITT/OUTPUT 
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PATCHES  CONTROL  DIRECTIVES 

TITLE,  ANGLE  PLY  LAMINA,  SUBCOM  TEST  CASE 

SUBTITLE,  SYf-WETRY  BOUNDARY  CONDITIONS 

TIME,  2 

LOAD,  5 

SDC,  10 

OUTPUT 

ALL 

BEGIN  BULK 


i 

i 
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PATCHES  DATA  DIRECTIVES 

CPDE3,  1,  1,  2,  3,  l,,,,,  4,  5,  6,  4 

DATAG,  1,,  2,  1.0,  3,  1.0,  5,  1.0,  6,  1.0 

DPATQ,  10,  1,  2,  3,  6,  5 

GRID,  1,,  0.0,  0.0 

GRID,  2,,  -1.0,  1.0 

GRID,  3,,  1.0,  1.0 

PATC11Q,  1,  1,  2,  3,  1 

HPN,  1,  1,  2.0 

MATC,  1,  1,  1,  1 

MTRX-1  20.35+6,  .5888+6,  .5888+6,  1.501+6,  .4621+6,  1.501+6,  0.7+6, 
0.0,  0.0,  0.7+6,  0.0,  0.7+6 
PPDE3,  1,  1,,  -45.0 

PLOAD3,  5,  1,  10,  1.0 

SDC20,  20,  1,  1,  1,  4,  6,  3 

SDC20,  20,  1,  1,  1,  4,  5,  2 

SDC20,  10,  1,  3,  1,  4,  6,  3 

SDC20,  10,  1,  3,  1,  4,  5,  2 

END  DATA 
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PATCHES  CONTROL  DIRECTIVES 


TITLE,  ANGLE  PLY  LAMINA,  SUBCOM  TEST  CASE 

SUBTITLE,  ASYMMETRY  BOUNDARY  CONDITIONS 

TIME,  2 

LOAD,  5 

SDC,  20 

SUBCASE,  2 

SUBTITLE,  COMBINED  SOLUTION 
SUBCOM,  1.0,  /,  1.0 


OUTPUT 


BEGIN  BULK 
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PATCHES  DATA  DIRECTIVES 


CPDE3, 

1, 

1,  2,  3, 

1, 

,,,,  4,  5,  6,  4 

DATAG, 

1», 

2,  1.0,  3, 

1.0 

,  5,  1.0,  6,  1.0 

DPATQ, 

10, 

1,  2,  3, 

6, 

5 

GRID, 

1,, 

0.0,  0.0 

GRID, 

2,, 

-1.0,  1.0 

GRID, 

3,, 

1.0,  1.0 

PATCHQ, 

1, 

1,  2,  3, 

1 

HPN, 

1, 

1,  2.0 

MATC, 

1, 

1 ,  1,  1 

MTRX-1 , 

20.35+6,  .5888+6 

9  • 

5888+6,  1.501+6,  .4621+6, 

0.0, 

0.0,  0.7+6 

,  o 

.0,  0.7+6 

PPDE3, 

1, 

1,,  -45.0 

PL0AD3, 

5, 

1,  10,  1.0 

SDC20, 

20, 

1,  1,  1, 

4, 

6,  3 

SDC20, 

20, 

1,  1,  1, 

4, 

5,  2 

SDC20, 

10, 

1,  3,  1, 

4, 

6,  3 

SDC20, 

10, 

1,  3,  1, 

4, 

5,  2 

END  DATA 

1.501+6, 


0.7+6, 
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CHAPTER  8 


JOB  CONTROL  OPTIONS 


8.1  Overview 


Only  one  control  card  is  necessary  to  execute  the  basic  PATCHES-III 
system  on  a  CDC  computer.  At  other  sites  only  two  cards  are  required. 
Versions  9.4  and  higher  of  the  program  operate  under  the  NOSBE  operating 
system.  The  program  was  compiled  using  the  FTN  4.7  compiler.  The 
following  section  shows  the  job  control  directives  necessary  for  each  of 
the  basic  PATCHES-III  executions  using  NOSBE  directives.  When  executing 
under  NOS  the  CATALOG  directives  must  be  changed  to  DEFINE  directives, 
and  other  minor  variations  may  occur,  depending  on  how  physical  tapes 
are  requested.  Other  minor  variations  may  also  be  necessary  at 
individual  computer  facilities  to  account  for  local  system  differences. 


8.2  Example  Deck  Setups 

The  following  group  of  deck  setups  assumes  that  the  user  has 
already  input  the  appropriately  formatted  job  and  accounting  cards.  The 
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CP  time  estimate  should  be  at  least  as  large  as  the  case  control  TIME 
directive.  The  maximum  field  length  parameter  should  accommodate  the 
maximum  open  core  requirement  for  the  problem  (try  200000  initially). 
The  PATCHES  file  may  require  an  ID=XXXX  entry  on  the  ATTACH  card  at  some 
sites. 


1.  Execution  of  a  new  model 


ATTACH,  PATCHES,  MR=1. 

PATCHES,  PL=9999999 . 

EOR  Card 

Case  and  bulk  data  decks 
EOF  Card 

2.  Execution  of  a  new  model  and  generation  of  a  checkpoint  file 

• 

REQUEST,  INPT,  *PF. 

ATTACH,  PATCHES,  MR=1. 

PATCHES,  PL=9999999 . 

CATALOG,  INPT,  RP=999. 

EOR 

Case  and  bulk  data  decks 
EOF 

3.  Restart  from  a  checkpoint  tape 

• 

ATTACH,  INPT,  PW=EXTEND,  MODIFY. 

ATTACH,  PATCHES,  MR=1. 

PATCHES,  PL=9999999. 

EOR 

Case  control  with  RESTART  request  and  bulk  data 
EOF 
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4.  Execution  of  a  new  model,  saving  the  PPDATA  file 


• 

REQUEST,  PPDATA,  *PF. 

ATTACH,  PATCHES,  MR=1 
PATCHES,  PL=9999999 . 

CATALOG,  PPDATA,  RP=999. 

EOR 

Case  and  bulk  data  decks 
EOF 

5.  Execution  of  a  SUBCOM  utilizing  a  previous  PPDATA  file, 
outputting  a  new  PPDATA  file 

• 

REQUEST,  PPDATA,  *PF. 

ATTACH,  INDATA,  old  PPDATA  file. 

ATTACH,  PATCHES,  MR=1. 

PATCHES,  PL=9999999. 

CATALOG,  PPDATA,  RP*999. 

EOR 

Case  control  with  SUKX3M  and  bulk  data  deck. 

EOF 

6.  Save  the  random  access  file  for  later  use. 

• 

REQUEST,  FNDM16,  *PF. 

ATTACH,  PATCHES,  MR-1. 

PATCHES,  PL=9999999 . 

CATALOG,  FNDM16,  RP=999. 

EOR 

Executive  deck  perhaps  with  BREAKPT  option  case  control 
and  bulk  data  decks 
EOF 
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7.  Execute  a  new  model  and  create  a  plot  file 

• 

REQUEST,  TAPE21 ,  *PF. 

ATTACH,  PATCHES,  MR=1. 

ATTACH,  PATPLOT,  MR=1. 

PATCHES,  PL=9999999. 

PATPLOT. 

CATALOG,  TAPE21,  plot  file  name,  RP=999. 
EOR 

Case  and  bulk  data  decks 
EOR 

Plotting  directives 
EOF 
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8.3  Program  Files 

The  PATCHES  file  names  in  order  are  as  follows: 


Default  File  Name 

INPUT 

OUTPUT 

TAPE1 

TAPE  2 

TAPE  3 

TAPE  4 

TAPE  8 

TAPE  9 

TAPE 10 

INPT 

PPDATA 

INDATA 

RNDM16 

TAPE 17 

TAPE  18 

TAPE 20 


Contents 

Input  directives 

Output  print  file 

Alternate  input  file 

Bulk  data  sequence 

Case  control  cards 

Line  sequence  network 

Bulk  data  output,  scratch  file 

Scratch  file 

Bulk  data  cards 

Checkpoint/restart  file 

Postprocessing  data  file 

Input  of  a  previous  PPDATA  file 

RASTUS  file,  random  access 

SUBCOM  request  file 

Scratch  file 

User  information  file 
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8.4  PATPLOT  Postprocessor 

PATCHES-III  output  can  be  plotted  using  the  PATPLOT  postprocessor 
which  operates  on  the  FNDM16  file  created  during  a  CDC  execution.  This 
file  on  the  VAX-11/780  is  labeled  RASTUS.DAT.  The  PATPLOT  program  is 
normally  used  in  an  interactive  mode  and  prompts  the  user  for 
information  about  what  data  are  to  be  plotted  on  each  frame.  The  user 
responds  from  his  terminal  to  each  prompt  with  a  line  of  free-form 
input.  The  syntax  is  exactly  as  described  in  Chapter  4  for  free-form 
bulk  data  input.  The  system  can  also  be  used  in  batch  mode  by  putting 
all  the  responses  in  one  sequential  input  record  following  the  CDC  job 
control  directive  that  executes  PATPLOT. 
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1.  PATPLOT  interactive  mode  plot  file  creation 

• 

Plotting  from  run  entitled  ...  DISK  THERMAL  STRESS  PROBLEM. 
Input  list  of  elements  to  plot. 

3(1),  3(2) 

Input  corresponding  list  of  faces. 

2(3, 4, 5) 

Input  the  plot  subtitle. 

CARTESIAN  STRESS  Sll 

Input  displacement  magnification  factor. 

0 

Input  data  component:  0  =  NONE,  1-3  =  DISP,  4-9  =  STRAIN, 
10-15  =  STRESS,  16-21  =  MATERIAL  FRAME  STRAIN, 

22-27  =  MATERIAL  FRAME  STRESS. 

10 

Input  viewing  angles  about  the  X,Y,Z  axes,  respectively. 

45,  0,  -45 

Input  hidden  lines  flag:  0  =  ALL,  1  =  HIDDEN,  2  =  BOTH. 

0 

Input  number  of  U/W  lines/patch. 

1 

Input  contours  flag:  0  =  NO,  -1  =  DELTA,  N  =  N  CONTOURS. 

-1 

10 

Input  carpet  flag:  0  =  NO,  1  =  YES. 

0 

Input  smoothing  flag:  0  =  NO,  1  =  YES. 

1 

Working. 

Input  list  of  elements  to  plot. 

0 


The  above  PATPLOT  session  will  plot  one  frame  showing  Cartesian 
stress  contours  on  faces  3,  4,  and  5  of  elements  1  and  2.  The  contours 
will  be  in  10  psi  intervals,  as  illustrated  in  Figure  8-1.  The  session 
was  ended  by  inputting  a  zero  in  response  to  the  last  prompt  for  the 
number  of  elements  in  the  next  frame.  The  same  session  executed  in 
batch  mode  on  a  CDC  computer  would  have  the  following  input. 
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2.  PATPLOT  batch  mode  plot  file  creation 

• 

ATTACH  (RNDM16,  ID  =  ...) 

REQUEST  (TAPE21 ,  *PF) 

RFL  (100000) 

PATPLOT 

CATALOG  (TAPE 21 ,  PLOTID,  ID  =  ...) 

*EOR 

3(1),  3(2) 

2(3,4, 5) 

CARTESIAN  STRESS  Sll 
0 

10 

45,  0,  -45 
0 
1 

-1 

10 

0 

1 

0 

*EOJ 

In  both  the  batch  mode  and  interactive  mode  sessions,  a  plot  file 
of  CALCOMP  compatible  instructions  has  been  generated  for  processing  on 
a  hardcopy  plotter.  The  interface  with  a  specific  hardware  device, 
i.c.,  CALCOMP,  TRILOG,  etc.,  is  a  site  dependent  operation.  The  format 
of  tiie  plot  file  is  industry  standard,  which  makes  this  a  rou'-ine 
operation.  There  is  also  a  version  of  PATPLOT  available  for  output 
directly  on  a  TEKTRONIX  4014,  which  is  interactive  in  a  truly  graphic 
sense.  The  prompting  is  exactly  the  same,  and  the  same  RNDM16  file  is 
used  by  the  interactive  graphics  version  of  PATPLOT. 
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An  examination  of  the  sample  session  file  shows  that  each  element 
face  plotted  must  be  defined  individually.  To  plot  all  six  faces  of  an 
element  then  requires  the  element  number  to  be  input  six  times  in  the 
element  list.  This  feature  places  complete  control  in  the  user's  hands 
for  composing  each  frame  plotted.  Any  or  all  faces  of  any  or  all 
elements  may  be  examined  for  any  response  mode.  This  allows  regions  as 
small  as  a  single  element  to  be  excised  and  viewed  in  great  detail 
(Figure  8-2),  or  the  overall  deformations  of  the  entire  model  can  be 
viewed  (Figure  8-3)  at  any  magnification. 


201 


OtSK  THERWRL  STRESS  PRBSLEfl 
RRTES1RN  STRESS  Sll 


Figure  8-1.  Disk  Thermal  Stress  Contour  Plot 


Figure  8-2.  Interlaminar  Normal  Stress  Contour  Plot 


UNOEFBRMCD  STRUCTURE 


Figure  8-3.  Undefonned  Structure  Plot 


CHAPTER  9 


DIAGNOSTICS 


9.1  Overview 


The  executive  and  data  processing  systems  in  PATCHES-III  have  been 
designed  to  diagnose  as  many  errors  as  possible  on  any  execution.  All 
cards  are  processed  and  cross  referenced,  and  v*ienever  possible  the 
remaining  fields  of  a  card  in  vrtiich  errors  have  been  detected  are  also 
scanned . 

In  the  case  data  region,  the  error  message  is  listed  on  the  line 
following  the  erroneous  card.  Errors  detected  in  the  OUTPUT  region  of 
case  control  are  not  fatal. 

In  the  bulk  data  region,  the  card  vrfiich  precipitated  the  error  is 
listed  along  with  a  description  of  the  error  condition.  A  run  in  v#iich 
errors  were  detected  will  terminate  at  the  DRY  location  unless  the 
ERRORS  PARAM  field  is  overridden. 
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9.2  Error  Messages 

The  error  messages  in  PATCHES-III  are  generally  self-explanatory, 
and  the  number  of  possible  messages  is  so  large  as  to  make  even  an 
enumeration  of  possible  error  conditions  unreasonable.  Therefore,  this 
section  will  identify  and  explain  only  the  most  comnon  or  more 
complicated  errors.  Letters  in  small  type  within  an  all-caps  error 
message  text  represent  variables. 

THE  CORNER  GRID  POINT  VALUES  RESULTING  FROM  TOE  CONSTRUCTION  OF 
HYPERPATCH  n  CANNOT  BE  ORIENTED  TO  MATCH  THE  PREVIOUSLY  DEFINED  GRID 
POINT  VALUES. 


GRID  POINTS  INVOLVED 

9l 

<?2 

93 

% 

95 

% 

g7 

9b 

ORIGINAL  COORDINATES 

X1 

x2 

X3 

X4 

X5 

x6 

X7 

X8 

Vl 

y2 

y3 

y« 

y5 

y6 

y7 

y8 

zl 

Z2 

z3 

Z4 

Z5 

Z6 

z7 

ZB 

COMPUTED  COORDINATES 

X? 

X* 

xS 

X? 

Xn 

X£ 

x7 

XS 

y? 

y2 

y3 

y? 

y§ 

y? 

y? 

y£ 

z? 

z3 

z? 

z? 

z2 

Z7 

Z8 

CONNECTIVITY  CONSISTENCY  CHECK 

nnnnnnn 

mmnntiiin 

This  message  indicates  that  the  8  corner  grid  points  associated  with 
element  n  (from  the  CPOE3  card)  as  generated  using  one  of  the  various 
HPAT  directives  are  not  consistent  with  the  prior  definitions  of  the 
grid  point  coordinates.  Those  grid  points  may  have  been  input  with  a 
set  of  GRID  cards,  or  they  could  be  the  result  of  prior  hyperpatch 
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constructions.  Hie  tolerance  on  this  test  is  as  defined  by  the  ACGG 
field  on  the  PARAM  directive  (default  is  5  significant  figures) .  Check 
the  connectivity  definition  on  the  CPDE3  directive.  Check  the 
construction  operations,  including  the  construction  of  associated 
patches  and  lines.  Originally  undefined  grid  points  (for  which 
consistency  is  not  required)  are  represented  on  output  as  a  .999E-99. 
The  "connectivity  consistency"  check  is  a  check  on  the  degeneracy  of  the 
connectivity  and  geometry  definitions  and  should  be  identical.  A  final 
item  to  check  stems  from  the  fact  that  an  element  needs  at  least  one 
grid  point  to  "hang  on  to"  to  get  its  bearings.  An  element  constructed 
first,  or  alone  in  space,  is  "consistent"  with  8  undefined  grid  points 
in  its  original  parameterization. 

UNRECOGNIZABLE  BULK  DATA  CARD. 

The  first  field  on  the  card  listed  before  this  message  does  not  have  a 
valid  mnemonic.  Check  the  spelling  of  the  card.  Check  that  the  card 
prior  to  this  one  did  not  require  a  continuation  character. 

VALUE  OUT  OF  RANGE.  EATA  FIELD  NUMBER  n. 

INPUT  =  fn  RANGE  =  Rmin  Rmax 

The  card  image  that  follows  this  message  has  a  value  that  is  not  within 
a  range  of  reasonable  values  Rmin  <  fn  <  Rmax*  The  position  of  the 
offending  parameter  is  n  vrtiere  fields  1  to  10  are  on  the  first  card,  11 
to  20  on  the  first  continuation  card,  etc. 
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ERROR  IN  A  RASTUS  REQUEST  FROM  TOE  xxxxxx  MODULE. 

CALLING  INDEX  SET  TO  n. 

This  message  is  generally  preceded  and  followed  by  other  messages  vhich 
serve  to  further  identify  the  problem.  However,  the  xxxxxx  parameter  is 
helpful  in  deciphering  the  cause  of  the  problem.  There  are  two  basic 
forms  for  this  parameter:  FETCHx  and  STOREx  representing  a  random 
access  fetch  or  store  request  for  volune  x.  The  request  was  for  page  n 
of  that  volume.  The  table  below  summarizes  the  various  volunes  and 


their  contents. 
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Vo  lime 
A 
B 
C 
D 
E 
F 
G 
H 
I 
J 
K 
L 
M 
N 
0 
P 
0 
R 
S 
T 
U 
V 
W 


Contents 

Auxiliary  storage 
Data  patches 
Material  properties 
Data  lines 
Element  matrices 
Function  tables 
Element  connectivity 
Geometric  hyperpatches 
Input  MAT  data 
Jacobi ans 

ED  tensor  for  elements 
Geometric  lines 
Matrices 

Mesh  point  numbers  for  elements 
Output  requests 
Geometric  patches 
Load  tables  by  elements 

R  table  of  imposed  displacements  for  mesh  points 

Rigid  body  transformation 

Load  vectors  by  elements 

Element  displacements,  strains,  stresses 

Not  active 

Element  matrices  F  set  and  partitioning 
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Volume  Contents 

X  Original  C  and  alpha  matrices 

Y  Mesh  point  vs.  F  table  ID's 

Z  Temperature  hyperpatches 

1  One-component  data  hyperpatches 

2  Element  tables  and  vectors 

3  Archive  element  matrices,  E  frame 

4  Packed  element  constraints  matrices 

5  Element  mesh  point  constraint  flags 

6  U*  Master  cartesian  displacements 

For  example,  a  message  concerning  the  FETCHP  module  with  an  index  of  7 
is  referring  to  a  problem  fetching  geometric  patch  7.  Therefore  check 
that  patch  7  was  properly  generated. 

THE  RESULTANT  GEOMETRY  FOR  HYPERPATCH  n  IS  NOT  GEOMETRICALLY  REASONABLE. 
The  construction  of  hyperpatch  n  resulted  in  local  regions  at  which  the 
Jacobians  at  the  Gaussian  points  were  negative.  This  implies  a  negative 
local  volume.  Check  the  connectivity  and  the  construction  operations. 
An  element  such  as  this  can  still  be  plotted,  and  this  can  serve  as  a 
clue  to  the  problem.  As  an  example  of  this  situation,  consider  a  ruled 
hyperpatch  in  vrtiich  one  of  the  patches  is  reversed  (physically  or  by 
connectivity) .  Then  the  lines  will  tend  to  cross  near  the  center  of 

the  element. 
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9.3  Debugging  Data 

Additional  information  can  tie  printed  to  aid  in  the  debugging  of  a 
model  through  the  use  of  the  PARAM  DEBUG  bulk  data  card.  The  user 
inputs  an  integer  value  for  the  DEBUG  parameter,  and  bits  of  the 
resultant  word  activate  the  debugging  options  listed  in  Table  9-1. 

A  DEBUG  value  of  138,  for  example,  would  activate  three  options 
since  138  =  128  +8+2.  Currently  only  two  bits  in  this  table  are 
active. 
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Table  9-1 

DEBUG  OPTIONS  AVAILABLE 


Base 

Value 

Bit 

Position 

Result 

1 

1 

Not  active 

2 

2 

Imposed  displacement  flags  and  YS  vector 

4 

3 

Not  active 

8 

4 

Not  active 

16 

5 

Not  active 

32 

6 

Not  active 

64 

7 

Not  active 

128 

8 

Imposed  displacement  transformations  (row  sort) 
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The  imposed  displacement  flags  printout  (bit  2)  gives  for  each  mesh 
point  the  code  numbers  representing  the  degrees  of  freedom.  Only  three 
distinct  values  exist  in  this  table  in  PATCHES-III:  50  for  constrained 
to  0.,  496  for  unconstrained,  and  1074  for  constrained  to  a  nonzero 
value.  The  YS  table  gives  the  magnitude  of  the  imposed  displacement  for 
each  component  of  each  mesh  point  consistent  with  USET. 

The  imposed  displacement  transformations  printout  (bit  8)  defines 
for  each  mesh  point  having  an  imposed  displacement  the  3X3 
transformation  matrix  (in  row  sort)  that  maps  displacement  components  in 
the  constraint  frame  to  components  in  the  reference  frame. 
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9.4  User  Information  File 


The  final  output  from  any  PATCHES-III  execution  is  the  User 
Information  File.  This  file  is  a  chronological  listing  of  the  major 
events  during  the  execution  of  the  program  and  the  resources  required. 
The  file  has  the  appearance  of  a  CDC  day  file  and  can  be  used  to 
estimate  CPU  costs  per  element  on  any  given  computer.  For  each  of  the 
categories  of  output,  the  following  definitions  apply. 


PROGRAM  REGION 
SUBREGION 

SSTIME 


DELTA 


The  major  region,  usually  the  link  name. 

The  subregion  name  or  start  or  end  for  a  program 
reg ion . 

The  number  of  system  seconds  (or  CRJ  seconds  on  CDC) 
up  to  this  point  in  the  job. 

NOTE:  A  job  will  begin  with  SSTIME  greater  than  0 
if  other  operations  are  performed  prior  to 
initiating  the  PATCHES  execution. 

The  difference  in  SSTIME  between  this  and  the  prior 
event. 


FIELD  LENGTH 

RASTUS  REQUESTS 

TOTAL  W0RD6 
TRANSFERRED 


Octal  total  field  length  at  this  time. 

Total  number  of  RASTUS  random  access  requests. 

Toual  number  of  words  transferred  (written  and  read) 
using  RASTUS. 


DELTA  WORDS 
TRANSFERRED 


SUB  INDEX  SORTS 


PERCENT  SORTS 


The  difference  in  TOTAL  WORDS  TRANSFERRED  between 
this  and  the  prior  event. 

The  number  of  RASTUS  subindex  sorts,  an  overhead 
requirement  for  the  adaptive  subindex  buffer 
technique. 

The  percent  of  RASTUS  REQUESTS  requiring  SUBINDEX 
SORTS.  Generally  under  3%. 
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